Nothing
rm(list=ls())
library(SEERaBomb)
load("~/data/SEER/mrgd/cancDef.RData")
load("~/data/SEER/mrgd/popsa.RData")
p=popsa%>%group_by(sex,age86)%>%summarise(py=sum(py)/1e5)
d=canc%>%select(cancer,sex,age86)%>%group_by(cancer,sex,age86)%>%summarise(cases=n())
(D=left_join(d,p))
d=D%>%mutate(incidence=cases/py)
graphics.off()
theme_set(theme_gray(base_size = 16))
# theme_update(legend.position = "top")
if(length(grep("linux",R.Version()$os))) windows <- function( ... ) X11( ... )
if(length(grep("darwin",R.Version()$os))) windows <- function( ... ) quartz( ... )
windows(width=12,height=7)
ggplot(aes(x=age86,y=incidence,col=sex),data=d)+ #log="y",
ylab("Incidence (Cases/100,000 Person-Years)")+
xlab("Age")+geom_line(size=1)+facet_wrap(~cancer) +
scale_y_log10(breaks=c(.01,1,100),labels=c(".01","1","100")) + theme( legend.position = c(.87, .0))
ggsave("~/Results/SEERaBomb/incidMat.eps")
d%>%filter(cancer%in%c("breast","lung","prostate"))%>%summarise(max(incidence))
(tb=sort(table(canc$cancer),dec=T)) #liver is top 5 world-wide, pancreas and ovary deadly, MDS since 2001
ggplot(aes(x=age86,y=incidence,col=sex),data=d%>%filter(cancer%in%names(tb)[1:9]))+
ylab("Incidence (Cases/100,000 Person-Years)")+
xlab("Age")+geom_line(size=1)+facet_wrap(~cancer) +
scale_y_log10(breaks=c(.01,1,100),labels=c(".01","1","100"))+theme(legend.position="top",strip.text=element_text(size=rel(1.4)))
ggsave("~/Results/SEERaBomb/incidMat3x3.eps")
HM=c("AML","AMLti","APL","MDS","CMML","CML","MPN","ALL","CLL","SLL","HCL","OL","NHL","MM","HL")
ggplot(aes(x=age86,y=incidence,col=sex),data=d%>%filter(cancer%in%HM))+
ylab("Incidence (Cases/100,000 Person-Years)")+
xlab("Age")+geom_line(size=1)+facet_wrap(~cancer) +
scale_y_log10(breaks=c(.01,1,100),labels=c(".01","1","100"))+theme(legend.position="top",strip.text=element_text(size=rel(1.4)))
ggsave("~/Results/SEERaBomb/incidHM.eps")
(p=popsa%>%group_by(sex,age86)%>%summarise(py=sum(py)/1e5))
(d=canc%>%select(cancer,sex,age86)%>%group_by(cancer,sex,age86)%>%summarise(cases=n()))
(D=left_join(d,p))
d=D%>%mutate(incidence=cases/py)
graphics.off()
theme_set(theme_gray(base_size = 16))
quartz(width=12,height=7)
ggplot(aes(x=age86,y=incidence,col=sex),data=d%>%filter(cancer%in%names(tb)[1:9]))+
ylab("Incidence (Cases/100,000 Person-Years)")+
xlab("Age")+geom_line(size=1)+facet_wrap(~cancer) +
scale_y_log10(breaks=c(.01,1,100),labels=c(".01","1","100"))+theme(legend.position="top",strip.text=element_text(size=rel(1.4)))
#now compare male and female incidences summed over all cancers
windows(width=7,height=7)
theme_set(theme_gray(base_size = 22))
d=canc%>%
select(sex,age86)%>%
group_by(sex,age86,add=F) %>%
summarise(cases=n())
(D=left_join(d,p))
d=D%>%mutate(incidence=cases/py)
ggplot(aes(x=age86,y=incidence,col=sex),data=d)+ #log="y",
ylab("Incidence (Cases/100,000 Person-Years)")+
xlab("Age")+geom_line(size=1)+geom_vline(xintercept = c(40,50,65),col="dark gray")+
scale_y_log10(breaks=c(.01,1,100,1000),labels=c(".01","1","100","1000")) + theme( legend.position = c(.87, .075))
ggsave("~/Results/SEERaBomb/sumMvsF.eps")
# now zoom in on screened cancers and note discontinuities at screening ages of 40, 50 and 65
windows(width=12,height=7)
theme_set(theme_gray(base_size = 22))
d=canc%>%
filter(cancer%in%c("colon","rectal","prostate","breast","CIS"))%>%
select(cancer,sex,age86)%>%
group_by(cancer,sex,age86) %>%
summarise(cases=n())
(D=left_join(d,p))
d=D%>%mutate(incidence=cases/py)
qplot(age86,incidence,col=sex,data=d,xlim=c(32,73),#log="y", #ylim=c(1,1e3),
ylab="Incidence (Cases/100,000 Person-Years)",
xlab="Age")+geom_line()+geom_vline(xintercept = c(40,50,65),col="dark gray")+
facet_wrap(~cancer) +theme(strip.text = element_text(size = rel(1.3))) +
scale_y_log10(breaks=c(10,100,1000),labels=c("10","100","1000"),limits=c(10,1000))
ggsave("~/Results/SEERaBomb/zoomScreenYrs.eps")
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.