## File Name: parmsummary_extend.R
## File Version: 0.12
####################################################
# extend parameter summary
parmsummary_extend <- function( dfr, level=.95,
est_label="est", se_label="se",df_label="df" )
{
dfr <- as.data.frame(dfr)
#-- compute t values and p values
dfr$t <- dfr[,est_label] / dfr[,se_label]
dfr$p <- 2 * stats::pnorm( - abs(dfr$t) )
#-- compute confidence intervals
if ( ! is.null(level) ){
if ( "df" %in% colnames(dfr) ){
quant <- - stats::qt( (1-level)/2, df=dfr$df)
} else {
quant <- - stats::qnorm( (1-level)/2 )
}
dfr[, paste0("lower", 100*level ) ] <-
dfr[,est_label] - quant * dfr[, se_label ]
dfr[, paste0("upper", 100*level ) ] <-
dfr[,est_label] + quant * dfr[, se_label ]
}
#-- output
return(dfr)
}
#####################################################
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.