Hi Fiona,

I think you can use the function rmronm form the package mnormt to simulate from a multivariate normal distribution with a given variance-covariance matrix. In case, I wrote a function for simulating a dataset with binary response variable, you can find it here:

http://mixedpsychophysics.wordpress.com

best
Alessandro Moscatelli

PS: for the moment, the some of the functions in the MERpsychophysics code are not compatible with lme4.1 - I hope I will upload the new code version in max 20 days…
  
Il giorno 04/dic/2013, alle ore 17:48, r-sig-mixed-models-request@r-project.org ha scritto:

> simulating datasets for mixed model analysis


	[[alternative HTML version deleted]]

