sim.step1: Simulate Data for the Predictive Distribution

In MarielleZZ/Replication: Test Replications by Means of the Prior Predictive p-Value

Description

Simulates data y.s using lavaan simulateData. The data y.s are based on samples from the posterior and represent samples from the predictive distribution.

Usage

sim.step1(n.r, nsim, post, pT, free.i, group=NULL, n.groups)

Arguments

n.r	Sample size for the new data (y.r)
nsim	Number of datasets y.s to be created.
post	Matrix with samples from the posterior distribution for the original data y.o
pT	The blavaan parameter table.
free.i	A vector with indices for the freely estimated parameters as given in pT.
group	If applicable, the name of the grouping variable.
n.groups	If applicable, the number of groups.

Value

y.s	A list containing all y.s.

Author(s)

M. A. J. Zondervan-Zwijnenburg

See Also

simulateData

Examples

## Not run: 
#the following example can be used, but takes >10 seconds

data <- data.frame(y=ChickWeight$weight,x=ChickWeight$Time)

model <- '
y ~ x     #regression
y ~1      #intercept not default in lavaan (but is in blavaan)
'

reg.post <- posterior.step1(y.o=data, model=model, nburnin=5000, nsample=5000)
reg.sim  <- sim.step1(n.r = 50, post=reg.post$post,
            pT=reg.post$pT, free.i=reg.post$free.i, nsim = 100)
head(reg.sim$y.s[[1]])

## End(Not run)

MarielleZZ/Replication documentation built on May 14, 2019, 7:55 p.m.