plot.ocm: Plot method for Continuous Ordinal Fits

In ordinalCont: Ordinal Regression Analysis for Continuous Scales

Description

Draws several summary and diagnostic plots, including the estimated g function, the estimated density function of the continuous ordinal score for the null model (no covariates), the histogram of the quantile residuals, the normal Q-Q plot and any smoother included in the model.

Usage

## S3 method for class 'ocm'
plot(
  x,
  plot.only = NULL,
  CIs = c("vcov", "no", "rnd.x.bootstrap", "fix.x.bootstrap", "param.bootstrap"),
  R = 100,
  main_gfun = "g function",
  main_density = "Density function when X=0",
  xlab = "Continuous ordinal scale [v]",
  CIcol = "lightblue",
  individual_plots = F,
  ...
)

Arguments

x	an object of class ocm
plot.only	either NULL, in which case all plots are displayed, or a value among "gfun", "quant_resid", "QQplot" or "smoother", in which case only the requested plot is displayed.
CIs	method used for confidence bands for the g function. "vcov" = Wald [default]; "no" = no CIS; "rnd.x.bootstrap" = random-x bootstrap; "fix.x.bootstrap" = bootstrap with fixed-x resampling; "param.bootstrap" = parametric bootstrap
R	the number of bootstrap replicates. Ignored if CIs="no"
main_gfun	title of the g function plot. Defauts to “g function (95% CIs)”
main_density	title of the density function plot. Defauts to “Density function when X=0”
xlab	label of the x axis for the g function and the density plots. Defaults to “Continuous ordinal scale [v]”
CIcol	color of the confidence interval bands. Defaults to “lightblue”
individual_plots	logical. If TRUE, every figure is drawn in a new window. If FALSE (default), the first four figures are drawn in a 2-by-2 array.
...	further arguments passed to or from other methods

Details

The estimated g function, quantile residual histogram and normal Q-Q plot of an ocm object are plotted. If smothers are included in the formula, the user has the option to plot them in the same graph or separately. If CIs is not "no", 95% confidence bands are also plotted.

Author(s)

Maurizio Manuguerra, Gillian Heller

See Also

ocm

Examples

fit.overall  <- ocm(overall  ~ cycleno + age + bsa + treatment, data=ANZ0001.sub, scale=c(0,100))
plot(fit.overall, CIs="vcov")
## Not run: 
plot(fit.overall, CIs="rnd.x.bootstrap", R=100)
plot(fit.overall, CIs="fix.x.bootstrap", R=100)
plot(fit.overall, CIs="param.bootstrap", R=100)

## End(Not run)

ordinalCont documentation built on Dec. 3, 2020, 1:06 a.m.