The R Project SVN R

Rev

Rev 74460 | Rev 86402 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 74460 Rev 74654
Line 134... Line 134...
134 
 
 134 
 
135 
\donttest{if(require("Matrix")) withAutoprint({
 135 
\donttest{if(require("Matrix")) withAutoprint({
136 
 ## similar to "nlme"s  'ergoStool' :
 136 
 ## similar to "nlme"s  'ergoStool' :
137 
 d.ergo <- data.frame(Type = paste0("T", rep(1:4, 9*4)),
 137 
 d.ergo <- data.frame(Type = paste0("T", rep(1:4, 9*4)),
138 
                      Subj = gl(9, 4, 36*4))
 138 
                      Subj = gl(9, 4, 36*4))
139 
 print(xtabs(~ Type + Subj, data = d.ergo)) # 4 replicates each
 139 
 xtabs(~ Type + Subj, data = d.ergo) # 4 replicates each
140 
 set.seed(15) # a subset of cases:
 140 
 set.seed(15) # a subset of cases:
141 
 print(xtabs(~ Type + Subj, data = d.ergo[sample(36, 10), ], sparse = TRUE))
 141 
 xtabs(~ Type + Subj, data = d.ergo[sample(36, 10), ], sparse = TRUE)
142 
 
 142 
 
143 
 ## Hypothetical two-level setup:
 143 
 ## Hypothetical two-level setup:
144 
 inner <- factor(sample(letters[1:25], 100, replace = TRUE))
 144 
 inner <- factor(sample(letters[1:25], 100, replace = TRUE))
145 
 inout <- factor(sample(LETTERS[1:5], 25, replace = TRUE))
 145 
 inout <- factor(sample(LETTERS[1:5], 25, replace = TRUE))
146 
 fr <- data.frame(inner = inner, outer = inout[as.integer(inner)])
 146 
 fr <- data.frame(inner = inner, outer = inout[as.integer(inner)])
147 
 print(xtabs(~ inner + outer, fr, sparse = TRUE))
 147 
 xtabs(~ inner + outer, fr, sparse = TRUE)
148 
})}% only if Matrix is available
 148 
})}% only if Matrix is available
149 
}
 149 
}
150 
\keyword{category}
 150 
\keyword{category}