Origin 也能做主成分分析?

最早了解到主成分分析(Principal Component Analysis熬词,PCA)是在一位師兄的畢業(yè)答辯上旁钧,當(dāng)時聽得是云里霧里,一頭霧水互拾。其實歪今,主成分分析是因子分析的一種常用方法,主要目的是為了減少變量數(shù)目颜矿,也就是降維寄猩。

在高通量測序中,主要基于基因表達(dá)量骑疆、種群豐度等進(jìn)行樣本的聚類田篇,下圖是一篇客戶文章的基于表達(dá)量的PCA結(jié)果替废。

(Plant biotechnology journal, 2018)


那么該如何進(jìn)行主成分分析呢?今天為大家介紹如何用Origin進(jìn)行分析并完成散點圖的繪制泊柬。

數(shù)據(jù)準(zhǔn)備

本文所用的數(shù)據(jù)是OmicShare PCA分析工具的示例文件(如下圖)椎镣,大家可隨意下載練習(xí)(http://www.omicshare.com/tools/Home/Soft/pca)。


示例為6個樣本的1300多個基因的表達(dá)量數(shù)據(jù)兽赁,接下來需要做得是把每個樣本的對應(yīng)的1300多個數(shù)據(jù)用2個數(shù)據(jù)表示状答,降到2維后,就可用二維的散點圖展示刀崖。

這里每個基因的名稱視作1個變量惊科,因此需要對數(shù)據(jù)進(jìn)行手動轉(zhuǎn)置,為了分組著色亮钦,在前面加1列分組數(shù)據(jù)馆截,如下圖(Sheet1 →Sheet3)。然后打開Origin或悲,將數(shù)據(jù)粘貼到Origin的表格中孙咪。

對于這一點OmicShare的 PCA工具就顯得智能和方便多了,直接上傳數(shù)據(jù)巡语,點提交任務(wù)就好啦翎蹈。

分析作圖

點窗口右側(cè)的Apps,在彈出的窗口中點今天要用的PCA插件男公。

然后選擇輸入數(shù)據(jù)(Input data)荤堪,具體方法如下圖,注意這里的數(shù)據(jù)添加小技巧枢赔。

在繪圖設(shè)置中這里選2D的散點圖澄阳,如下圖,當(dāng)然你也可以選3D的散點圖(使用3個主成分)踏拜,點OK 即可完成分析和作圖碎赢。如果樣本較多,也可以選擇顯示置信區(qū)間橢圓(Show Confidence Ellipse)速梗。

雙擊Score Plot中的圖表肮塞,放大散點圖,然后就可以對圖表進(jìn)行個性化調(diào)整姻锁,比如改變標(biāo)記的形狀枕赵、顏色、大小位隶,顯示數(shù)據(jù)標(biāo)簽等拷窜。

最終的的效果如下:


APP安裝

Origin的強(qiáng)大之處在于具有很多官方的拓展工具(Apps),安裝方法與之前介紹的Google Map Import、Heat Map Dendrogram等一樣篮昧,只要將.opx的文件拖拽到打開的Origin界面即可完成安裝赋荆,注意最好安裝Origin 2017之后的版本。


本文用到的PCAC插件和示例數(shù)據(jù)已經(jīng)上傳到論壇恋谭,點擊?這里?可以下載糠睡。


關(guān)于PCA分析挽鞠,更詳細(xì)的教程請查閱周老師的帖子《Omicshare tools——二維PCA分析使用教程》疚颊。


今天的內(nèi)容就到這里啦~

拓展閱讀

《群體結(jié)構(gòu)圖形三劍客——PCA圖》

《OmicShare再添新工具 高級三維PCA圖》

參考文獻(xiàn)

Guo J, Qi J, He K, et al. The Asian cornborer Ostrinia furnacalis feeding increases the direct and indirect defense ofmid‐whorl stage commercial maize in the field[J]. Plant biotechnologyjournal, 2018.

為了更好的閱讀體驗,建議關(guān)注基迪奧生物公眾號~

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末信认,一起剝皮案震驚了整個濱河市材义,隨后出現(xiàn)的幾起案子,更是在濱河造成了極大的恐慌嫁赏,老刑警劉巖其掂,帶你破解...
    沈念sama閱讀 206,013評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件,死亡現(xiàn)場離奇詭異潦蝇,居然都是意外死亡款熬,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,205評論 2 382
  • 文/潘曉璐 我一進(jìn)店門攘乒,熙熙樓的掌柜王于貴愁眉苦臉地迎上來贤牛,“玉大人,你說我怎么就攤上這事则酝⊙臭ぃ” “怎么了?”我有些...
    開封第一講書人閱讀 152,370評論 0 342
  • 文/不壞的土叔 我叫張陵沽讹,是天一觀的道長般卑。 經(jīng)常有香客問我,道長爽雄,這世上最難降的妖魔是什么蝠检? 我笑而不...
    開封第一講書人閱讀 55,168評論 1 278
  • 正文 為了忘掉前任,我火速辦了婚禮挚瘟,結(jié)果婚禮上叹谁,老公的妹妹穿的比我還像新娘。我一直安慰自己刽沾,他們只是感情好本慕,可當(dāng)我...
    茶點故事閱讀 64,153評論 5 371
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著侧漓,像睡著了一般锅尘。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 48,954評論 1 283
  • 那天藤违,我揣著相機(jī)與錄音浪腐,去河邊找鬼。 笑死顿乒,一個胖子當(dāng)著我的面吹牛议街,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播璧榄,決...
    沈念sama閱讀 38,271評論 3 399
  • 文/蒼蘭香墨 我猛地睜開眼特漩,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了骨杂?” 一聲冷哼從身側(cè)響起涂身,我...
    開封第一講書人閱讀 36,916評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎搓蚪,沒想到半個月后蛤售,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,382評論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡妒潭,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點故事閱讀 35,877評論 2 323
  • 正文 我和宋清朗相戀三年悴能,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片雳灾。...
    茶點故事閱讀 37,989評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡漠酿,死狀恐怖,靈堂內(nèi)的尸體忽然破棺而出佑女,到底是詐尸還是另有隱情记靡,我是刑警寧澤,帶...
    沈念sama閱讀 33,624評論 4 322
  • 正文 年R本政府宣布团驱,位于F島的核電站摸吠,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏嚎花。R本人自食惡果不足惜寸痢,卻給世界環(huán)境...
    茶點故事閱讀 39,209評論 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一處隱蔽的房頂上張望紊选。 院中可真熱鬧啼止,春花似錦、人聲如沸兵罢。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,199評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽卖词。三九已至巩那,卻和暖如春,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背即横。 一陣腳步聲響...
    開封第一講書人閱讀 31,418評論 1 260
  • 我被黑心中介騙來泰國打工噪生, 沒想到剛下飛機(jī)就差點兒被人妖公主榨干…… 1. 我叫王不留,地道東北人东囚。 一個月前我還...
    沈念sama閱讀 45,401評論 2 352
  • 正文 我出身青樓跺嗽,卻偏偏與公主長得像,于是被迫代替她去往敵國和親页藻。 傳聞我的和親對象是個殘疾皇子桨嫁,可洞房花燭夜當(dāng)晚...
    茶點故事閱讀 42,700評論 2 345

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