【包】SummarizedExperiment

SummarizedExperiment類

該類常被用來儲存表達矩陣向楼,比如TCGAbiolinks下載的TCGA表達矩陣。
該類數(shù)據(jù)將原始信息(metadata)和表達矩陣(assays)聯(lián)合到了一起谐区,可同步對他們取子集湖蜕。
其實和bioconductor中的ExpreesionSet很接近,只是行信息比較靈活宋列,允許使用GRanges或DataFrame對象昭抒。

1、基本結構

如下圖所示,每個對象可以存儲多個表達矩陣(assays)灭返,行一般是基因或引物盗迟,列一般是樣本編號。每個表達矩陣都有兩個關聯(lián)的表(Row data, Column data)婆殿,分別通過基因名和樣本名與表達矩陣相關聯(lián)诈乒,存儲著基因以及樣本的相關信息。此外婆芦,還有單獨的一個文件(meta data)描述整個實驗的相關信息怕磨。


SummarizedExperiment

2、操作匯總

代碼僅展示了格式

test<-SummarizedExperiment(assays= list(counts), rowData= DataFrame(), 
colData= DataFrame(Treatment= rep(c("flame","non_flame"),3)))
# 創(chuàng)建一個對象消约,有無rowData都行
assays(test) # 取出表達矩陣肠鲫,取出的是元素為列表的列表
assay(test,1) # 取出的是列表
#assays(test)[[1]][1:5,1:5] == assay(test)[1:5,1:5]==assay(test,1)[1:5,1:5]
assays(test)$counts # 取出特定表達矩陣
rowData(test) # 取出行信息
colData(test) # 取出列信息
metadata(test) # 取出實驗信息
test[,]可以像數(shù)據(jù)框一樣取子集
test[,test$Treatment=="flame"] # 可用這種形式根據(jù)列信息對表達矩陣取子集,同樣可以用在行信息取子集

RangedSummarizedExperiment類

該類是SummarizeddExperiment的子類或粮,繼承了其所有的方法导饲。
主要區(qū)別在于:該類的行是基因范圍,其關聯(lián)表是GRanges類或GRangesList類氯材。

rowRangs() # 取得行信息
rang<-GRanges(seqnames="1", ranges=100000:1100000)
subsetByOverlaps(test,rang)
#使用Granges類的方法取chromosom 1的100000:1100000的行
最后編輯于
?著作權歸作者所有,轉載或內容合作請聯(lián)系作者
  • 序言:七十年代末渣锦,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子氢哮,更是在濱河造成了極大的恐慌袋毙,老刑警劉巖,帶你破解...
    沈念sama閱讀 206,839評論 6 482
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件冗尤,死亡現(xiàn)場離奇詭異听盖,居然都是意外死亡,警方通過查閱死者的電腦和手機裂七,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,543評論 2 382
  • 文/潘曉璐 我一進店門皆看,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人背零,你說我怎么就攤上這事腰吟。” “怎么了徙瓶?”我有些...
    開封第一講書人閱讀 153,116評論 0 344
  • 文/不壞的土叔 我叫張陵蝎困,是天一觀的道長。 經常有香客問我倍啥,道長,這世上最難降的妖魔是什么澎埠? 我笑而不...
    開封第一講書人閱讀 55,371評論 1 279
  • 正文 為了忘掉前任虽缕,我火速辦了婚禮,結果婚禮上蒲稳,老公的妹妹穿的比我還像新娘氮趋。我一直安慰自己伍派,他們只是感情好,可當我...
    茶點故事閱讀 64,384評論 5 374
  • 文/花漫 我一把揭開白布剩胁。 她就那樣靜靜地躺著诉植,像睡著了一般。 火紅的嫁衣襯著肌膚如雪昵观。 梳的紋絲不亂的頭發(fā)上晾腔,一...
    開封第一講書人閱讀 49,111評論 1 285
  • 那天,我揣著相機與錄音啊犬,去河邊找鬼灼擂。 笑死,一個胖子當著我的面吹牛觉至,可吹牛的內容都是我干的剔应。 我是一名探鬼主播,決...
    沈念sama閱讀 38,416評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼语御,長吁一口氣:“原來是場噩夢啊……” “哼峻贮!你這毒婦竟也來了?” 一聲冷哼從身側響起应闯,我...
    開封第一講書人閱讀 37,053評論 0 259
  • 序言:老撾萬榮一對情侶失蹤纤控,失蹤者是張志新(化名)和其女友劉穎,沒想到半個月后孽锥,有當?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體嚼黔,經...
    沈念sama閱讀 43,558評論 1 300
  • 正文 獨居荒郊野嶺守林人離奇死亡,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內容為張勛視角 年9月15日...
    茶點故事閱讀 36,007評論 2 325
  • 正文 我和宋清朗相戀三年惜辑,在試婚紗的時候發(fā)現(xiàn)自己被綠了唬涧。 大學時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點故事閱讀 38,117評論 1 334
  • 序言:一個原本活蹦亂跳的男人離奇死亡盛撑,死狀恐怖碎节,靈堂內的尸體忽然破棺而出,到底是詐尸還是另有隱情抵卫,我是刑警寧澤狮荔,帶...
    沈念sama閱讀 33,756評論 4 324
  • 正文 年R本政府宣布,位于F島的核電站介粘,受9級特大地震影響殖氏,放射性物質發(fā)生泄漏。R本人自食惡果不足惜姻采,卻給世界環(huán)境...
    茶點故事閱讀 39,324評論 3 307
  • 文/蒙蒙 一雅采、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧,春花似錦婚瓜、人聲如沸宝鼓。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,315評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽愚铡。三九已至,卻和暖如春胡陪,著一層夾襖步出監(jiān)牢的瞬間沥寥,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 31,539評論 1 262
  • 我被黑心中介騙來泰國打工督弓, 沒想到剛下飛機就差點兒被人妖公主榨干…… 1. 我叫王不留营曼,地道東北人。 一個月前我還...
    沈念sama閱讀 45,578評論 2 355
  • 正文 我出身青樓愚隧,卻偏偏與公主長得像蒂阱,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子狂塘,可洞房花燭夜當晚...
    茶點故事閱讀 42,877評論 2 345

推薦閱讀更多精彩內容