r - 根据R中的间隔(子集函数)选择矩阵

Question

我想基于X2列值对一个矩阵进行子集化。我试了一下:

on <- subset(mat.num, X2 <= -3)
un <- subset(mat.num, X2 >= -1.50000 & X2 <= -0.3599999)
dn <- subset(mat.num, X2 >= -0.3599998 & X2 <= 0.5)

但是我得到这个错误:

Error in subset.matrix(mat.num, X2 <= -3) : object 'X2' not found.

ps:我有一个 X2列:
mat.num头:

     T_EBV     X2
[1,] 0.09 -0.00777840
[2,] 0.26  0.03600431
[3,] 0.20 -0.06191900
[4,] 0.25  0.13423752
[5,] 0.42  0.06354759
[6,] -0.20  0.06303164

Answer 1

矩阵方法没有像使用数据框那样引用列名。您可能想要:

subset(mat.num, mat.num[,2] <= -3)

如果查看 subset.matrix的代码，您会发现它没有在任何特殊环境中评估子集条件:

if (missing(subset)) 
        subset <- TRUE
else if (!is.logical(subset)) 
        stop("'subset' must be logical")
x[subset & !is.na(subset), vars, drop = drop]

与使用 subset.data.frame和 eval的 substitute相反。