# NOT RUN {
# load whitening library
library("whitening")
# load TGCA LUSC data set
data(lusc)
names(lusc)
#"rnaseq2" "methyl" "sex" "packs"
#"survivalTime" "censoringStatus"
dim(lusc$rnaseq2) # 130 206 gene expression
dim(lusc$methyl) # 130 234 methylation level
# }
# NOT RUN {
library("survival")
s = Surv(lusc$survivalTime, lusc$censoringStatus)
plot(survfit(s ~ lusc$sex), xlab = "Years", ylab = "Probability of survival", lty=c(2,1), lwd=2)
legend("topright", legend = c("male", "female"), lty =c(1,2), lwd=2)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab