1

我正在使用一些使用数据表类而不是数据框类的 R 代码。

我将如何在 R 中执行以下操作而无需将 map.dt 转换为 map.df?

map.dt = data.table(chr = c("chr1","chr1","chr1","chr2"), ref = c(1,0,3200,3641), pat = c(1,3020,3022, 3642), mat = c(1,0,3021,0))
parent = "mat"
chrom = "chr1"
map.df<-as.data.frame(map.dt);
parent.block.starts<-map.df[map.df$chr == chrom & map.df[,parent] > 0,parent];

注意:父母需要动态分配,它是用户的输入。在此示例中,我选择了“垫子”,但它可以是任何列。

Note1: parent.block.starts 应该是一个整数向量。

注 2:map.dt 是一个数据表,其中列名是 c("chr","re​​f","pat","mat")。

问题是在数据表中我无法按名称访问给定列,或者至少我不知道如何访问。

如果您有任何建议,请告诉我!

谢谢!

4

1 回答 1

3

这里的最终目标有点不清楚,特别是没有样本数据,但如果你想按字符名称访问行,有两种方法可以做到这一点:

Columns = c("A", "B")
# .. means "look up one level"
dt[,..Columns]

dt[,get("A")]
dt[,list(get("A"), get("B"))]

但是,如果您发现自己需要经常使用这种技术,那么您可能使用data.table得不好。

编辑

根据您的编辑,此行将返回相同的结果,而无需进行任何as.data.frame转换:

> map.dt[chr==chrom & get(parent) > 0, get(parent)]
于 2013-09-20T14:52:31.930 回答