IP屬地:山東
用這個(gè)代碼做出來(lái)文件有點(diǎn)小,而生信自學(xué)網(wǎng)中用的perl腳本大概175m左右,這個(gè)代碼跑出來(lái)90m左右。這是什么問(wèn)題呢,????
TCGA基因表達(dá)數(shù)據(jù)和生存數(shù)據(jù)合并打開臨床數(shù)據(jù)excel,分別自定義排序futime和fustat,futime中去掉生存天數(shù)為0和沒(méi)有的數(shù)據(jù),fustat中去掉無(wú)生存狀態(tài)的數(shù)據(jù)(1為死亡,0為存活) 將上圖...
為什么輸出的文件不全呀,還是有哪里不對(duì)
基因表達(dá)矩陣和生存時(shí)間的合并通常來(lái)說(shuō),GEO和TCGA的基因表達(dá)數(shù)據(jù)和臨床數(shù)據(jù)是分開的,為了合并這兩個(gè)數(shù)據(jù),這里使用R進(jìn)行操作。 基因表達(dá)數(shù)據(jù)如下: 生存數(shù)據(jù)如下: 在R中讀入上述兩個(gè)數(shù)據(jù),然后使用代碼...