
Categorical (visual) predictive check plots.
cat.pc(object,
dv=xvardef("dv",object),
idv=xvardef("idv",object),
level.to.plot=NULL,
subset=NULL,
histo=T,
median.line=F,
PI.lines=F,
xlb=if(histo){
paste("Proportion of ",dv)
} else {
paste(idv)
},
ylb=if(histo){
paste("Percent of Total")
} else {
paste("Proportion of Total")
},
main=xpose.create.title.text(NULL,dv,
"Predictive check of",object,subset=subset,...),
strip="Default",
...)
Xpose data object.
The dependent variable (e.g. "DV"
or "CP"
.)
The indenpent variable (e.g. "TIME"
.)
The levels to plot.
Subset of data.
If FALSE
then a VPC is created, given that idv
is defined.
Make a median line?
Make prediction interval lines?
Label for x axis.
label for y axis.
Main title.
Defining how the strips should appear in the conditioning plots.
Extra arguments passed to the function.
## read in table files
runno <- 45
xpdb <- xpose.data(runno)
## create proportion (visual) predictive check
cat.pc(xpdb,idv=NULL)
cat.pc(xpdb,idv="DOSE")
cat.pc(xpdb,idv="DOSE",histo=F)
cat.pc(xpdb,idv="TIME",histo=T,level.to.plot=1)
Run the code above in your browser using DataLab