[R-sig-ME] Fw: Fw: simulating a linear random intercept model with exogeneity assumption

[Thierry Onkelinx]

Dear Mervat,

This the r-sig-mixed-models mailing list. Not the write-code-for-me mailing
list. Please show us what you have tried. This can easily be solved using
expand.grid() and rnorm().

Best regards,

ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie & Kwaliteitszorg / team Biometrics & Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey

2017-07-11 13:28 GMT+02:00 mervat mohamed via R-sig-mixed-models <
r-sig-mixed-models op r-project.org>:

>
>
>
> ----- Forwarded Message -----From: mervat mohamed <
> mervat_moh2006 op yahoo.com>To: mailman-owner op r-project.org <
> mailman-owner op r-project.org>Sent: Sunday, July 2, 2017, 5:39:42 AM
> GMT+2Subject: Fw: [R-sig-ME] Fw: simulating a linear random intercept model
> with exogeneity assumption
>
>
>      On Wednesday, June 21, 2017 8:27 PM, mervat mohamed via
> R-sig-mixed-models <r-sig-mixed-models op r-project.org> wrote:
>
>
>
>
>   Show original message    On Monday, June 19, 2017 1:36 PM, mervat
> mohamed via R-sig-mixed-models <r-sig-mixed-models op r-project.org> wrote:
>
>
>  Hi r-sig group
> I want to know how to write the code of simulating a linear random
> intercept model with following assumptions using R programing:
> the model:
> yij=b.+b1 xij + u.j + eij,  where j refer to the group number and i refer
> to the observation number in the j group
> model assumptions:
>   1- xij ~ N (3,1.5)  2- u.j ~ N (0,1)  3 - eij ~ N (0,1)  4 - cov (xij ,
> u.j)=0  5 - cov (xij , eij)=0
>   6- cov (u.j , eij)= 0
> thnx for help
> mervat
>     [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models op r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>   _______________________________________________
> R-sig-mixed-models op r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>
>
> _______________________________________________
> R-sig-mixed-models op r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

	[[alternative HTML version deleted]]