r - 试图通过单索引模型进行预测时出现两个R错误消息

Question

我正在尝试使用单一指数模型来估算3只股票和TSX的alpha beta和sigma2_ei。这是我当前的R代码

导入数据

setwd("~/Desktop/R Data Sets") Ass2Data1 <- read.csv("~/Desktop/R Data     Sets/Ass2DataSheet1.csv") View(Ass2Data1)

将数据转换为矩阵形式

b <- as.matrix(Ass2Data1)

生成初始 vector 和矩阵

x <- rep(0,60) 
xx <- matrix(x, ncol=4, nrow=3)
stock <- rep(0,3) 
alpha <- rep(0,3) 
beta <- rep(0,3) 
mse <- rep(0,3) 
Rbar <- rep(0,3) 
Ratio <- rep(0,3)

col1 <- rep(0,3) 
col2 <- rep(0,3) 
col3 <- rep(0,3) 
col4 <- rep(0,3) 
col5 <- rep(0,3)

使指数上的每只股票回归并记录结果

for(i in 1:3){
alpha[i] <- lm(data=Ass2Data1,formula=Ass2Data1[,1] ~ Ass2Data1[,4]$coefficients[1])
beta[i] <- lm(data=Ass2Data1,formula=Ass2Data1[,2] ~ Ass2Data1[,4]$coefficients[2])
Rbar[i] <- alpha[i]+beta[i]*mean(b[,4])
mse[i] <-sum(lm(data=Ass2Data1,formula=Ass2Data1[,i] ~ Ass2Data1[,4])$residuals2)/(nrow(b)-2)
Ratio[i] <- (Rbar[i]/beta[i])
stock[i] <- i }

创建表

xx <- (cbind(stock,alpha,beta,Rbar,mse,Ratio))

但是，我不断收到以下错误消息:

Ass2Data1 [，4] $ coefficients中的错误:$运算符对原子 vector 无效

和

beta [i] * mean(b [，4])中的错误:二进制运算符的非数字参数

如果有人能就我在这里做错的事情为我指明正确的方向，将不胜感激。

Answer 1

检查?lm，更改结束)

alpha[i] <- lm(data=Ass2Data1,formula=Ass2Data1[,1] ~ Ass2Data1[,4])$coefficients[1]
beta[i] <- lm(data=Ass2Data1,formula=Ass2Data1[,2] ~ Ass2Data1[,4])$coefficients[2]

更好的方法是通过 coef获得系数:

alpha[i] <- coef(lm(data=Ass2Data1,formula=Ass2Data1[,1] ~ Ass2Data1[,4]))[1]
 beta[i] <- coef(lm(data=Ass2Data1,formula=Ass2Data1[,2] ~ Ass2Data1[,4]))[2]