这个真的让我很难过。我找不到这行代码有什么问题,我在函数中使用它,它使整个函数失败。
这是“矩阵”类的向量,result.1
:
SPY TLT VGK EEM BIV RWX IYR EWJ DBC GLD PCY LQD TIP
-0.08 0.09 -0.07 -0.07 -0.09 -0.07 -0.07 -0.07 -0.066 -0.04 -0.08 -0.08 -0.08
我通过运行以下代码来计算此向量的平均值:
mean.momentum.1 <- mean(result.1)
结果是一个值为-0.062309的数字。我想得到大于平均值的列的名称,所以我运行这段代码:
names(result.1[,(result.1 > mean.momentum.1)])
输出正如我所料,返回一个字符向量“TLT”,“GLD”。
问题在于:当我使用第二个几乎相同的矩阵result.2
执行此操作时,每次得到NULL
结果时,我应该得到结果“TLT”。
这是矩阵result.2
:
SPY TLT VGK EEM BIV RWX IYR EWJ DBC GLD TIP PCY LQD
-0.08 0.15 -0.07 -0.06 -0.054 -0.07 -0.07 -0.07 -0.06 -0.06 -0.07 -0.07 -0.06
我使用与上面相同的方法计算平均值(-0.05298)并将其命名为mean.momentum.2
然后,我运行这一行:
names(result.2[,(result.2 > mean.momentum.2)])
每当我期待“TLT”时,我都会回复NULL
。我的子集化方式出了什么问题?如果我运行该行:
result.2 > mean.momentum.2
我得到一个逻辑向量,其中TLT是TRUE
,但是当我尝试用这个方法进行子集时它不起作用。它在第一个实例中完美地工作,但从不在第二个实例中工作,并且因为我得到NULL
,所以我的整个函数都失败了。
谢谢...
答案 0 :(得分:2)
在引用矩阵时使用colnames()
和rownames()
而不是names()
:
mean.momentum.1 <- mean(result.1)
colnames(result.1)[which(result.1 > mean.momentum.1)]
## [1] "TLT" "GLD"
另请注意上面colnames()
的语法。您使用colnames(result1)
返回向量which()
的子集,以选择result.1
大于mean.momentum.1
的向量元素。
这适用于result.2
:
mean.momentum.2 <- mean(result.2)
colnames(result.2)[which(result.2 > mean.momentum.2)]
## [1] "TLT"