r - GAM() error: unused argument(s) (bs = "ps", k = 8)? -

after loading mgcv , running below model - returns error below.. code has worked in previous times (i.e. yesterday). appreciated.

> aa1<-gam(bin~s(mud,bs="ps",k=8),family=binomial, gamma=1,data=abaren)  error in s(mud, bs = "ps", k = 8) : unused argument(s) (bs = "ps", k = 8) 

--

i'm getting error below.. seems rather strange given fact of code running 2 days ago..

## uq abundance 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 && == 0) return(to) :                  missing value 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,       lab="density per 0.0132 m2", xlab="mud content (%)") 

are using wrong package? both mgcv , gam packages have function s(), former's takes bs , k argument influence basis functions.

for example, can reproduce error using gam::gam():

> require("gam") > data(kyphosis) > gam(kyphosis ~ s(age, k = 4, bs = "ps") + number, family = binomial,  +     data=kyphosis, trace = true) error in s(age, k = 4, bs = "ps") : unused arguments (k = 4, bs = "ps") > gam(kyphosis ~ s(age) + number, family = binomial, data=kyphosis, trace = true) gam s.wam loop 1: deviance = 66.42095  gam s.wam loop 2: deviance = 63.77252  gam s.wam loop 3: deviance = 63.25199  gam s.wam loop 4: deviance = 63.13399  gam s.wam loop 5: deviance = 63.11016  gam s.wam loop 6: deviance = 63.10748  gam s.wam loop 7: deviance = 63.10727  gam s.wam loop 8: deviance = 63.10725  gam s.wam loop 9: deviance = 63.10725  call: gam(formula = kyphosis ~ s(age) + number, family = binomial,      data = kyphosis, trace = true)  degrees of freedom: 80 total; 75.00002 residual residual deviance: 63.10725 

but not mgcv package , mgcv::gam():

> require(mgcv) > require("mgcv") > set.seed(2) ## simulate data...  > dat <- gamsim(1,n=400,dist="normal",scale=2) gu & wahba 4 term additive model > b <- gam(y ~ s(x0, k = 5, bs = "ps") + s(x1) + s(x2) + s(x3), data=dat) loading required package: splines