| Line 105... |
Line 105... |
| 105 |
%% esophNA <- esoph; esophNA$ncases[66] <- NA
|
105 |
%% esophNA <- esoph; esophNA$ncases[66] <- NA
|
| 106 |
%% with(esophNA, {
|
106 |
%% with(esophNA, {
|
| 107 |
%% interaction.plot(agegp, alcgp, ncases/ncontrols, col = 2:5)
|
107 |
%% interaction.plot(agegp, alcgp, ncases/ncontrols, col = 2:5)
|
| 108 |
%% # doesn't show *last* group either
|
108 |
%% # doesn't show *last* group either
|
| 109 |
%% interaction.plot(agegp, alcgp, ncases/ncontrols, col = 2:5, type = "b")
|
109 |
%% interaction.plot(agegp, alcgp, ncases/ncontrols, col = 2:5, type = "b")
|
| 110 |
%% ## alternative take non-NA's {"cheating"}
|
110 |
%% ## alternative take non-NAs {"cheating"}
|
| 111 |
%% interaction.plot(agegp, alcgp, ncases/ncontrols, col = 2:5,
|
111 |
%% interaction.plot(agegp, alcgp, ncases/ncontrols, col = 2:5,
|
| 112 |
%% fun = function(x) mean(x, na.rm = TRUE),
|
112 |
%% fun = function(x) mean(x, na.rm = TRUE),
|
| 113 |
%% sub = "function(x) mean(x, na.rm=TRUE)")
|
113 |
%% sub = "function(x) mean(x, na.rm=TRUE)")
|
| 114 |
%% })
|
114 |
%% })
|
| 115 |
%% rm(esophNA) # to clear up
|
115 |
%% rm(esophNA) # to clear up
|