[R-sig-ME] within subjects, 3 conditions, 3 lines
marKo
mtoncic at ffri.hr
Wed Jan 21 09:30:24 CET 2015
-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1
I think that you have correctly identified the problem.
Cheers,
Marko
On 01/21/2015 09:16 AM, Stanislav Aggerwal wrote:
> I will attempt an answer to my own question. Possibly the reason I
> got the same answer with the 3 methods is that my simulation was
> wrong. It simulated the responses as:
>
> y<- (cond=="a")*(a1 + a4*x) + (cond=="b")*(a1+a2 + (a4+a5)*x) +
> (cond=="c")*(a1+a3 + (a4+a6)*x) + rnorm(5*3*nsubj,mean=0,sd=.5)
>
> This contains the intercepts and slopes which are unique to each
> subject and randomly sampled from some population values. However
> the term: rnorm(5*3*nsubj,mean=0,sd=.5) is wrong because it says
> the errors are independent. Instead they should be correlated
> within each subject.
>
> I will give that a go.
>
> Stan
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1
iQIcBAEBAgAGBQJUv2OgAAoJEJcj4KySkkQsUbgP/jhOyGuwPpAyS8I6uFElxztg
sQdiGeUaVKjutcJJsHisezEcRxcf84k09hWnR1bnc/fuFQgDqS58qnoUa2JJqiZh
BrnCjIV7Cj///DLC0MLzyzvmZLq21medi6vdUwdAF5gc1GOJHoA9IWgoEACucqEc
SbpWiaqUWpNK2rCOZ3crTEBZNHEHO72qKBq670L5zDTfRfx4IXsjGMyY68Lf/Tz/
jECjDGM3byO2pw5KginLWl3hA+0ChcKNNmti3rY8p9JVZ8c4iEDrT0Hjvx/VN8zn
JPcRSJFR263haozcENATtAgUmBCT66j+IfB/7aQgAOthTdEWWjtnnnVxTbDvBkik
tVB617PXdUL5l7kcmxZYcwiigOQerJrEDcUhoL8SMtL54JluWg9MZIZxXmjGHMlm
mtrvbv6rfU3XzB80ck8zm6SSgPwP1IYoqJZ7QWBx0XTuU057tjyG8rA6F/jkYc3t
Uoe0NwY1Ox9Q62paJQ4sysyf2dgyW2RYRVVQuK1UjlH9ZfmETtq36J6yyEXjuhir
ztVsKeofNgRKo4FrmhcOLcBPofzQBnzTI1RjuqCxcbt9vsr4Wwn1uvup20OqDYb7
VKbzoKr8HR26YDS7Oc3vUiOzjszNEALnKAWFoUfW/BQRQGzIZl0lR1g52SzRphL5
u/shRpd5ZYjTJUgT9IQ+
=/B+Q
-----END PGP SIGNATURE-----
More information about the R-sig-mixed-models
mailing list