问题描述
我有一个包含 20 名参加交叉研究的受试者的数据集。 我们在 3 次不同的干预(20 x3 = 60 个样本)后测量了 3 次相同的代谢物(n=288)。
我想创建一个循环来对每个代谢物的单向重复测量方差分析执行事后分析。
我尝试了以下方法;然而,每一个都给了我同样的错误。 我附上了脚本和相关的错误信息。
谢谢
for (i in (3:ncol(data))){
#variable = data[,i]
variable = print(colnames(data[i]))
pwc <- data %>% pairwise_t_test(variable ~ Diet,paired = TRUE,p.adjust.method = "holm") #default
outFile <- paste(variable,"_one-way_repeated_ANOVA_post-hoc.txt",sep="")
#write.table(pwc,file=outFile,row.names =FALSE,col.names=T,quote=F,sep="\t")
}
##Error message:
# Error: Can't extract columns that don't exist.
# Column `variable` doesn't exist.
#### OR ####
ColNames <- colnames(data[3:ncol(data)])
for (i in ColNames){
variable = print(i)
pwc <- data %>% pairwise_t_test(variable ~ Diet,sep="\t")
}
##Error message:
# Error: Can't extract columns that don't exist.
# Column `variable` doesn't exist.
#### OR ####
ColNames <- colnames(data[3:ncol(data)])
for (i in ColNames){
pwc <- data %>% pairwise_t_test(i ~ Diet,p.adjust.method = "holm") #default
outFile <- paste(i,sep="\t")
}
##Error message:
# Error: Can't extract columns that don't exist.
# Column `i` doesn't exist.
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)