pyscenic的結(jié)果不可重復(fù)

一、pyscenic的結(jié)果不可重復(fù)

最近在運(yùn)行pyscenic的結(jié)果的時(shí)候發(fā)現(xiàn)了兩次相同參數(shù)相同數(shù)據(jù)的結(jié)果存在不一致霜运,回溯發(fā)現(xiàn)scenic的第一步GRN的結(jié)果已經(jīng)出現(xiàn)了差異挑童,轉(zhuǎn)錄因子和目標(biāo)基因之間的importance發(fā)生了變化序六,同時(shí)整體的importance排序也已經(jīng)不一樣了,如下:


GRN結(jié)果實(shí)例

同樣的環(huán)境悠就、數(shù)據(jù)和腳本千绪,結(jié)果不一樣,就會(huì)懷疑是隨機(jī)因素導(dǎo)致的梗脾,雖然中文檢索沒有人解答荸型,但是在github,pyscenic的團(tuán)隊(duì)還是給出了解釋:


不可重復(fù)性的解釋

大概意思是運(yùn)行scenic第一個(gè)部分的GENIE3/GRNBoost本身就會(huì)存在隨機(jī)性炸茧,作者也是推薦大家為了保證的可靠性瑞妇,可以通過多次運(yùn)行取均值的方法來解決。

好梭冠,既然不是自己的問題辕狰,那就相當(dāng)于問題解決。不過控漠,其實(shí)仔細(xì)發(fā)現(xiàn)蔓倍,每一次運(yùn)行GRN結(jié)果雖然不一樣,但是差別不會(huì)說特別巨大盐捷,基本上顯著相關(guān)的基因都還是在的偶翅,只是具體的importance存在一些差異,但其實(shí)后面通過過濾碉渡,許多低質(zhì)量的GRN關(guān)系都會(huì)被過濾掉的倒堕,這也就是說為何許多高質(zhì)量文章的結(jié)果都能大概的復(fù)現(xiàn)出來,尤其是核心的結(jié)果還是比較靠譜的爆价。

二、關(guān)于運(yùn)行矩陣的要求

其實(shí)另一個(gè)比較困惑的問題就是媳搪,在運(yùn)行pyscenic或者scenic的時(shí)候铭段,到底用counts矩陣還是標(biāo)準(zhǔn)化處理之后的矩陣呢?這個(gè)帖子剛好提及到了秦爆,我也就順便總結(jié)一下序愚。

如果是單個(gè)數(shù)據(jù),答案肯定是最好用counts矩陣等限,這也是官方建議的爸吮,因?yàn)閟cenic內(nèi)部也會(huì)進(jìn)行相似的數(shù)據(jù)預(yù)處理芬膝,所以最好用最原始的counts矩陣粉楚,保留的生物學(xué)信息也最原始疹吃。
那么就會(huì)涉及到一個(gè)問題,如果我有多個(gè)樣本的數(shù)據(jù)呢甜橱?counts矩陣之間必然會(huì)存在批次效應(yīng)桐早,而這種樣本間的批次其實(shí)對(duì)相關(guān)性的判斷影響還蠻大的(當(dāng)然GRNBoost存在一定的魯棒性)癣缅,這個(gè)時(shí)候我是應(yīng)該繼續(xù)用counts還是用整合之后去除了批次效應(yīng)的表達(dá)矩陣呢?
hh官方的建議還是蠻中肯的 yes or no哄酝,他回答or友存。建議大家不要用整合的數(shù)據(jù)(存在大批次效應(yīng)的數(shù)據(jù))一次性地去做scenic,最好是一個(gè)個(gè)樣本跑陶衅,最后取交集/均值屡立,從源頭解決問題。當(dāng)然搀军,他也承認(rèn)膨俐,如果非要用整合的樣本,適當(dāng)?shù)男?zhǔn)批次效應(yīng)也是可以的(例如Seurat的SCTranform)奕巍,但對(duì)結(jié)果肯定會(huì)有一定的影響吟策。同時(shí),其實(shí)我們其實(shí)可以自己嘗試一下的止,整合樣本的counts/處理后的矩陣/各個(gè)樣本單獨(dú)運(yùn)行的結(jié)果檩坚,三種情況做一個(gè)綜合考慮,才最合理诅福。(其實(shí)基本上匾委,按我個(gè)人的經(jīng)驗(yàn)來說,如果只選top的話氓润,應(yīng)該都大差不差赂乐,因?yàn)轱@著的永遠(yuǎn)是最突出的,但如果要focus細(xì)微差異的基因的話咖气,大家就需要嚴(yán)謹(jǐn)且謹(jǐn)慎對(duì)待自己的分析結(jié)論了挨措。)

