我如何告诉 R 包 Limma 在 read.idat() 中使用什么作为 "targets"？

Question

我正在分析一些微阵列数据。对于每个捐赠者，我都有一个“干预前”和“干预后”的 idat 文件。我已经使用带有 read.idat() 函数的 Limma 包成功地将这些读入 R。但是，生成的对象在目标中只有一列：“IDATfile”。我相信，如果我使用 read.ilmn()，我会指定一个 targets.txt 文件，但在使用 read.idat() 时我看不到这个选项。例如。在 Limma user guide Illumina 示例中，目标是“Donor”、“Age”和“Cell Type”。 我如何告诉 Limma 将什么作为目标？我想要“捐助者”和“干预”。

我的意思的一个例子：

idatfiles <- dir(pattern="idat")
bgxfile <- dir(pattern="bgx")
x <- read.idat(idatfiles, bgxfile)
colnames(x$targets)
[1] "IDATfile"

我希望这是“捐赠者”和“干预”，而不是“IDATfile”。我可以通过 read.idat(..., dateinfo=TRUE) 将原始 IDAT 文件的其他一些列作为进一步的目标，但我不知道如何编辑这些列以使它们成为“Donor”和“干预”：

[1] "IDATfile"   "ScanInfo"   "DecodeInfo"

如果需要更多信息，请告诉我，非常感谢您的帮助！

Answer 1

如果您只想编辑别名，您可以使用：

colnames(x$targets) <- c("Donor")

但我认为行是目标数据框中的样本，所以这真的是您想要的吗？

http://web.mit.edu/~r/current/arch/i386_linux26/lib/R/library/limma/html/EList.html

targets data.frame containing information on the target RNA samples. Rows correspond to samples. May have any number of columns.

我如何告诉 R 包 Limma 在 read.idat() 中使用什么作为 "targets"？

How do I tell R package Limma what to use as "targets" in read.idat()?

r

bioinformatics

data-analysis

bioconductor

limma