R/print.qtest.mvmeta.R

In mvmeta: Multivariate and Univariate Meta-Analysis and Meta-Regression

###
### R routines for the R package mvmeta (c)
#
print.qtest.mvmeta <-
function(x, digits=3, ...) {
#
################################################################################
# 
  # FORMAT
  signif <- symnum(x$pvalue,corr=FALSE,na=FALSE,cutpoints=c(0, 
    0.001,0.01,0.05,0.1,1),symbols=c("***","**","*", ".", " "))
  Q <- formatC(x$Q,digits=digits,format="f")
  pvalue <- formatC(x$pvalue,digits=digits,format="f")
  table <- cbind(Q,df=x$df,"p-value"=pvalue,signif)
  colnames(table)[4] <- ""
#
  # PRINT
  cat(if(x$k==1L) "Uni" else "Multi","variate ","Cochran Q-test for ",
    if(x$residual) "residual ", "heterogeneity","\n",sep="")
  if(x$k>1L) cat("\nOverall test:","\n")
  cat("Q = ",Q[1]," (df = ",x$df[1],"), p-value = ",pvalue[1],"\n",sep="")
  if(x$k>1L) {
    cat("\nTests on single outcome parameters:","\n")
    print(table[-1,,drop=FALSE],quote=FALSE,right=TRUE,print.gap=2)
    cat("---\nSignif. codes: ", attr(signif, "legend"), "\n")
  }
  cat("\n")
}