Корреляционный тест петли в R

Question

Я пытаюсь создать фрейм данных со значениями p и оценками, которые сравнивают один ген со многими различными маркерами экспрессии.Мой тест 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")