1 列与样本大小不等的所有其他列之间的协方差

Question

我正在尝试为每个公司找到 beta 值（协方差（x_i，y_i）/sd（y_i）），但由于我的两个数据框的大小不同，因此出现“尺寸不兼容”错误。

第一个数据框：（从 2009 年到 2020 年）

bist <- structure(list(Date = structure(c(14245, 14252, 14259, 14266, 
14273, 14280, 14287, 14294, 14301, 14308), class = "Date"), Adj_Close = c(27892.699219, 
25055.699219, 24736.099609, 26705.599609, 26141.099609, 25398.400391, 
23844.699219, 24152.300781, 23940.699219, 23282.599609)), row.names = c(NA, 
10L), class = "data.frame")

第二个数据框：

weekly_format <- structure(list(tarih = structure(c(1349308800, 1337212800, 1342051200, 
1345075200, 1340236800, 1345680000, 1337817600, 1347494400, 1344470400, 
1338422400), tzone = "UTC", class = c("POSIXct", "POSIXt")), 
    stock_code = c("ACIBD", "ACIBD", "ACIBD", "ACIBD", "ACIBD", 
    "ACIBD", "ACIBD", "ACIBD", "ACIBD", "ACIBD"), adj_close_price = c(24.85, 
    23.8, 24.55, 24.5, 24.3, 24.85, 23.75, 24.75, 24.85, 23.7
    )), row.names = c(NA, -10L), class = c("tbl_df", "tbl", "data.frame"
))

我所做的：

bist <- bist %>% 
    mutate(bist_return = Adj_Close/lag(Adj_Close)-1)

weekly_format_return <-  weekly_format %>% #weekly return
    group_by(stock_code) %>% 
    mutate(week_return = adj_close_price/lag(adj_close_price)-1) %>% 
    filter(n() >= 450) #i need at least 450 weeks of data

> head(weekly_format_return)
# A tibble: 6 x 5
# Groups:   stock_code [1]
 tarih      stock_code   adj_close_price     mon_return
   <date>     <chr>         <dbl>             <dbl>
1 2020-03-31 ACIBD           7.1               NA     
2 2020-03-24 ACIBD           6.9             -0.0282
3 2020-03-17 ACIBD           6.2             -0.101 
4 2020-03-10 ACIBD           7.96             0.284 
5 2020-03-03 ACIBD           8.88             0.116 
6 2020-02-25 ACIBD           9.26             0.0428

我试图找到指数值 (bist$bist_return) 和不同公司的股票收益 (weekly_format_return$week_return) 之间的协方差，然后将其除以 bist$bist_return 的标准值，但由于某些公司的值只有一段时间我找不到协方差。

covariance <- cov(bist$bist_return, weekly_format_return$week_return)/sd(bist$bist_return, na.rm = TRUE)

给出“不兼容的尺寸”错误。

也许我应该使用“use”参数来处理缺失值（它们实际上并没有缺失，它们只是不存在）但我不知道该怎么做。

感谢并感谢您的回复。