如何解决矩阵的t检验没有足够的“ x”观察结果
我正在尝试对矩阵进行t检验,并将结果输出到最后一列。我在df中创建了示例数据,并在脚本中复制了运行df的错误。我收到一个错误,即没有足够的“ x”观测值。关于代码中有什么错误的任何想法?
df <- data.frame(
Sample_S_0_rep1 = c(0,2478.55,3.91,0.69,0.63,0),Sample_S_0_rep2 = c(0,789.63,2.98,0.15,Sample_S_0_rep3 = c(0,802.76,2.99,0.02),Sample_S_24_rep1 = c(0,993.88,3.71,0.35,Sample_S_24_rep2 = c(0,896.23,3.48,0.09,Sample_S_24_rep3 = c(0,947.83,4.02,0)
)
group1 <- df[i,grep('^Sample_S_0',names(df))]
group2 <- df[i,grep('^Sample_S_24',names(df))]
group1 <- log2(as.matrix(group1)+1)
group2 <- log2(as.matrix(group2)+1)
Result <- apply(df,1,function(df){t.test(df[group1],df[group2])$p.value})
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。