0

我正在做安捷伦双通道微阵列分析的差异表达，我得到了 MA 列表，我想将其转换为表达矩阵，我可以使用 MAlist$A 来生成它，因为它是 log2 值。M矩阵有很多负值，真的很低，不知道好不好？这是箱线图 在此处输入图像描述

我的代码

  require('limma')
    
    targetinfo <- readTargets('targets.txt', sep = '\t')
    project <- read.maimages(targetinfo, source = 'agilent')
    # Perform background correction on the fluorescent intensities
   MA <- backgroundCorrect(project, method = 'normexp', offset = 16)

这是我的目标.txt

show(targets)
                    FileName disease
1  SampleFiles/GSM176570.txt   FX     
2  SampleFiles/GSM176575.txt   FX     
3  SampleFiles/GSM176579.txt   FX     
4  SampleFiles/GSM176583.txt   FX     
5  SampleFiles/GSM176586.txt   C     
6  SampleFiles/GSM176589.txt   C
4

1 回答 1

1

简短的回答是您应该使用 M 值。我看不到你的目标文件，所以我不知道背景校正完成了什么。

转换的一个好方法是使用本手册中convert显示的包作为表达式集，所以下面我使用了一个与你类似的示例数据：

library(limma)
targetinfo <- readTargets("TargetBeta7.txt")
project <- read.maimages(targetinfo,source="genepix")
project <- backgroundCorrect(project, method="normexp", offset=50)
MA <- normalizeWithinArrays(project)

然后转换：

#install using BiocManager::install("convert")
library(convert)
eset = as(MA, "ExpressionSet")

您可以检查这些值：

head(exprs(eset))
    6Hs.195.1    6Hs.168     6Hs.166    6Hs.187.1     6Hs.194   6Hs.243.1
1 -0.14680077  0.1104113 -0.01172555 -0.057836495 -0.06372548  0.85761577
2 -0.08776962  0.1893185 -0.02600634 -0.072301864  0.19329923  1.39021146
3  0.55834060  0.3051494  0.33223516 -0.008061431  0.46835596  0.41179532
4  0.21925587  0.4453152 -0.10680881 -0.026872543  0.09459589 -0.01676835
5 -0.33695592 -0.2091978 -0.13209201 -0.093929553 -0.16276967 -0.19276659
6 -0.46437638 -0.2647629 -1.04651233 -0.067690463 -0.20793054 -0.39403884

它与以下内容相同：

head(MA$M)
      6Hs.195.1    6Hs.168     6Hs.166    6Hs.187.1     6Hs.194   6Hs.243.1
[1,] -0.14680077  0.1104113 -0.01172555 -0.057836495 -0.06372548  0.85761577
[2,] -0.08776962  0.1893185 -0.02600634 -0.072301864  0.19329923  1.39021146
[3,]  0.55834060  0.3051494  0.33223516 -0.008061431  0.46835596  0.41179532
[4,]  0.21925587  0.4453152 -0.10680881 -0.026872543  0.09459589 -0.01676835
[5,] -0.33695592 -0.2091978 -0.13209201 -0.093929553 -0.16276967 -0.19276659
[6,] -0.46437638 -0.2647629 -1.04651233 -0.067690463 -0.20793054 -0.39403884
于 2021-02-20T11:18:39.763 回答