R/bigssa.R

bigssa <-
  function(formula,data=NULL,type=NULL,nknots=NULL,rparm=NA,
           lambdas=NULL,skip.iter=TRUE,se.fit=FALSE,rseed=1234,
           gcvopts=NULL,knotcheck=TRUE,gammas=NULL,weights=NULL,
           random=NULL,remlalg=c("FS","NR","EM","none"),
           remliter=500,remltol=10^-4,remltau=NULL) {
    ###### Fits Smoothing Spline ANOVA models
    ###### Nathaniel E. Helwig (helwig@umn.edu)
    ###### Last modified: August 10, 2015
    
    if(class(data)[1]=="makessa"){
      ssafit <- ssawork(formula,data)
    } else{
      ssamk <- makessa(formula,data,type,nknots,rparm,
                       lambdas,skip.iter,se.fit,rseed,
                       gcvopts,knotcheck,gammas,weights,
                       random,remlalg,remliter,remltol,remltau)
      ssafit <- ssawork(formula,ssamk)
    }
    ssafit <- c(ssafit,list(call=formula))
    class(ssafit) <- "bigssa"
    return(ssafit)
    
  }
taylerablake/thin-plate-splines documentation built on May 8, 2019, 11:16 p.m.