ExpressionSet 数据结构

ExpressionSet 是 Biobase 提供的将多种来源信息组合成方便的单个数据结构,在芯片数据分析很常用。
ExpressionSet 对象结构总结如下图。其中 assayData 保存表达数据;metaData 分别包含 featureData 和 phenoData, 保存探针注释信息和样品信息;experimentData 保存实验信息。


ExpressionSet

当我们用 GEOquery 处理 GEO 数据时可以用 exprs, pData 等函数,就是因为数据被处理成 ExpressionSet 对象了。

使用 ALL 包可以进行 ExpressionSet 对象操作练习,熟悉数据结构和方法函数。

library(Biobase)
library(ALL)

data(ALL)

exprs 函数提取表达矩阵(assayData)。

> exprs(ALL)[1:3, 1:3]
             01005    01010    03002
1000_at   7.597323 7.479445 7.567593
1001_at   5.046194 4.932537 4.799294
1002_f_at 3.900466 4.208155 3.886169

行名为探针名;列名为样品名。可以分别用 featureNamessampleNames 提取。

> featureNames(ALL)[1:5]
[1] "1000_at"   "1001_at"   "1002_f_at" "1003_s_at" "1004_at"  
> sampleNames(ALL)[1:5]
[1] "01005" "01010" "03002" "04006" "04007"

pData 函数取得样品信息。

> pData(ALL)[1:3, 1:8]
       cod diagnosis sex age BT remission CR   date.cr
01005 1005 5/21/1997   M  53 B2        CR CR  8/6/1997
01010 1010 3/29/2000   M  19 B2        CR CR 6/27/2000
03002 3002 6/24/1998   F  52 B4        CR CR 8/17/1998

相应地,用 featureData 函数取得探针注释数据。

> featureData(ALL)
An object of class 'AnnotatedDataFrame': none

(然而,什么都没有发生)
没有注释数据时,可以自己去 GEO 找到相应芯片平台注释文件,或者寻找相应的注释 R 包。

experimentData 提取实验信息。

> experimentData(ALL)
Experiment data
  Experimenter name: Chiaretti et al. 
  Laboratory: Department of Medical Oncology, Dana-Farber Cancer Institute, Department of Medicine, Brigham and Women's Hospital, Harvard Medical School, Boston, MA 02115, USA. 
  Contact information:  
  Title: Gene expression profile of adult T-cell acute lymphocytic leukemia identifies distinct subsets of patients with different response to therapy and survival. 
  URL:  
  PMIDs: 14684422 16243790 

  Abstract: A 187 word abstract is available. Use 'abstract' method.

ExpressionSet 可以进行索引切片。

> ALL[1:3, 1:3]
ExpressionSet (storageMode: lockedEnvironment)
assayData: 3 features, 3 samples 
  element names: exprs 
protocolData: none
phenoData
  sampleNames: 01005 01010 03002
  varLabels: cod diagnosis ... date last seen (21 total)
  varMetadata: labelDescription
featureData: none
experimentData: use 'experimentData(object)'
  pubMedIds: 14684422 16243790 
Annotation: hgu95av2 

参考资料
ExpressionSet
An end to end workflow for differential gene expression using Affymetrix microarrays

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 226,679评论 6 526
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 97,610评论 3 411
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 174,336评论 0 372
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 62,146评论 1 306
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 71,001评论 6 405
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 54,522评论 1 318
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 42,633评论 3 433
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 41,777评论 0 283
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 48,291评论 1 329
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 40,272评论 3 352
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 42,410评论 1 363
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 37,977评论 5 354
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 43,671评论 3 342
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 34,086评论 0 25
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 35,299评论 1 278
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 50,991评论 3 385
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 47,397评论 2 369

推荐阅读更多精彩内容