View source: R/plot.oceNPMLE.R
plot.oceNPMLE | R Documentation |
oceNPMLE
object.Plot oceNPMLE
object.
## S3 method for class 'oceNPMLE'
plot(
x,
xlab = "Ordering Score",
ylab = "Proportion with a larger ordering score",
ylim = c(0, 1),
col = c("red", "blue"),
mark.time = TRUE,
...
)
x |
oceNPMLE object (see |
xlab |
x label |
ylab |
y label |
ylim |
limits for the y axis, defaults to c(0,1) |
col |
color vector, col[1] for group=0 and col[2] for group=1. |
mark.time |
logical, should censored values be plotted? |
... |
Extra arguments (e.g., lwd=2) added to |
No return value, called for side effects.
data(simScenario5)
dataFormt<-oceFormat(data=simScenario5, oceTime=c("T1","T2","T3"),
oceStatus=c("I1","I2","I3"), group=c("Z"),
oceNames = c("Death","Stroke/MI","Bleed"))
npmleOutput<- oceNPMLE(dataFormt)
plot(npmleOutput, xlab="Custom x label", mark.time=FALSE, lwd=2)
# can add lines from coxph output
coxOutput<- oceCoxph(dataFormt)
plot(coxOutput,linesonly=TRUE, col=c("orange","purple"),lwd=2)
legend("bottomleft",
legend=c("grp=0, NPMLE","grp=1, NPMLE","grp=0, coxph","grp=1, coxph"),
col=c("red","blue","orange","purple"),lty=c(1,1,1,1),lwd=2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.