Advertisement
Guest User

Untitled

a guest
Apr 27th, 2017
92
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.58 KB | None | 0 0
  1. #!/usr/bin/octave -qf
  2. filename = input("Please enter the input file:", "s")
  3. A = load ("-ascii", filename)
  4. ch = menu ("Please choose from one of the following:", "linear", "log-x", "log-y", "log-log")
  5. switch (ch)
  6. case 1
  7. x = A(:, 1);
  8. y = A(:, 2);
  9. case 2
  10. x = log (A(:, 1));
  11. y = A(:, 2);
  12. case 3
  13. x = A(:, 1);
  14. y = log (A(:, 2));
  15. case 4
  16. x = log (A(:, 1));
  17. y = log (A(:, 2));
  18. endswitch
  19. Nobs = rows(A);
  20. X = [ones(Nobs,1) x];
  21.  
  22. [beta sigma] = ols(y,X)
  23. yest = X * beta;
  24.  
  25. p = length(beta);
  26. varb = sigma * ((X'*X)\eye(p));
  27. se = sqrt(diag(varb))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement