我目前正在为学校的一个项目运行方差分析,该项目有大量可能的运行(1400 次左右),但其中一些无法在 R 中运行方差分析。我编写了一个脚本来运行所有方差分析,但其中一些它们不会运行,并且 Rout 文件给了我
错误于contrasts<-
(*tmp*
, 值 = "contr.treatment") :
对比只能应用于具有 2 个或更多级别的因素
调用:aov ... model.matrix -> model.matrix.default ->对比
有没有办法编写脚本,让 R 在出现错误的情况下继续执行脚本?
我的整个脚本,除了文件加载、附加、创建接收器、库加载等之外......
ss107927468.model<-aov(Race.5~ss107927468, data=snp1)
summary(ss107927468.model)
任何帮助,将不胜感激。
查看功能try()
这是帮助页面(?try
)。您将 R 表达式包装在try()
调用,如果成功,则生成的对象包含(在本例中)拟合的模型。如果失败,则使用具有类的对象"try-error"
被返回。这使您可以轻松检查哪些模型有效,哪些无效。
您可以进行测试来决定是打印模型摘要还是仅打印失败消息,例如:
ss107927468.model <- try(aov(Race.5~ss107927468, data=snp1))
if(isTRUE(all.equal(class(ss107927468.model), "try-error"))) {
writeLines("Model failed")
} else {
summary(ss107927468.model)
}
本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)