[R] nested hierarchical design

From: Jelmer Elzinga <JelmerAnne.Elzinga_at_unil.ch>
Date: Tue 16 Jan 2007 - 17:00:41 GMT


Dear R-Helpers,

I would like to know what syntax I need to use to do a nested anova for 1. a continuous variable and 2. count data (x out of y) 1. The first I used to do in SPSS and I would like to be able to do it in R as well.
This is the hierarchical model I would like to use: a continuous variable explained by factor A(fixed) + factor B(random) nested in A + factor C (random) nested in factor B (which is nested in A).
The data is unbalanced at all levels, meaning different nrs of B per A, different nrs of C for B, different nrs of values per C In SPSS this seems to be quite simple with a small writing in the syntax (see below). But how do I get the same results in R with lme?

SPSS syntax:
UNIANOVA
  var BY A B C
/RANDOM = B C
/METHOD = SSTYPE(3)
/INTERCEPT = INCLUDE
/CRITERIA = ALPHA(.05)
/DESIGN = A B(A) C(B(A))).
2. the same model but than for count data (like 15 out of 30, 23 out of 60) instead of the continous variable(I know the basics of glm in R)

Thanks a lot for your help
Jelmer



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Wed Jan 17 04:29:10 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Wed 17 Jan 2007 - 10:30:24 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.