小果的单日小技巧DAY4:预测某样本亚型对免疫治疗的响应step2
尔云间 一个专门做科研的团队
欢迎点赞+收藏+关注
生信人R语言学习必备
立刻拥有一个Rstudio账号
开启升级模式吧
(56线程,256G内存,个人存储1T)

上一期小果分享了利用生信云平台进行某样本亚型对免疫治疗的响应的预测,废话不多说,我们继续吧(仍旧是在服务器进行操作哦)。
食用步骤
Step1:Rscript preTIDE.R 亚型转录组数据进行标准化。
Step2:Rscript readTIDE.R 表型文件 TIDE输出 样本转录组 样本信息 TCGA中该类癌症全部样本转录组数据(剔除低表达)
Step3: Rscript readsubmap.R submap网站得到的矩阵中的16个pvalu
加餐~
1. 点击上传数据 输入***文件,输入的文件名必须和示例文件名一致,且都是英文无特殊字符。如果使用自有数据时,需注意文件名和对应的内容要匹配,以及输入数据里的行名列名等等。
2. 如有非文件类参数,点击 输入参数对应的值,默认比如 pvalue值是0.05 gene值是 TP53
【输入文件名】

所需文件展示
easy_input_anno.txt这是记录亚型样本的分型信息的文本文件↓

TIDE网站生成的输出文件↓

这是自有样本的转录组文件↓

记录自有转录组信息的文件,不改变列名,程序需要利用列名sample与label进行检索↓

这是TCGA中该类癌症的全部样本数据,应为fpkm格式,并且剔除了均值表达小于1的低表达基因↓

【生成文件展示】
skcm.immunotherapy.for.SubMap.cls

skcm.immunotherapy.for.SubMap.gct

Immune2.for.SubMap.cls

mmune2.for.SubMap.gct

以上这四个文件用于后期submap
嗝~吃得差不多了,我们下期见~

“生信果”,生信入门、R语言、生信图解读与绘制、软件操作、代码复现、生信硬核知识技能、服务器、生物信息学的教程,以及基于R的分析和可视化等原创内容,一起见证小白和大佬的成长。