rbsurv方法筛选基因

 try http:// if https:// URLs are not supported
source("https://bioconductor.org/biocLite.R")
biocLite("rbsurv")
library(rbsurv)
data(gliomaSet)
gliomaSet
x <- exprs(gliomaSet)
x <- log2(x)
time <- gliomaSet$Time
status <- gliomaSet$Status

x是行代表基因，列代表样本

image.png

###################################################

code chunk number 2: rbsurv.Rnw:200-201

###################################################
fit <- rbsurv(time=time, status=status, x=x, method="efron", max.n.genes=20)

###################################################

code chunk number 3: rbsurv.Rnw:213-214

###################################################
fit$model