1

我对 R 真的很陌生,并且在浏览 SO 和其他 R 资源的两个晚上一直在努力解决这个问题......我有一个数据集,其中有两个矩阵,一个用于股市观察(股票收益)和一个用于权重(指数市值的一部分)。在这两种情况下,列是该股票代码的数据,行是该月的数据。

我希望使用 density() 在第二步中使用 Kullback Leibler 散度来比较行(每月收益)。我的问题是我无法找出一种语法来移动使用的权重行 w[i,](因为我一直在使用 apply() / lapply(),这可能是问题所在)然后存储 density() 输出矩阵整齐地列在一个列表中。我的结果只是最终成为第一个重量行而不是整体。

我也尝试过使用:

mylist <- list()
for (x in 1:n) {mylist <- density(R[x,]], 1, weights = (W[x,])}

由于某种原因没有输出密度()的概率/ y值输出......

我将不胜感激任何帮助!

4

0 回答 0