我看了一下计算普通最小二乘回归的函数ols:
ols(Y, X, [intercept=true], [mode=0])
这里因变量Y必须是一个是向量,不能是矩阵,所以我来请假一下应该怎么实现?
可以用高阶函数cross(func, X, [Y])以列为单位遍历计算:
pcross(ols, X, matrixY)
如果x是向量,可以如下实现:
each(ols{,x}, matrixY)