R/getSmo-18-6-13.R

In gamlss: Generalized Additive Models for Location Scale and Shape

#-----------------------------------------------------------------------------------
getSmo <- function(object, what = c("mu", "sigma", "nu", "tau"),  parameter= NULL,  which=1)
{
  if (!is.gamlss(object)) stop("this is design for gamlss objects only")
  what <- if (!is.null(parameter))  {
    match.arg(parameter, choices=c("mu", "sigma", "nu", "tau"))} else  match.arg(what)               
  if (!what%in%object$par) stop(paste(what,"is not a parameter in the object","\n")) 
  #evalq(paste("object$",what,"[", which, "]",sep=""))
  AllSmo <- object[[paste(what,".coefSmo",sep="")]]
  Smo <- if (which==0) AllSmo else AllSmo[[which]]
  Smo
}