affy包查看全匹配探针与错配探针

affymetrix表达谱芯片 查看全匹配探针PM与错配探针MM

affymetrix表达谱芯片在设计时 除了全匹配(perfect match, pm) 探针(25-mer探针)进行进行荧光信号检测外,有时会设计错配(mismatch, mm) 探针(一个位点错配),基于PM-MM探针设计可以使检测低浓度靶序列的特异性和灵敏度大大提高,也便于进行背景校正等。而由于芯片原始数据处理的方法不同,有时会利用错配探针信号,而有时并不会参考错配探针信息。譬如,常见的 RMA算法进行数据标准化并未利用到mm信号值,mas5 方法则要参考mm进行计算

因此在实际分析的过程中需要明确对应的affymetrix芯片 是不是包含两种探针信号,可选择的原始数据标准化的分析方法是不是相互配套。

分析过程中可以对读取的affymetrix芯片原始数据,所对应的AffyBatch 对象 利用pm() 函数和mm()查看是否包含对应的两种探针信号值

以 GSE66196 数据为例  https://www.omicsclass.com/article/609 

> PM=pm(AffyBatch)
> head(PM,1)
       GSM1616746_Pool_1_Hilscher_26-07-2012.CEL.gz GSM1616747_Pool_1_Reif-Hilscher_8-2-2012.CEL.gz
369707                                          390                                             205
       GSM1616748_Pool_1wdh_Reif-Hilscher_9-2-2012.CEL.gz GSM1616749_Pool_2_Hilscher_26-07-2012.CEL.gz
369707                                                181                                          497
       GSM1616750_Pool_2_Reif-Hilscher_8-2-2012.CEL.gz GSM1616751_Pool_2wdh_Reif-Hilscher_9-2-2012.CEL.gz
369707                                             223                                                190
       GSM1616752_Pool_3_Hilscher_26-07-2012.CEL.gz GSM1616753_Pool_3_Reif-Hilscher_8-2-2012.CEL.gz
369707                                          525                                             240
       GSM1616754_Pool_3wdh_Reif-Hilscher_9-2-2012.CEL.gz GSM1616755_Pool_4_Hilscher_26-07-2012.CEL.gz
369707                                                192                                          446
       GSM1616756_Pool_4_Reif-Hilscher_8-2-2012.CEL.gz GSM1616757_Pool_4wdh_Reif-Hilscher_9-2-2012.CEL.gz
369707                                             214                                                155
> MM=mm(AffyBatch)
> head(MM,1)
       GSM1616746_Pool_1_Hilscher_26-07-2012.CEL.gz GSM1616747_Pool_1_Reif-Hilscher_8-2-2012.CEL.gz
370871                                          174                                             157
       GSM1616748_Pool_1wdh_Reif-Hilscher_9-2-2012.CEL.gz GSM1616749_Pool_2_Hilscher_26-07-2012.CEL.gz
370871                                                159                                          215
       GSM1616750_Pool_2_Reif-Hilscher_8-2-2012.CEL.gz GSM1616751_Pool_2wdh_Reif-Hilscher_9-2-2012.CEL.gz
370871                                             148                                                137
       GSM1616752_Pool_3_Hilscher_26-07-2012.CEL.gz GSM1616753_Pool_3_Reif-Hilscher_8-2-2012.CEL.gz
370871                                          213                                             180
       GSM1616754_Pool_3wdh_Reif-Hilscher_9-2-2012.CEL.gz GSM1616755_Pool_4_Hilscher_26-07-2012.CEL.gz
370871                                                152                                          234
       GSM1616756_Pool_4_Reif-Hilscher_8-2-2012.CEL.gz GSM1616757_Pool_4wdh_Reif-Hilscher_9-2-2012.CEL.gz
370871                                             151                                                134





相关课程:GEO芯片数据挖掘GEO芯片数据标准化


更多生物信息课程:

1. 文章越来越难发?是你没发现新思路,基因家族分析发2-4分文章简单快速,学习链接:基因家族分析实操课程基因家族文献思路解读

2. 转录组数据理解不深入?图表看不懂?点击链接学习深入解读数据结果文件,学习链接:转录组(有参)结果解读转录组(无参)结果解读

3. 转录组数据深入挖掘技能-WGCNA,提升你的文章档次,学习链接:WGCNA-加权基因共表达网络分析

4. 转录组数据怎么挖掘?学习链接:转录组标准分析后的数据挖掘转录组文献解读

5. 微生物16S/ITS/18S分析原理及结果解读OTU网络图绘制cytoscape与网络图绘制课程

6. 生物信息入门到精通必修基础课:linux系统使用perl入门到精通perl语言高级R语言入门R语言画图

7. 医学相关数据挖掘课程,不用做实验也能发文章:TCGA-差异基因分析GEO芯片数据挖掘GEO芯片数据标准化GSEA富集分析课程TCGA临床数据生存分析TCGA-转录因子分析TCGA-ceRNA调控网络分析

8.其他,二代测序转录组数据自主分析NCBI数据上传二代测序数据解读




  • 发表于 2019-02-15 15:57
  • 阅读 ( 3762 )
  • 分类:R

0 条评论

请先 登录 后评论
Daitoue
Daitoue

167 篇文章

作家榜 »

  1. omicsgene 700 文章
  2. 安生水 348 文章
  3. Daitoue 167 文章
  4. 生物女学霸 120 文章
  5. xun 82 文章
  6. 红橙子 78 文章
  7. rzx 75 文章
  8. CORNERSTONE 72 文章