##----eval=FALSE,echo=TRUE----------------------------------------------------#如果(!require('PDATK'))BiocManager::install('PDATAK')##----加载包-------------------------------------------------------------库(PDATK)##----生存实验_结构-------------------------------------------#--创建一些伪数据#化验分析1<-基质(rnorm(100),nrow=10,ncol=10,dimnames=列表(粘贴0('gene_',seq_len(10)),粘贴0('示例_',seque_len#列和行元数据rowMData<-DataFrame(gene_name=行名称(分析1),id=seq_len(10),row.names=行名称(分析1))colMData<-DataFrame(sample_name=列名(分析1),overall_survival=样本.int(1000,10),os_status=样本(c(0L,1L),10,替换=真),row.names=列名(分析1))#——用它来构建生存实验surv实验<-生存实验(分析=简单列表(rna=分析1),rowData=rowMData,colData=colMData、metadata=list(a='Some metadata'),survival_time='overall_survival',event_occurred='os_status')## -----------------------------------------------------------------------------#——建立总结性实验sumExperience<-总结实验(分析=简单列表(rna=分析1),rowData=rowMData,colData=colMData、metadata=list(a='Some meta data'))#--将其转换为生存实验#使用sumExp参数,该参数必须命名survExperiment<-生存实验(sumExp=sumExperiment,survival_time='overall_survival',event_occurred='os_status')##----队列列表_构造函数---------------------------------------------------cohortList<-Cohort列表(列表(cohort1=survExperiment,cohort2=survExperiment),mDataTypes=c('rna_seq','rna_micro')##----生存模型构造器------------------------------------------------set.seed(1987)surveModel<-生存模型(surveExperiment,randomSeed=1987)