random.effects: Extract the random effects of the model

In Bayesthresh: Bayesian thresholds mixed-effects models for categorical data

Description

Extract the random effects of the model of the class Bayesthres

Usage

random.effects(object, HPDinterval=FALSE, prob=0.95)

Arguments

object	is an object of class "Bayesthres"
HPDinterval	If HPDinterval=TRUE, return the credibility interval for fixed effects model. By default is FALSE. See package coda
prob	A numeric scalar in the interval (0,1) giving the target probability content of the intervals. The nominal probability content of the intervals is the multiple of 1/nrow(obj) nearest to prob

Examples

# Not run
data(sensory)

Consumer <- factor(sensory$consumer)
Sacarose <- factor(sensory$sacarose)

#### Model 
# Not run
dex1 <- Bayesthresh(cor ~ (1|Consumer) + Sacarose, burn = 0, Write=TRUE,
																			jump = 1, ef.iter = 10, data=sensory) 
random.effects(dex1)
random.effects(dex1, HPDinterval=TRUE)

Bayesthresh documentation built on May 30, 2017, 6:24 a.m.