我有一个(kxn)矩阵。我最初设法将第 1 列与其他列进行线性回归(使用 lm 函数),并仅提取系数。
fore.choose <- matrix(0, 1, NCOL(assets))
for(i in seq(1, NCOL(assets), 1))
{
abc <- lm(assets[,1]~assets[,i])$coefficients
fore.choose[1,i] <- abc[2:length(abc)]
}
系数被放置在fore.choose
矩阵中。
我现在需要做的是将第 2 列与其他列进行线性回归,然后将第 3 列等以此类推,并仅提取系数。
输出将是 OLS 单变量系数的方阵。有点类似于相关矩阵,但它是我感兴趣的 beta 系数。
fore.choose <- matrix(0, 1, NCOL(assets))
最初需要成为
fore.choose <- matrix(0, NCOL(assets), NCOL(assets))