Re: [R] lme cant get parameter estimated correctly

2008-04-06 Thread Bill.Venables
PM To: [EMAIL PROTECTED] Subject: [R] lme cant get parameter estimated correctly I am caught in a mental trap. Why isn't the between groups variance estimated (0.0038) to be around the value with which I generated the data (0.0002)? Thanks Toby set.seed(76589437887) fph = 0.4 Sigh = sqrt

Re: [R] lme cant get parameter estimated correctly

2008-04-06 Thread toby909
Bill.Venables at csiro.au writes: Here is a demo you may like to consider. (I can see what you are trying to do with your loops, but I prefer to do it this way.) This is just for pedagogical purpose, so I like to keep the simple-minded 'for' loop. But what I really wonder is why do I not

[R] lme cant get parameter estimated correctly

2008-04-05 Thread toby909
I am caught in a mental trap. Why isn't the between groups variance estimated (0.0038) to be around the value with which I generated the data (0.0002)? Thanks Toby set.seed(76589437887) fph = 0.4 Sigh = sqrt(0.0002) Sigi = sqrt(0.04) ci = 1 fpi = matrix(,7200,3) for (i in 1:90) { fph =