回歸投資組合超額收益

Question

我的數據：

dput(head(mydata))
structure(list(DATE = structure(c(-315619200, -312940800, -310435200, 
-307756800, -305164800, -302486400), tzone = "UTC", class = c("POSIXct", 
"POSIXt")), RF = c(0.33, 0.29, 0.35, 0.19, 0.27, 0.24), RMRF = c(-6.99, 
0.99, -1.46, -1.7, 3.08, 2.09), SMB = c(2.13, 0.71, -0.65, 0.32, 
1.42, -0.24), UMD = c(-3.28, 3.59, 1.85, 2.6, 4.77, 1.03), HML = c(2.65, 
-2.15, -2.69, -2.22, -3.83, -0.3), JANDUM = c(1, 0, 0, 0, 0, 
0), R4 = c(-4.57, 1.5, -2.83, -1.98, 3.54, 2.15)), row.names = c(NA, 
-6L), class = c("tbl_df", "tbl", "data.frame"))

所以，我的數據包含：

R4 是投資組合的百分比回報率，RF 是無風險商品的回報率（無風險利率），RMRF 是投資組合 Market Portfolio 的超額回報，SMB、UMD 和 HML 是 3 個因子，JANDUM 是虛擬變量一月（一月虛擬）。

數據是從 1/1960 到 12/2003 的每月頻率（總共有 528 個觀測值）。

那就是我試圖建立並且我正在努力將投資組合超額收益（R4-RF）回歸到一個常數和所有其他變量（RMRF、SMB、UMD、HML 和 JANDUM）。

我怎樣才能做到這一點？

Answer 1

也許這會讓你開始？

mydata$PER <-mydata$R4 - mydata$RF 
mydata$JANDUM <- as.factor(mydata$JANDUM)
model <- lm(PER ~ DATE + RMRF + SMB + UMD + HML + JANDUM, data = mydata)
summary(model)