load('GSE38942_eSet.Rdata')
class(gset)
gset[[1]] ##降级提取gset
class(gset[[1]])
dat1<-exprs(gset[[1]])##使用函数exprs获取样本表达矩阵
pd1<-pData(gset[[1]]) ##使用函数pData获取样本临床信息（如性别、年龄、肿瘤分期等等）
group_list1<-ifelse(pd1$source_name_ch1=="colon_NH_placebo"|pd1$source_name_ch1=="colon_NMS_placebo",
'CTR','NMS')
group <- names(table(pd1$source_name_ch1) )
group
mygroup <- c("CTRL","NMS","NMS+V")
group_list1 <- 1:length(pd1$source_name_ch1)
for (j in 1:length(group)){ group_list1[which(pd1$source_name_ch1==group[j] )]<- mygroup[j] }