2013-08-15 2 views
0

Я запускаю приведенный ниже код, и он возвращает ошибку ниже: любые решения будут очень признателен.** Ошибка в if (del == 0 && to == 0) return (to): отсутствует значение, где требуется TRUE/FALSE **

## 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 (%)") 
+1

Похоже, что это не генерируется 'gam'. У вас есть NA в Antho $ Bin или Antho $ MudUQ? 'min' и' max' могут быть заданы аргументы na.rm. –

+0

Да, у меня есть NA в Antho $ Bin и Antho $ MudUQ - как бы эти аргументы были закодированы? И что представляет собой [Antho $ Bin == 1] в xmin и xmax на самом деле? Извините, я очень зелёный, когда дело доходит до кода R – brober

ответ

1

Возможно:

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

Я думаю, что ошибка броска функции seq.

Смежные вопросы