nice_mcmcglmm_icc {BiostatsUHNplus} | R Documentation |
Nice table of intraclass correlation coefficients from MCMCglmm::MCMCglmm() model output
Description
Nice table of intraclass correlation coefficients from MCMCglmm::MCMCglmm() model output
Usage
nice_mcmcglmm_icc(mcmcglmm_object, prob = NULL, decimals = NULL)
Arguments
mcmcglmm_object |
returned output from MCMCglmm() |
prob |
probability for highest posterior density interval, similar to a confidence interval. Default is 0.95 (if provided) |
decimals |
number of decimal places to use in estimates |
Value
grouped_df
Examples
## Not run:
data(ae)
ae$AE_SEV_GD <- as.numeric(ae$AE_SEV_GD);
ae$Drug_1_Attribution <- 0;
ae$Drug_1_Attribution[ae$CTC_AE_ATTR_SCALE %in% c("Definite", "Probable", "Possible")] <- 1;
ae$Drug_2_Attribution <- 0;
ae$Drug_2_Attribution[ae$CTC_AE_ATTR_SCALE_1 %in% c("Definite", "Probable", "Possible")] <- 1;
prior2RE <- list(R = list(V = diag(1), fix = 1), G=list(G1=list(V=1, nu=0.02),
G2=list(V=1, nu=0.02)));
model1 <- MCMCglmm::MCMCglmm(Drug_1_Attribution ~ AE_SEV_GD + Drug_2_Attribution,
random=~ae_detail + Subject, family="categorical", data=ae, saveX=TRUE,
verbose=FALSE, burnin=2000, nitt=10000, thin=10, pr=TRUE, prior=prior2RE);
mcmcglmm_icc <- nice_mcmcglmm_icc(model1);
## End(Not run)
[Package BiostatsUHNplus version 0.0.10 Index]