Good afternoon, Given lme4 package in R, I often show a random-intercept model (i.e., y~1 + (1|grouping)) using a picture like: https://github.com/hkil/m/blob/master/mlm.PNG. I wonder how to extend this picture to show a fully crossed random-intercepts model of the form: y~1 + (1|grouping_1) + (1|grouping_2) ? Specifically, do I need to have 2 red, level-2 distributions (side-by-side) each for showing the random intercepts of each grouping_j and if so, is there any way to show the "crossing*"* of the two random-effects in the picture? Many thanks, Simon
graphics to show crossed random-effects
3 messages · Simon Harmel, Thierry Onkelinx
Dear Simon, I think you need to plot the crossed random effects on different axes to illustrate their independence. E.g. the first on the x-axis, the second on the y-axis and the density on the z-axis. Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance thierry.onkelinx at inbo.be Havenlaan 88 bus 73, 1000 Brussel www.inbo.be /////////////////////////////////////////////////////////////////////////////////////////// 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 /////////////////////////////////////////////////////////////////////////////////////////// <https://www.inbo.be> Op do 16 jul. 2020 om 20:34 schreef Simon Harmel <sim.harmel at gmail.com>:
Good afternoon, Given lme4 package in R, I often show a random-intercept model (i.e., y~1 + (1|grouping)) using a picture like: https://github.com/hkil/m/blob/master/mlm.PNG. I wonder how to extend this picture to show a fully crossed random-intercepts model of the form: y~1 + (1|grouping_1) + (1|grouping_2) ? Specifically, do I need to have 2 red, level-2 distributions (side-by-side) each for showing the random intercepts of each grouping_j and if so, is there any way to show the "crossing*"* of the two random-effects in the picture? Many thanks, Simon [[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
Dear Thierry, Thanks, what do you think about this one: https://github.com/hkil/m/blob/master/Cros.PNG Simon On Fri, Jul 17, 2020 at 8:32 AM Thierry Onkelinx <thierry.onkelinx at inbo.be> wrote:
Dear Simon, I think you need to plot the crossed random effects on different axes to illustrate their independence. E.g. the first on the x-axis, the second on the y-axis and the density on the z-axis. Best regards, ir. Thierry Onkelinx Statisticus / Statistician Vlaamse Overheid / Government of Flanders INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND FOREST Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance thierry.onkelinx at inbo.be Havenlaan 88 bus 73, 1000 Brussel www.inbo.be /////////////////////////////////////////////////////////////////////////////////////////// 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 /////////////////////////////////////////////////////////////////////////////////////////// <https://www.inbo.be> Op do 16 jul. 2020 om 20:34 schreef Simon Harmel <sim.harmel at gmail.com>:
Good afternoon, Given lme4 package in R, I often show a random-intercept model (i.e., y~1 + (1|grouping)) using a picture like: https://github.com/hkil/m/blob/master/mlm.PNG. I wonder how to extend this picture to show a fully crossed random-intercepts model of the form: y~1 + (1|grouping_1) + (1|grouping_2) ? Specifically, do I need to have 2 red, level-2 distributions (side-by-side) each for showing the random intercepts of each grouping_j and if so, is there any way to show the "crossing*"* of the two random-effects in the picture? Many thanks, Simon [[alternative HTML version deleted]]
_______________________________________________ R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models