祝让飞
祝让飞 - 生物信息工程师 生信一期

性别: 浙江 - 杭州 注册于 2017-03-25

偏执又爱折腾的精选伪IT男,求互粉!!! 招募医学数据分析朋友一起探讨肿瘤信息学,微信:13456826965

向TA求助
2272金币数
15620 经验值
441个粉丝
主页被访问 29217 次,10,5.10

1139 个回答

2 赞同

如何用R语言做综合2或3个基因的生存分析(TCGA数据)?

sigExpressionProfile <- profile[, sigIds]fmla <- as.formula(paste("Surv(time, status) ~", paste0(sigGenes, collapse = "+")))cox <- coxph(fmla, data = sigExpressionProfile)coef <- cox$coefficientsriskScore <- as.matrix(sigExpressionProfile) %*% coef 解释: sigExpressionProfile 表示一...

回答于 2017-09-21 14:03

0 赞同

WGCNA得到的基因模块如何做生存分析?

因为同一个模块中的基因存在共表达,可以选择该模块的特征向量作为该模块的整体表达水平,然后去做预后即可

回答于 2017-09-21 13:59

0 赞同

TCGAbiolinks安装问题求助

TCGAbiolinks貌似现在已经不能用了。。。。。

回答于 2017-09-21 13:58

0 赞同

通过简易下载器下载了TGCA的ClinicalFull_matrix.txt,怎么分析...

这个跟你的问题对胃:https://www.shengxin.ren/article/145

回答于 2017-09-19 15:50

0 赞同

批量生存分析TCGA数据不匹配的问题

第一个表达,行为基因,列为样本,你弄反了。。。。。

回答于 2017-09-19 15:48

0 赞同

GSE2034这个数据库从哪里可以看到样本的乳腺癌亚型,看别人文章...

作者没有把这部分分享出来,建议去文章的补充材料中找找,没有的话邮件作者

回答于 2017-09-18 15:23

0 赞同

TCGA下载工具检索的RNA-seq-HTseq-FPKM与RNA-seq-HTseq-FPKM-cou...

看这里:https://www.shengxin.ren/question/20

回答于 2017-09-18 15:22

2 赞同

TGCA下载的ClinicalFull_matrix.txt下来的哪一列代表生存时间?

follow_ups-follow_up 为随访记录,一般不作为最终的生存时间,你找找看有没有单独一列vital_status,一般以这列为最终的生存状态,然后其中有对应的OS,如果最终状态是Alive时,再去follow_ups里面找随访时间

回答于 2017-09-18 15:21

0 赞同

TCGA 重复样本处理

TCGA中一般重复样本还是比较少的,如果出现重复简单处理为直接选择其中一个即可

回答于 2017-09-18 15:19

0 赞同

GSE100878 这个为什么用GEO转换器导不出数据,只有样本信息的,...

你这个数据属于Illumina HiSeq 2000测序数据,并非芯片数据,解决办法是: 下载以下数据即可得到表达谱:ftp://ftp.ncbi.nlm.nih.gov/geo/series/GSE100nnn/GSE100878/suppl/GSE100878_Breast_Cancer_cell_lines_FPKM.txt.gz

回答于 2017-09-15 20:12