Dear sir,



I am applying the package $B!H(Bgamm4$B!I(B .  I have one question about $B!F(Bgamm$B!G(B. 



I can specify smoothing parameters in $B!F(Bgam$B!G(B  such as sp=c(0.001,0.001):, 



kfit<-gam(V4~V1+s(V2,fx=F,k=10)+s(V3,fx=F,k=10),sp=c(0.001,0.001),

data=train,family=binomial(link=logit),control=gam.control(maxit=10000))



Please let me know how to specify them in $B!F(Bgamm$B!G(B or $B!F(Bgamm4$B!G(B in the below statements:



kfit<-gamm(V4~V1+s(V2,fx=F,k=10)+s(V3,fx=F,k=10),data=train,method="REML",family=binomial(link=logit),niterPQL=1000,control=lmeControl(opt="optim"),random=list(V5=~1))



kfit<-gamm4(V4~V1+s(V2,fx=F,k=10)+s(V3,fx=F,k=10),data=train,method="REML",

family=binomial(link=logit),niterPQL=2000,control=lmeControl(opt="optim",tmaxIter=500,msmaxIter=500),random=~(1|V5))



Best regards,



Prof. M.Tsujitani

Osaka electro-Communication University.



	[[alternative HTML version deleted]]

