当前位置: 首页>>代码示例 >>用法及示例精选 >>正文


R esoph 吸烟、饮酒和食管癌


R语言 esoph 位于 datasets 包(package)。

说明

来自法国 Ille-et-Vilaine 的食管癌 case-control 研究的数据。

用法

esoph

格式

包含 88 种年龄/酒精/烟草组合记录的 DataFrame 。

[,1]"agegp"年龄阶层1 25--34岁
2 35--44
3 45--54
4 55--64
5 65--74
6 75+
[,2]"alcgp"酒精消耗1 0--39克/天
2 40--79
3 80--119
4 120+
[,3]"tobgp"烟草消费1 0-- 9 克/天
2 10--19
3 20--29
4 30+
[,4]"ncases"案件数量
[,5]"ncontrols"控制数量

例子

require(stats)
require(graphics) # for mosaicplot
summary(esoph)
## effects of alcohol, tobacco and interaction, age-adjusted
model1 <- glm(cbind(ncases, ncontrols) ~ agegp + tobgp * alcgp,
              data = esoph, family = binomial())
anova(model1)
## Try a linear effect of alcohol and tobacco
model2 <- glm(cbind(ncases, ncontrols) ~ agegp + unclass(tobgp)
                                         + unclass(alcgp),
              data = esoph, family = binomial())
summary(model2)
## Re-arrange data for a mosaic plot
ttt <- table(esoph$agegp, esoph$alcgp, esoph$tobgp)
o <- with(esoph, order(tobgp, alcgp, agegp))
ttt[ttt == 1] <- esoph$ncases[o]
tt1 <- table(esoph$agegp, esoph$alcgp, esoph$tobgp)
tt1[tt1 == 1] <- esoph$ncontrols[o]
tt <- array(c(ttt, tt1), c(dim(ttt),2),
            c(dimnames(ttt), list(c("Cancer", "control"))))
mosaicplot(tt, main = "esoph data set", color = TRUE)

作者

Thomas Lumley

来源

Breslow, N.E. 和 Day, N.E. (1980) 癌症研究的统计方法。第一卷:Case-Control 研究分析。 IARC 里昂/牛津大学出版社。

相关用法


注:本文由纯净天空筛选整理自R-devel大神的英文原创作品 Smoking, Alcohol and (O)esophageal Cancer。非经特殊声明,原始代码版权归原作者所有,本译文未经允许或授权,请勿转载或复制。