Nothing
library(SEERaBomb)
load("~/data/SEER/mrgd/cancDef.RData") #loads in canc
d=canc%>%filter(cancer=="CML")
d=d%>%mutate(status=as.numeric(COD>0),surv=(surv+0.5)/12)
d=d%>%select(yrdx,agedx,sex,surv,status) #SEER CML survival times
# mkMrt("username", "password") #use your Human Mortality Database account info
load("~/data/usMort/mrt.RData")
brkst=c(0,0.5,1,2,3,4,5,6,8)
brksy=c(1972,1990,2002,2015)
D=msd(d,mrt,brkst,brksy) #msd() = Mortality Since Diagnosis
head(D,2)
quartz(width=6,height=5)
# quartz(width=6,height=5)
theme_update(legend.position = c(.85, .37),axis.text=element_text(size=rel(1.4)),
axis.title=element_text(size=rel(1.4)),strip.text = element_text(size = rel(1.5)))
# theme_update(legend.position=c(.48,.43),axis.text=element_text(size=rel(1.4)),
# axis.title=element_text(size=rel(1.4)),strip.text=element_text(size=rel(1.5)))
g=qplot(x=t,y=RR,data=D,col=Years,geom=c("line","point"), ylim=c(0,18),
xlab="Years Since CML Diagnosis",ylab="Relative Risk of Mortality")
g=g+geom_line(size=1)+facet_grid(sex~.)+geom_abline(intercept=1, slope=0)
g+geom_errorbar(aes(ymin=rrL,ymax=rrU,width=.15))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.