gamfit {spgam}R Documentation

~~function to do ... ~~

Description

~~ A concise (1-5 lines) description of what the function does. ~~

Usage

gamfit(form, gam.data, pts, region, h = 0, hvals = NA, ngrid = 4000)

Arguments

form ~~Describe form here~~
gam.data ~~Describe gam.data here~~
pts ~~Describe pts here~~
region ~~Describe region here~~
h ~~Describe h here~~
hvals ~~Describe hvals here~~
ngrid ~~Describe ngrid here~~

Details

~~ If necessary, more details than the description above ~~

Value

~Describe the value returned If it is a LIST, use

comp1 Description of 'comp1'
comp2 Description of 'comp2'

...

Warning

....

Note

~~further notes~~

~Make other sections like Warning with section{Warning }{....} ~

Author(s)

~~who you are~~

References

~put references to the literature/web site here ~

See Also

~~objects to See Also as help, ~~~

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--    or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function(form,gam.data,pts,region,h=0,hvals=NA,ngrid=4000){
 formul=formula(form)
 termos<-terms(formul)
 at<-as.character(attr(termos,"variables"))[-1]
 data<-as.data.frame(gam.data[,at])
 why<-data[,at[1]]
 ny<-nrow(data)
 g2=rep(0,ny)
 fit=glm(formul,family=binomial,data=data)
 g1=predict(fit)
 fhat <- g1 + g2
 phat <- exp(fhat)/(1+exp(fhat))
 w <- phat*(1-phat)
 z <- fhat + (why-phat)/w
  
if(h!=0){
new=waddmod(formul,data,pts,g2=g2,z=z,w=w,hopt=h,hvals=hvals,region=region,sameh=T);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=h,hvals=hvals,region=region,sameh=T);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=h,hvals=hvals,region=region,sameh=T);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=h,hvals=hvals,region=region,sameh=T);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=h,hvals=hvals,region=region,sameh=T);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=h,hvals=hvals,region=region,ngrid=ngrid,sameh=T)
  }
else{
new=waddmod(formul,data,pts,g2=g2,z=z,w=w,hvals=hvals,region=region);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=old$h,hvals=hvals,region=region);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=old$h,hvals=hvals,region=region);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=old$h,hvals=hvals,region=region);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=old$h,hvals=hvals,region=region,sameh=T);old=new
new=waddmod(formul,data,pts,g2=old$g2,z=old$z,w=old$w,hopt=old$h,hvals=hvals,region=region,ngrid=ngrid,sameh=T)
  }
new
  }

[Package spgam version 1.0 Index]