Description Usage Arguments Methods Examples
Three plots (selectable by which
) are currently available:
a plot for the fitted model (which
=1), a plot for the individual mean profile (which
=2)
and a plot for the observed data and the corresponded mean profile
(which
=3).
1 2 3 4 |
x |
an object of class |
which |
if a subset of the plots is required, specify a subset of the numbers 1:3. Default is ( |
xlab |
label to plots. |
ylab |
label to plots. |
factor |
identify the factor for |
main |
title to plots in addition to the caption. |
subSET |
logical expression indicating elements to keep in individual mean profile plots: missing values are taken as FALSE
for |
ident |
logical expression indicating whether or not to add the number of the subject to individual mean profile plots.
The |
caption |
captions to appear above the plots. |
cex.caption |
controls the size of caption. |
signature(x="ANY", y="ANY")
:Generic function.
signature(x="cold", y="missing")
:Plot diagnostics for cold
object.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 | ##### data= datacold
### AR1R
mod1R <- cold(z ~ Time * Treatment, data = datacold, time = "Time",
id = "Subject", dependence = "AR1R")
plot(mod1R, which = 1, xlab = "Time(weeks)", ylab = "Count",
factor = Treatment, main = "Model AR1R")
plot(mod1R, which = 2, xlab = "Time(weeks)", ylab = "Count",
main = "Model AR1R")
par(mfrow = c(1, 2))
plot(mod1R, which = 2, ident = TRUE, subSET = Treatment == "0",
ylab = "Count", main = "0")
plot(mod1R, which = 2, ident = TRUE, subSET = Treatment == "1",
ylab = "Count", main = "1")
par(mfrow = c(1, 1))
par(mfrow = c(2, 2))
plot(mod1R, which = 3, subSET = (id == c(2)), xlab = "Time (weeks)",
ylab = "count", main = "0_Subject2")
plot(mod1R, which = 3, subSET = (id == c(10)), xlab = "Time (weeks)",
ylab = "Count", main = "0_Subject10")
plot(mod1R, which = 3, subSET = (id == c(26)), xlab = "Time (weeks)",
ylab = "Count", main = "1_Subject26")
plot(mod1R, which=3, subSET=(id == c(30)), xlab = "Time (weeks)",
ylab = "Count", main = "1_Subject32")
par(mfrow = c(1, 1))
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.