# NOT RUN {
# }
# NOT RUN {
library(xpose4)
## Read in the data
runno <- "57"
xpdb <- xpose.data(runno)
####################################
# here are the real data plots
####################################
kaplan.plot(x="TIME",y="DV",object=xpdb)
kaplan.plot(x="TIME",y="DV",object=xpdb,
events.to.plot=c(1,2),
by=c("DOSE==0","DOSE!=0"))
kaplan.plot(x="TIME",y="DV",object=xpdb,
events.to.plot=c(1,2),
by=c("DOSE==0","DOSE==10",
"DOSE==50","DOSE==200"))
## make a PDF of the plots
pdf(file=paste("run",runno,"_kaplan.pdf",sep=""))
kaplan.plot(x="TIME",y="DV",object=xpdb,
by=c("DOSE==0","DOSE==10",
"DOSE==50","DOSE==200"))
dev.off()
####################################
## VPC plots
####################################
kaplan.plot(x="TIME",y="DV",object=xpdb,VPC=T,events.to.plot=c(1))
kaplan.plot(x="TIME",y="DV",object=xpdb,VPC=T,
events.to.plot=c(1,2,3),
by=c("DOSE==0","DOSE!=0"))
kaplan.plot(x="TIME",y="DV",object=xpdb,VPC=T,
events.to.plot=c(1),
by=c("DOSE==0","DOSE==10","DOSE==50","DOSE==200"))
## make a PDF of all plots
pdf(file=paste("run",runno,"_kaplan.pdf",sep=""))
kaplan.plot(x="TIME",y="DV",object=xpdb,VPC=T,
by=c("DOSE==0","DOSE==10","DOSE==50","DOSE==200"))
dev.off()
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab