SimComp (version 2.2)

plot.SimCi: Plot function for SimCi-objects

Description

A plot of the results of SimCiDiff and SimCiRat, respectively.

Usage

# S3 method for SimCi
plot(x, xlim, xlab, ylim, ...)

Arguments

x

an object of class "SimCi" as obtained by calling SimCiDiff or SimCiRat

xlim

a numeric vector of length 2, giving the x coordinate range

xlab

a title for the x axis

ylim

a numeric vector of length 2, giving the y coordinate range

…

arguments to be passed to plot

Value

A plot of the confidence intervals of a "SimCi" object.

See Also

SimCiDiff, SimCiRat

Examples

# NOT RUN {
# Example 1:
# Simultaneous confidence intervals related to a comparison of the groups
# B and H against the standard S, on endpoint Thromb.count, assuming unequal
# variances for the groups. This is an extension of the well-known Dunnett-
# intervals to the case of heteroscedasticity.

data(coagulation)

interv1 <- SimCiDiff(data=coagulation, grp="Group", resp="Thromb.count",
  type="Dunnett", base=3, alternative="greater", covar.equal=FALSE)
interv1
plot(interv1)

# Example 2:
# Simultaneous confidence intervals related to a comparisons of the groups
# B and H against the standard S, simultaneously on all endpoints, assuming
# unequal covariance matrices for the groups. This is an extension of the well-
# known Dunnett-intervals to the case of heteroscedasticity and multiple
# endpoints.

data(coagulation)

interv2 <- SimCiDiff(data=coagulation, grp="Group", resp=c("Thromb.count","ADP","TRAP"),
  type="Dunnett", base=3, alternative="greater", covar.equal=FALSE)
summary(interv2)
par(mfrow=c(1,3)); plot(interv2)

# Example 3:
# Simultaneous confidence intervals for ratios of means, related to an all-pair
# comparison of the groups B, H and S, simultaneously on all endpoints, assuming unequal
# covariance matrices for the groups. This is an extension of the well-known Tukey-
# intervals to the case of heteroscedasticity and multiple endpoints, and in terms of
# ratios of means instead of differences.

data(coagulation)

interv3 <- SimCiRat(data=coagulation, grp="Group", resp=c("Thromb.count","ADP","TRAP"),
  type="Tukey", alternative="two.sided", covar.equal=FALSE)
summary(interv3)
par(mfrow=c(3,1)); plot(interv3)
# }