R build status


You can install the released version of Bayesrel from CRAN with:


or install the latest version of Bayesrel from [github] (https://github.com) with the help of the remotes-package:



Unidimensional data

This is a basic example which shows you how to compute alpha, lambda2, the glb, and omega for an example real data set:

## basic example code
## load example data set from the package
## run the main reliability function
res <- strel(asrm)
## get a full result output
## return the probability that coefficient alpha is larger than .70
pStrel(res, estimate = "alpha", low.bound = .70)

## get the posterior median of, e.g., alpha instead of the mean:

Multidimensional data

This is a basic example which shows you how to compute omega_t and omega_h for an example real data set. The data follow a second-order factor model with no crossloadings (required):

## basic example code
## run the Bayesian omegas, specify 5 group factors
res <- bomegas(upps, n.factors = 5, missing = "listwise")
## get a full result output
## return the probability that coefficient omega_t is larger than .70
pOmegas(res, cutoff.t = .70)
## plot posterior predictive check for the higher-order (second-order) factor model
secoFit(res, upps)

In the example above we implicitly assumed that the items of the data set were ordered so that, with 5 group factors, the first four items load on the first factor, items 5-8 load on the second factor and so on. When the data is not organized this way and/or the items cannot be distributed among the factors evenly, one can specify a model syntax relating the items to the group factors in lavaan style. The item names need to equal the variable names in the data:

model <- "
  f1 =~ U17_r + U22_r + U29_r + U34_r
  f2 =~ U4 + U14 + U19 + U27
  f3 =~ U6 + U16 + U28 + U48
  f4 =~ U23_r + U31_r + U36_r + U46_r
  f5 =~ U10_r + U20_r + U35_r + U52_r

The reliability is then estimated as follows:

res <- bomegas(upps, n.factors = 5, model = model, missing = "listwise")