假设表达矩阵data,行为基因,列为样本,则代码如下: data=data[which(apply(data,1,function(x){return(sum(x>10))})>ncol(data)*0.25),]
回答于 2019-02-18 20:21
T1,T2,T3,T4可以分别做四个0,1的变量来处理,也可以做1,2,3,4来处理
回答于 2019-02-18 20:17
每个CpG位点对应TSS 都是固定的,没必要合并,你随便找一个没有合并的甲基化谱文件里面就有这个对应的信息
回答于 2019-02-18 20:16
这个网上教程很多,看这里: http://weixin.shengxin.ren/web/wei_xin/examples/Html/10000/13151.html http://api.sangerbox.com/weixin.html?word=NCBI%2B%E6%8F%90%E4%BA%A4%E8%87%AA%E5%B7%B1%E7%9A%84%E6%B5%8B%E5%BA%8F%E6%95%B0%E6%8D%AE&page=0&ToUserName=gh_2dc29180a56b&sort=0&MAXpage=10
回答于 2019-01-23 12:03
这是因为 clinical名称只是病人的ID,RNASeq的名称后面多了三位数表示哪种组织,比如-01表示tumor,-10表示normal; 所以 随访信息的单位是 病人,RNASeq的单位是病人的某个组织,ID肯定是不同的,你自己需要做一下变换
回答于 2019-01-23 11:58