基于其他两个因素水平的测量值的组合创建一个新的因素水平

Question

我正在分析腹部CT扫描中不同身体部位的自动脂肪体积测量。在每个患者的扫描中，在连续的椎体水平上进行测量，每个患者有多个单独测量的隔室（皮下和内脏）。先前的研究已经确定内脏/皮下脂肪测量的比率特别重要。

我很难在我的数据集中计算这个比率。在这个示例代码中，每个患者有六个条目。每个条目都与椎骨水平的隔间测量的脂肪体积相关联。

我想做的是创建一个新的测量类型-“vat/sat”-这只是三个椎体水平中每一个的两个测量值的比率。本质上，我试图为每个患者插入三个新的观察值，这些观察值与一个新的因子水平和值相关，该因子水平和值是对其他值的操作。非常感谢您的帮助。

library(data.table)
data <- data.table(ID = rep(c(1:4),each = 6), value = rnorm(24, 1000, 500),
                   level = rep(c('l1','l2','l3')), 
                   measure = rep(c(rep('vat',3),rep('sat',3)),4))

编辑：我一直在这个项目中使用data.table，熟悉基本操作，但似乎无法解决这个问题。

Answer 1

我会考虑使用更自然的宽格式：

res = dcast(data, ID + level ~ measure)[, rat := vat/sat][]

回到long

melt(res, id=c("ID", "level"))

由于数据的特殊性，最后需要的是[]。表格打印。如果没有它，当你打字时。。。

> res = dcast(data, ID + level ~ measure)[, rat := vat/sat][]
> res 
# nothing happens
> res 
# now it prints

我不确定它是否在函数留档中，但您可能需要使用BrowseVignettes（"data.table"）查看小插曲，因为它们涵盖了这样的怪癖并有助于建立语法的直觉。