待解决问题:
做代谢组学分析,一共有147种代谢物,需要每个单独和其他选定的协变量一起放入线性回归模型中,依次替换147种代谢物,完成147次回归分析。
实现步骤:
R语言
1、 数据预处理:excel和spss
2、 for循环
3、 multinom() 线性回归
4、 导出结果分析
步骤1、2、4参见我之前写的《批量线性回归并整理导出结果到word文档》,这里只讲批量multinom() 线性回归的核心步骤。
具体操作:
3. multinom() 线性回归
timo <- data.frame() #先创建一个空数据框,否则后面<- 赋值会提示找不到timo
for (q in c(222:368)){
result <- multinom(tdata$Y <