# Applied Exercises ## 1. Writing a custom "regress" function In this exercise you will write your custom function which accepts a vector of dependent variables and a matrix of regressors and computes the OLS estimator as well as its standard error. {% tabs %} {% tab title="Exercise" %} Write a function which accepts a vector `y` and a matrix `X` (without the column of ones) as inputs and which returns a $$2 \times N$$ matrix where the first row carries the estimated regression coefficients and the second row carries their standard errors. Check your function works appropriately by comparing to MATLAB's `regress`function. **Extension 1 - Constant in the regression:** Extend the function to accept a boolean argument `constant`which determines whether a constant should be included in the regression. **Extension 2** - **Using structure array objects to return estimated values):** Read the [MATLAB reference page on structure arrays](https://www.mathworks.com/help/matlab/ref/struct.html) to learn about `struct` objects. Modify your OLS function such that it returns a struct object `estimates` with sub-objects `estimates.params` and `estimates.se` which contain the parameter estimates and standard errors, respectively. {% endtab %} {% tab title="Getting started" %} * Create a new folder for this exercise and copy the *olsdata.mat* file that you downloaded from this page into it. * In MATLAB, create a new script and save it into the same folder naming it e.g. `exercise01.m`. * Start your script with the following commands ``` clear all; close all; clc; load olsdata.mat ``` This will clear everything and load the data into your workspace. * Inspect the workspace window. You should see a matrix `X` and a vector `y`. * The goal of this exercise is to estimate `beta_true` by ordinary least squares (OLS) and to compute its standard errors using the formula above. * Like in exercise one of unit one, code up the ols estimator using matrix algebra. ``` beta_hat = ...formula goes here... ``` * Then compute the estimate $$\hat{\sigma}^2$$ using the formula from the theory section ``` sigma2_hat = ...formula goes here... ``` * Finally, estimate the variance-covariance matrix of the estimator $$\hat{Var}(\hat{\beta} \lvert X)$$ using the formula given above. Note that the expression should involve your estimate `sigma2_hat`. ``` V_beta_hat = ...formula involving sigma2_hat... ``` * Finally, extract the estimated variances of $$\hat{\beta}$$ from the diagonal of the variance-covariance matrix and take the square root to get standard errors. ``` se_beta = ...formula involving diagonal of V_beta_hat... ``` * Now you have two column vectors, `beta` which is $$3 \times 1$$ and `se_beta` which is $$3\times 1$$. Think about how to concatenate them into a $$2 \times 3$$ matrix where the first row carries `beta` and the second row carries `se_beta`. * After successfully testing your script, move all code after loading the data into a new file where you declare a function that accepts `y` and `X` as inputs and returns the concatenated matrix. Save this function as `run_ols.m` in your main folder. * Modify your script so that it looks like this: {% code title="exercise01.m" %} ``` clear all; close all; clc; load olsdata.mat ols_result = run_ols(y, X); ols_result ``` {% endcode %} * Run your script and verify that your OLS function is being called and that the correct result is delivered. {% endtab %} {% endtabs %} {% file src="" %} Datasets from 'Crashcourse' section {% endfile %} ### Theory Let $$y$$ be a $$T \times 1$$ vector of data on the dependent variable and let $$X$$ be a $$T \times N$$ matrix with data on the regressors (where the first column is a vector of ones). The OLS estimator of the regression coefficients is defined as. $$\hat{\beta} = (X'X)^{-1} (X'y)$$ Under the assumption of spherical errors, and estimate of the variance-covariance matrix of the estimator can be obtained as $$\hat{Var}(\hat{\beta} \vert X) = \hat{\sigma}^2 ;(X'X)^{-1}$$ where $$\hat{\sigma}^2 = \frac{(y-X\hat{\beta})'(y-X\hat{\beta})}{T-N}$$ The standard errors of the elements of $$\hat{\beta}$$ are defined as the square root of the diagonal elements of $$\hat{Var}(\hat{\beta}\vert X)$$.
## 2. Simulating capital evolution in the canonical OLG model In this exercise you will simulate a capital evolution in a simple canonical OLG model. {% tabs %} {% tab title="Exercise" %} Write a script that calculates a vector that contains the evolution of capital for $$T=30$$ periods when initial capital per worker is $$k\_0=0.1$$ for the following parameter values. | **Parameter** | **Value** | | ------------- | --------- | | $$A$$ | 4 | | $$\alpha$$ | 0.4 | | $$\beta$$ | 0.9 | | $$n$$ | 0.02 | **Tip:** Create a vector of zeros of appropriate size, then loop through the elements of the vector and fill in the values of capital per worker. When you have finished you can plot the resulting vector with the command `plot(k)` where `k` is the vector holding the evolution of capital. How many periods does it take the economy to reach its steady state (i.e. capital per worker does not change anymore)? Now try out to simulate the economy for different parameter values. Does a steady state exist and how long does it take to get there? {% endtab %} {% tab title="Getting started" %} * Create a new folder for this exercise and create a new script in MATLAB which you save into the same folder naming it e.g. 'exercise02.m'. * Start your script with the following commands ``` clear all; close all; clc; ``` This will clear everything. * We would like to store the values of $$k\_t$$ for $$t=1, ..., T$$ in a vector. Initialise the vector at the beginning of your code by creating a vector of `NaN` or `zero` values. ``` k = nan(...appropriate dimensions here...) ``` * Think about how you can use a loop to calculate the values of $$k$$ iteratively. How does this loop look like? Where does it start and where does it end. How do you save a single value that you computed for $$k\_{t+1}$$ in the vector you created in the step before? {% endtab %} {% tab title="Check" %} If you have implemented the exercise correctly, your plot of the evolution of the capital stock should look like this. ![Evolution of the capital stock in the simulated OLG economy](https://1418968082-files.gitbook.io/~/files/v0/b/gitbook-legacy-files/o/assets%2F-Lp2LO4uQaxs9KN1kWAs%2F-LpyU4KStVTTnR1hgxJG%2F-LpyULXCdID7jql2LCR2%2Folg-capital-stock.png?alt=media\&token=f3015471-1a1d-453d-b936-8f82633ff923) {% endtab %} {% endtabs %} ### Theory In the canonical overlapping-generations (OLG) model, the evolution of the capital stock per worker is described by the following law of motion. $$k\_{t+1} = \frac{1}{1+n} \frac{\beta}{1+\beta} ; (1-\alpha) ; A ; k\_t^\alpha$$ where $$k$$ is the capital per worker, $$n$$ is the population growth rate, $$\beta$$ is the discount factor, $$\alpha$$ is the capital share in the production function and $$A$$ is the TFP parameter from the production function. ## 3. Minimization of a non-differentiable function by bisection In this exercise you will implement a simple algorithm to find the minimum of a non-differentiable function by bisection minimization. {% tabs %} {% tab title="Exercise" %} Code up the algorithm that is described in the theory section below to find the minimum of $$f(x)$$ on the interval $$\[0, 5]$$. **Tips:** Use a while loop which iterates on the algorithm steps 1-4 until convergence is achieved. Use variables to keep track of the old and new candidate values for $$x^{min}$$. {% endtab %} {% endtabs %} ### Theory In this exercise, you will find the minimum of the following non-differentiable function on the interval $$\[a, b]$$. $$f(x) = \lvert x - 1 \rvert + (x-2)^2$$ Here is a plot of this function with the respective minimum ![Function to be minimized](https://1418968082-files.gitbook.io/~/files/v0/b/gitbook-legacy-files/o/assets%2F-Lp2LO4uQaxs9KN1kWAs%2F-LpyUWCvOG99dBOj8FHF%2F-LpyV9u38Afg1mOKcBPG%2Fbisection-min-function.gif?alt=media\&token=8961cebc-7155-4f77-b746-7ff5aef4a84b) To find the minimum of $$f(x)$$ on the interval, use the following algorithm: 1. Pick the point in the middle of the interval $$x\_M$$. This splits the interval into two intervals, the left interval $$\[a, x\_M]$$ and the right interval $$(x\_M,b]$$ 2. In each half, pick as a candidate for a minimum $$x\_L$$ and $$x\_R$$ as the middle point of the respective interval 3. Evaluate $$f(x\_L)$$ and $$f(x\_R)$$ 4. Now two cases are possible * (a) If $$f(x\_L)>f(x\_R)$$ i.e. the minimizer is in the right interval. If this is true, set the new interval equal to the right interval. * (b) If $$f(x\_L)<=f(x\_R)$$, the minimizer is in the left interval. If this is true, set the new interval equal to the left interval. 5. Using the new interval $$\[a^*, b^*]$$, repeat steps 1-4 until the absolute difference between the new and old candidate is smaller than `1e-5`. --- # Agent Instructions: Querying This Documentation If you need additional information that is not directly available in this page, you can query the documentation dynamically by asking a question. Perform an HTTP GET request on the current page URL with the `ask` query parameter: ``` GET https://coding-courses.gitbook.io/matlab/matlab-programming/exercises.md?ask= ``` The question should be specific, self-contained, and written in natural language. The response will contain a direct answer to the question and relevant excerpts and sources from the documentation. Use this mechanism when the answer is not explicitly present in the current page, you need clarification or additional context, or you want to retrieve related documentation sections.