R/print.pencopula.R

In pencopula: Flexible Copula Density Estimation with Penalized Hierarchical B-Splines

print.pencopula <- function(x,...) {
  d <- get("d",x)
  D <- get("D",x)
  p <- get("p",x)
  ddb <- get("ddb",x)
  base <- get("base",x)
  q <- get("q",x)
  cat("\nNumber of marginal observations: ",get("n",x),"\n",sep="")
  cat("\nNumber of covariates: ",p,"\n",sep="")
  cat("\nNumber of estimated coefficients: ",get("DD",x),"\n",sep="")
  if(get("base",x)=="B-spline") cat("\nHierarchy level of marginal hierarchical B-Spline d: ",d,"\n",sep="")
  if(get("base",x)=="Bernstein") cat("\nNumber of marginal Bernstein Polynomials: ",ddb,"\n",sep="")
  if(base=="B-spline") cat("\nMaximum hierarchy level of hierarchical B-Spline basis D: ",D,"\n",sep="")
  if(base=="Bernstein") cat("\nUsed basis: ",base,"\n",sep="")
  if(base=="B-spline") cat("\nUsed basis: ",base," of degree ",q,"\n",sep="")
  names <- c()
  for(i in 1:p) names <- c(names,paste("lambda.",i,"=",get("lambda",x)[i]," ",sep=""))
  cat("\nlambda: ",names,"\n",sep="")
  cat("\nlog-likelihood: ",get("log.like",x),"\n",sep="")
  cat("\nAIC: ",get("AIC",x),"\n",sep="")
  cat("\ncAIC: ",get("cAIC",x),"\n",sep="")
  cat("\nBIC: ",get("BIC",x),"\n",sep="")
}