Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- #!/usr/bin/octave -qf
- filename = input("Please enter the input file:", "s")
- A = load ("-ascii", filename)
- ch = menu ("Please choose from one of the following:", "linear", "log-x", "log-y", "log-log")
- switch (ch)
- case 1
- x = A(:, 1);
- y = A(:, 2);
- case 2
- x = log (A(:, 1));
- y = A(:, 2);
- case 3
- x = A(:, 1);
- y = log (A(:, 2));
- case 4
- x = log (A(:, 1));
- y = log (A(:, 2));
- endswitch
- Nobs = rows(A);
- X = [ones(Nobs,1) x];
- [beta sigma] = ols(y,X)
- yest = X * beta;
- p = length(beta);
- varb = sigma * ((X'*X)\eye(p));
- se = sqrt(diag(varb))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement