一个10*10的矩阵matrixname,
把每一行的进行1:6和7:10,做t.test。
在此之前,对每行做了f-test,
T_equ <- matrixname[F_equ, ]
T_unequ <- matrixname[-F_equ, ] #我看了T_unequ ,没问题
然后把方差相等的行拿出来做var.equal=TRUE的t.tes
t,没问题:
T_equ_te <- apply(T_equ, 1, function(x) t.test(x[1:6],x[7:10], var.equal=TRUE)[[3]]) #索引3就是想把每个P-value提出来给变量)
再把方差不相等的行拿出来做方差不等的t.test,报错了
T_unequ_te <- apply(T_unequ, 1, function(x) t.test(x[1:6],x[7:10])[[3]])#不等就啥也不用加
Error in t.test.default(x[1:6], x[7:10]) :
data are essentially constant