setwd(“G:\Bcell_bladder\GEO��Ϊѵ����\GSE31684”)
lnc <- read.csv(“Bcell����������ϸ��ϵ���������.csv”,stringsAsFactors = F)
data <- read.table(“GSE31684_lnc.txt”,stringsAsFactors = F)
linchuang <- read.csv(“GSE31684_clinical.csv”,stringsAsFactors = F,row.names = 1)
linshi <- apply(data,1,scale)
linshi <- t(linshi)
colnames(linshi) <- colnames(data)
data <- linshi
linchuang <- linchuang[colnames(data),]
name <- intersect(rownames(data),lnc[,1])
data <- data[name,]
cox <- matrix(,,4)
colnames(cox) <- c(‘exp(coef)’, ‘exp(-coef)’, ‘lower .95’, ‘upper .95’)
name_cox <- c()
for (i in 1:dim(data)[1]){
a <- summary(coxph(Surv(linchuangevent) ~ data[i,]))
if (a%7B%0A%20%20%20%20name_cox%20%3C-%20c(name_cox%2Crownames(data)%5Bi%5D)%0A%20%20%20%20cox%20%3C-%20rbind(cox%2Ca#card=math&code=logtest%5B3%5D%20%3C%200.05%29%7B%0A%20%20%20%20name_cox%20%3C-%20c%28name_cox%2Crownames%28data%29%5Bi%5D%29%0A%20%20%20%20cox%20%3C-%20rbind%28cox%2Ca)conf.int)
}
else {
name_cox <- name_cox
cox <- cox
}
}