Question

我正在尝试创建一个具有p值的数据框，并估计将一个基因与许多不同的表达标记进行比较。当我仅在一个表达式上使用cor.test时，我的cor.test可以工作，但是当我尝试循环运行它时，它就会中断，并给我以下错误：“'x'和'y'必须具有相同的长度”。我想知道如何使此循环正常工作并构建数据框架。

下面是我正在循环中运行的内容以及该循环的代码。

M3 <- ads$mean

Expression <- c("Exp1","Exp2","Exp3")

for (i in seq_along(Expression))
{
  corr<-cor.test(M3, Expression[i], method = "pearson")
  cor_df<-data.frame(Expression = Expression[i],pvalue = corr$p.value, 
 cor = corr$estimate)

 }

Answer 1

根据您的评论，如果Exp1，Exp2和Exp3是数据框（df）中的列，则可以使用以下内容：

corr <- cor.test(M3, df[ ,Expression[i]], method = "pearson")

R中的相关性测试循环

1 个答案: