iccstan: intra-class correlations in rstan
In norimune/glmmstan: Generalized liner mixed model in Rstan

Description Usage Arguments Author(s) See Also Examples

calculating intra-class correlations(ICC) in rstan.

1	iccstan(y, group, model=NULL, iter=2000, warmup = NULL, chains= 2, thin=1, cores=NULL, parallel=TRUE)

`y`	list.
`group`	vector.
`model`	compiled stanmodel.
`iter`	rstan parameter: total number of samples, including warmup.
`warmup`	rstan parameter: number of adaptation samples.
`chains`	rstan parameter: number of chains.
`thin`	rstan parameter: A positive integer specifying the period for saving sample.
`parallel`	If TRUE, run the stanmodel in parallel using doparallel package.
`cores`	Number of cores to use when executing the chains in parallel.

Hiroshi Shimizu

glmmstan,stan

data(baseball)

#intra-class correlations
model <- iccmodel()
y <- subset(baseball,select=c(HR,HIT))
fit <- iccstan(y,baseball$team,model=model)
output_icc(fit)