uvsdSim {hbmem}R Documentation

Function uvsdSim

Description

Simulates data from a hierarchical UVSD model.

Usage

uvsdSim(I = 30, J = 200, K = 6, muN = -0.5, s2aN = 0.2, s2bN = 0.2,
    muS = 0.5, s2aS = 0.2, s2bS = 0.2, muS2 = 0, s2aS2 = 0.2,
    s2bS2 = 0.2, crit = matrix(rep(c(-1.5, -0.5, 0, 0.5, 1.5),
        each = I), ncol = (K - 1)))

Arguments

I Number of participants.
J Number of items.
K Number of response options.
muN Mean of new-item distribution
s2aN Variance of participant effects on mean of new-item distribution.
s2bN Variance of item effects on mean of new-item distribution.
muS Mean of studied-item distribution
s2aS Variance of participant effects on mean of studied-item distribution.
s2bS Variance of item effects on mean of studied-item distribution.
muS2 Mean variance of studied-item distribution, sigma2
s2aS2 Variance of participant effects sigma2.
s2bS2 Variance of item effects on sigma2.
crit Matrix of criteria (not including -Inf or Inf). Columns correspond to criteria, rows correspond to participants.

Value

The function returns an internally defined "uvsdSim" structure.

Author(s)

Michael S. Pratte

References

See Pratte, Rouder, & Morey (2009)

See Also

hbmem

Examples

#Data from hiererchial model
sim=uvsdSim(s2aS2=0,s2bS2=0) #one Sigma2 
slotNames(sim)
table(sim@resp,sim@cond)

#Usefull to make data.frame for passing to functions
dat=as.data.frame(cbind(sim@subj,sim@item,sim@cond,sim@lag,sim@resp))
colnames(dat)=c("sub","item","cond","lag","resp")

table(dat$resp,dat$cond)

[Package hbmem version 0.1 Index]