我目前正在为学校的一个项目运行方差分析,该项目有大量可能的运行(1400左右),但其中一些无法在R中运行方差分析。我写了一个脚本来运行所有的方差分析,但其中一些将不会运行,并且路由文件在contrasts<-(*tmp*,value = "contr.treatment")中给我错误:对比只能应用于具有2个或更多级别调用的因子:aov...model.matrix -> model.matrix.default ->对比<-执行暂停
有没有办法编写脚本,让R在出现错误的情况下继续执行脚本?
我的整个脚本,除了文件加载,附加,创建接收器,库加载等是...
ss107927468.model<-aov(Race.5~ss107927468, data=snp1)
summary(ss107927468.model)任何帮助都将不胜感激。
发布于 2010-12-05 23:23:44
请参阅函数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)
}发布于 2010-12-06 05:30:01
我在plyr包中使用failwith。您可以将其与llply结合使用,并将您的函数封装在它周围。
https://stackoverflow.com/questions/4359310
复制相似问题