| Line 120... |
Line 120... |
| 120 |
## Nice for taking margins ...
|
120 |
## Nice for taking margins ...
|
| 121 |
xtabs(Freq ~ Gender + Admit, DF)
|
121 |
xtabs(Freq ~ Gender + Admit, DF)
|
| 122 |
## And for testing independence ...
|
122 |
## And for testing independence ...
|
| 123 |
summary(xtabs(Freq ~ ., DF))
|
123 |
summary(xtabs(Freq ~ ., DF))
|
| 124 |
|
124 |
|
| 125 |
## with NA's
|
125 |
## with NAs
|
| 126 |
DN <- DF; DN[cbind(6:9, c(1:2,4,1))] <- NA
|
126 |
DN <- DF; DN[cbind(6:9, c(1:2,4,1))] <- NA
|
| 127 |
DN # 'Freq' is missing only for (Rejected, Female, B)
|
127 |
DN # 'Freq' is missing only for (Rejected, Female, B)
|
| 128 |
(xtNA <- xtabs(Freq ~ Gender + Admit, DN)) # NA prints 'invisibly'
|
128 |
(xtNA <- xtabs(Freq ~ Gender + Admit, DN)) # NA prints 'invisibly'
|
| 129 |
print(xtNA, na.print = "NA") # show NA's better
|
129 |
print(xtNA, na.print = "NA") # show NAs better
|
| 130 |
xtabs(Freq ~ Gender + Admit, DN, na.rm = TRUE) # ignore missing Freq
|
130 |
xtabs(Freq ~ Gender + Admit, DN, na.rm = TRUE) # ignore missing Freq
|
| 131 |
## Use addNA = TRUE to tabulate missing factor levels:
|
131 |
## Use addNA = TRUE to tabulate missing factor levels:
|
| 132 |
xtabs(Freq ~ Gender + Admit, DN, addNA = TRUE)
|
132 |
xtabs(Freq ~ Gender + Admit, DN, addNA = TRUE)
|
| 133 |
xtabs(Freq ~ Gender + Admit, DN, addNA = TRUE, na.rm = TRUE)
|
133 |
xtabs(Freq ~ Gender + Admit, DN, addNA = TRUE, na.rm = TRUE)
|
| 134 |
## na.action = na.omit removes all rows with NAs right from the start:
|
134 |
## na.action = na.omit removes all rows with NAs right from the start:
|