R/summary.ame.R

Defines functions summary.ame

Documented in summary.ame

#' Summary of an AME object
#' 
#' Summary method for an AME object
#' 
#' @param object the result of fitting an AME model
#' @param ... additional parameters (not used)
#' @return a summary of parameter estimates and confidence intervals for an AME
#' fit
#' @author Peter Hoff
#' @S3method summary ame 
summary.ame <-
function(object, ...)
{ 
  fit<-object
  require(amen)
  tmp<-cbind(apply(fit$BETA,2,mean), apply(fit$BETA,2,sd) ,
       apply(fit$BETA,2,mean)/apply(fit$BETA,2,sd) , 
       2*(1-pnorm( abs(apply(fit$BETA,2,mean)/apply(fit$BETA,2,sd)))))
  colnames(tmp)<-c("pmean","psd","z-stat","p-val") 
  cat("\nRegression coefficients:\n")
  print(round(tmp,3))


  tmp<-cbind(apply(fit$VC,2,mean), apply(fit$VC,2,sd) )
  colnames(tmp)<-c("pmean","psd") 
  cat("\nVariance parameters:\n")
  print(round(tmp,3))
} 

 
 
pdhoff/amen documentation built on Dec. 21, 2020, 6:10 a.m.