plot.residuals {rmutil}R Documentation

Plot Residuals

Description

plot.residuals is used for plotting residuals from models obtained from dynamic models for given subsets of the data.

Usage

plot.residuals(z, x=NULL, subset=NULL, ccov=NULL, nind=NULL,
        recursive=TRUE, pch=20, ylab="Residual", xlab=NULL,
        main=NULL, ...)

Arguments

z An object of class recursive, from carma, gar, kalcount, kalseries, kalsurv, or nbkal.
x Vector of of values for the x-axis. If missing, time is used. It can also be specified by the strings "response" or "fitted".
subset A logical vector defining which observations are to be used.
ccov If the name of a time-constant covariate is supplied, separate plots are made for each distinct value of that covariate.
nind Observation number(s) of individual(s) to be plotted.
recursive If TRUE, plot recursive residuals, otherwise ordinary residuals.
others Plotting control options.

Author(s)

J.K. Lindsey

See Also

carma, gar, kalcount, kalseries, kalsurv, nbkal plot.iprofile, plot.mprofile.

Examples

library(repeated)
times <- rep(1:20,2)
dose <- c(rep(2,20),rep(5,20))
mu <- function(p) exp(p[1]-p[3])*(dose/(exp(p[1])-exp(p[2]))*
        (exp(-exp(p[2])*times)-exp(-exp(p[1])*times)))
shape <- function(p) exp(p[1]-p[2])*times*dose*exp(-exp(p[1])*times)
conc <- matrix(rgamma(40,2,scale=mu(log(c(1,0.3,0.2)))/2),ncol=20,byrow=TRUE)
conc[,2:20] <- conc[,2:20]+0.5*(conc[,1:19]-matrix(mu(log(c(1,0.3,0.2))),
        ncol=20,byrow=TRUE)[,1:19])
conc <- ifelse(conc>0,conc,0.01)
z <- gar(conc, dist="gamma", times=1:20, mu=mu, shape=shape,
        preg=log(c(1,0.4,0.1)), pdepend=0.1, pshape=log(c(1,0.2)))
plot.residuals(z, subset=1:20, main="Dose 1")
plot.residuals(z, x="fitted", subset=1:20, main="Dose 1")
plot.residuals(z, x="response", subset=1:20, main="Dose 1")

[Package rmutil version 1.0 Index]