p17 oxtut2e shortcut.ox

时间序列分析中著名的OxMetrics软件包

#include <oxstd.h>

MyOls(const mY, const mX)
{
//decl b;
//return olsc(mY, mX, &b);
return (mX'mX)^-1 * (mX'mY);
} 

main()
{ decl b;
// mY, mX both just random
b = MyOls( rann(10, 1), ranu(10, 2) );
// b now holds the result
print("in main(): b=", b);
}