總而言之,本身基因網(wǎng)絡(luò)分析GRN的結(jié)果就是挺玄學(xué)的崩溪,受到諸多因素的影響浅役,所以大家可以結(jié)合其他結(jié)論,例如文獻(xiàn)和實(shí)驗(yàn)來確保最終結(jié)果是可靠的伶唯,大膽分析小心驗(yàn)證即可觉既。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個(gè)濱河市,隨后出現(xiàn)的幾起案子瞪讼,更是在濱河造成了極大的恐慌钧椰,老刑警劉巖,帶你破解...
    沈念sama閱讀 222,681評(píng)論 6 517
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件符欠,死亡現(xiàn)場(chǎng)離奇詭異嫡霞,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)背亥,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 95,205評(píng)論 3 399
  • 文/潘曉璐 我一進(jìn)店門秒际,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人狡汉,你說我怎么就攤上這事娄徊。” “怎么了盾戴?”我有些...
    開封第一講書人閱讀 169,421評(píng)論 0 362
  • 文/不壞的土叔 我叫張陵寄锐,是天一觀的道長(zhǎng)。 經(jīng)常有香客問我尖啡,道長(zhǎng)橄仆,這世上最難降的妖魔是什么? 我笑而不...
    開封第一講書人閱讀 60,114評(píng)論 1 300
  • 正文 為了忘掉前任衅斩,我火速辦了婚禮盆顾,結(jié)果婚禮上,老公的妹妹穿的比我還像新娘畏梆。我一直安慰自己您宪,他們只是感情好,可當(dāng)我...
    茶點(diǎn)故事閱讀 69,116評(píng)論 6 398
  • 文/花漫 我一把揭開白布奠涌。 她就那樣靜靜地躺著宪巨,像睡著了一般。 火紅的嫁衣襯著肌膚如雪溜畅。 梳的紋絲不亂的頭發(fā)上捏卓,一...
    開封第一講書人閱讀 52,713評(píng)論 1 312
  • 那天,我揣著相機(jī)與錄音慈格,去河邊找鬼怠晴。 笑死,一個(gè)胖子當(dāng)著我的面吹牛浴捆,可吹牛的內(nèi)容都是我干的蒜田。 我是一名探鬼主播,決...
    沈念sama閱讀 41,170評(píng)論 3 422
  • 文/蒼蘭香墨 我猛地睜開眼汤功,長(zhǎng)吁一口氣:“原來是場(chǎng)噩夢(mèng)啊……” “哼!你這毒婦竟也來了溜哮?” 一聲冷哼從身側(cè)響起滔金,我...
    開封第一講書人閱讀 40,116評(píng)論 0 277
  • 序言:老撾萬榮一對(duì)情侶失蹤色解,失蹤者是張志新(化名)和其女友劉穎,沒想到半個(gè)月后餐茵,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體科阎,經(jīng)...
    沈念sama閱讀 46,651評(píng)論 1 320
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡,尸身上長(zhǎng)有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 38,714評(píng)論 3 342
  • 正文 我和宋清朗相戀三年忿族,在試婚紗的時(shí)候發(fā)現(xiàn)自己被綠了锣笨。 大學(xué)時(shí)的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 40,865評(píng)論 1 353
  • 序言:一個(gè)原本活蹦亂跳的男人離奇死亡道批,死狀恐怖错英,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情隆豹,我是刑警寧澤椭岩,帶...
    沈念sama閱讀 36,527評(píng)論 5 351
  • 正文 年R本政府宣布,位于F島的核電站璃赡,受9級(jí)特大地震影響判哥,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜碉考,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 42,211評(píng)論 3 336
  • 文/蒙蒙 一塌计、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧侯谁,春花似錦锌仅、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 32,699評(píng)論 0 25
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽(yáng)。三九已至嫩痰,卻和暖如春剿吻,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背串纺。 一陣腳步聲響...
    開封第一講書人閱讀 33,814評(píng)論 1 274
  • 我被黑心中介騙來泰國(guó)打工丽旅, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人纺棺。 一個(gè)月前我還...
    沈念sama閱讀 49,299評(píng)論 3 379
  • 正文 我出身青樓榄笙,卻偏偏與公主長(zhǎng)得像,于是被迫代替她去往敵國(guó)和親祷蝌。 傳聞我的和親對(duì)象是個(gè)殘疾皇子茅撞,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 45,870評(píng)論 2 361

推薦閱讀更多精彩內(nèi)容