我有这样的文件
Age.Range Average Probability
1 0 to 04 400 0.00400
2 05 to 09 221 0.00221
3 10 to 14 216 0.00216
4 15 to 19 409 0.00409
X [个人年龄; 0至80年之间的整数]
Y [监视个人的持续时间; 1至之间的整数 50年或“终身”]我需要计算X年龄(例3)的人在从今天开始直到Y(例7)的时间间隔内患上癌症的可能性。在R中,我需要在第一列中找到X值和X + Y值,并对这两个范围之间的第三列中的所有值求和:
X= 3
x+y=10
probability= 0.004 + 0.00221 + 0.00216
答案 0 :(得分:1)
以下功能可满足您的需求。它获取年龄范围的起点,然后使用findInterval
在“概率”列中找到索引。然后,只需添加这些概率即可。
sumProbs <- function(DF, X, Y){
DF[["Age.Range"]] <- as.character(DF[["Age.Range"]])
Age.Start <- strsplit(DF[["Age.Range"]], " to ")
Age.Start <- as.integer(sapply(Age.Start, '[[', 1))
i <- findInterval(c(X, X + Y), Age.Start)
p <- DF[["Probability"]][i[1]:i[2]]
sum(p)
}
sumProbs(df1, 3, 7)
#[1] 0.00837
dput
格式的数据。
df1 <-
structure(list(Age.Range = c("0 to 04", "05 to 09",
"10 to 14", "15 to 19"), Average = c(400L, 221L,
216L, 409L), Probability = c(0.004, 0.00221, 0.00216,
0.00409)), row.names = c("1", "2", "3", "4"),
class = "data.frame")