## Installation

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

`install.packages("Bayesrel")`

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

`remotes::install_github("juliuspf/Bayesrel")`

## Example

### 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:

```
library(Bayesrel)
## basic example code
## load example data set from the package
## run the main reliability function
res <- strel(asrm)
## get a full result output
summary(strel)
## 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:
median(res$Bayes$samp$Bayes_alpha)
```

### 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):

```
library(Bayesrel)
## basic example code
## run the Bayesian omegas, specify 5 group factors
res <- bomegas(upps, n.factors = 5, missing = "listwise")
## get a full result output
summary(res)
## 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")`