我正在运行以下代码并返回以下错误 - 任何解决方案都将非常感谢。
## UQdata = Positive count abundance only data (i.e. positive values only)
aa2<-gam(UQdata~s(MudUQ,bs="ps", k=15) ,family=Gamma(link=log),data=Antho)
xmin <- ceiling(min(Antho$MudUQ[Antho$Bin==1]))
xmax <- floor(max(Antho$MudUQ[Antho$Bin==1]))
Mudnew <- seq(from=xmin, to=xmax, by=0.1)
**Error in if (del == 0 && to == 0) return(to) :
missing value where TRUE/FALSE needed**
pred.dat <- data.frame(Mudnew)
names(pred.dat) <- "MudUQ"
pred.aa2 <- data.frame(predict.gam(aa2, pred.dat, se.fit=TRUE, type="response"))
pred.aa2.comb <- data.frame(pred.dat, pred.aa2)
names(pred.aa2.comb)
plot(fit ~ MudUQ, data=pred.aa2.comb, type="l", lwd=2, col=1, ylab="Density per 0.0132 m2", xlab="Mud content (%)")
答案 0 :(得分:1)
也许:
xmin <- ceiling( min( Antho$MudUQ[Antho$Bin==1] , na.rm=TRUE) )
xmax <- floor( max( Antho$MudUQ[Antho$Bin==1] , na.rm=TRUE) )
我认为seq
函数会抛出错误。