**if(del==0&;&;to==0)返回(to)时出错:缺少需要TRUE/FALSE的值**

**if(del==0&;&;to==0)返回(to)时出错:缺少需要TRUE/FALSE的值**,r,R,我正在运行下面的代码,它返回下面的错误-任何解决方案将不胜感激 ## 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]

我正在运行下面的代码,它返回下面的错误-任何解决方案将不胜感激

## 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 (%)")
##UQdata=仅正计数丰度数据(即仅正值)
aa2也许:

    xmin <- ceiling( min( Antho$MudUQ[Antho$Bin==1] , na.rm=TRUE) )
    xmax <- floor( max( Antho$MudUQ[Antho$Bin==1] , na.rm=TRUE) )

xmin似乎不是由
gam
生成的。在Antho$Bin或Antho$MudUQ中是否有NA
min
max
可以指定na.rm参数。是的,我在ANTO$Bin和ANTO$MudUQ中都有na参数-这些参数将如何编码?那么xmin和xmax中的[Antho$Bin==1]实际实现了什么呢?对不起,我对R代码很生疏