summarize {brxx}

summarize: Summarize Stan output as median, SD, and HPD quantiles

Description

This function converts raw MCMC sample data into matrix formatted summaries

Usage

summarize(Samples, nrow, ncol, CI)

Arguments

Samples

S by theta matrix of sampled parameter estimates.

nrow

Number of rows of target summary matrix

ncol

Number of columns of target summary matrix

CI

Creddible interval quantile, as a decimal (ie, for 95 percent, 0.95)

Value

Returns median, SD, and HPD CI limits

Examples

## Not run: 
your_data_s=standardize(your_data)
formatted_data=prep(your_data_s,nfactors=3)
out=sampling(model, data=formatted_data, iter=5000, seed=999)
res=as.matrix(out)
unpacked=unpack(Samples=res,Format=formatted_data)
processed=process(Loading_Matrix=unpacked$Loading_Matrix,
                 Format=formatted_data,
                 Rotate="oblimin")

summarize(processed$Loadings,
         nrow=Formatted_data$P,
         ncol=Formatted_data$Q)$Table
summarize(processed$Communality,
         nrow=Formatted_data$P,
         ncol=1)$Table
summarize(processed$Uniqueness,
         nrow=Formatted_data$P,
         ncol=1)$Table
summarize(processed$G_Factor,
         nrow=Formatted_data$P,
         ncol=1)$Table
summarize(processed$Interfactor_Correlations,
         nrow=Formatted_data$Q,
         ncol=Formatted_data$Q)$Table
summarize(processed$Omega,
         nrow=1,
         ncol=1)$Table
summarize(unpacked$Tau_Matrix,
         nrow=Formatted_data$P,
         ncol=1)$Table
## End(Not run)