qrLMM (version 2.3)

group.plots: Plot function for grouped data

Description

Functions for plotting a profiles plot for grouped data.

Usage

group.plot(x,y,groups,...)
group.lines(x,y,groups,...)
group.points(x,y,groups,...)

Arguments

y

the response vector of dimension \(N\) where \(N\) is the total of observations.

x

vector of longitudinal (repeated measures) covariate of dimension \(N\). For example: Time, location, etc.

groups

factor of dimension \(N\) specifying the partitions of the data over which the random effects vary.

...

additional graphical arguments passed to matplot. See par.

Author

Christian E. Galarza <> and Victor H. Lachos <>

See Also

Orthodont, Cholesterol, QRLMM

Examples

if (FALSE) {
#A full profile plot for Soybean data

data(Soybean,package = "qrNLMM")
attach(Soybean)

group.plot(x = Time,y = weight,groups = Plot,type="b",
             main="Soybean profiles",xlab="time (days)",
             ylab="mean leaf weight (gr)")

#Profile plot by genotype

group.plot(x = Time[Variety=="P"],y = weight[Variety=="P"],
             groups = Plot[Variety=="P"],type="l",col="blue",
             main="Soybean profiles by genotype",xlab="time (days)",
             ylab="mean leaf weight (gr)")

group.lines(x = Time[Variety=="F"],y = weight[Variety=="F"],
              groups = Plot[Variety=="F"],col="black")
}