通过箱图绘制生命表四分位数,用于R中的右删失数据

时间:2013-02-02 03:55:33

标签: r boxplot survival-analysis

我有回顾性调查的数据。在调查时间之前没有经历过特定感兴趣事件的个人被置于审查观察类别中,其余未经审查。如何绘制这个权利审查数据的箱线图,显示生命表四分位数,同时考虑到审查和未审查的观察结果?

(我感兴趣的变量' fbi'是一个持续时间变量,因此对于未经审查的障碍持续时间是可用的,对于被删除的我已经用原点之间的时间间隔替换了持续时间' fbi'调查日期和另一个二分变量" cens"是否可以识别审查和未经审查的案件。)

可以使用以下方式模拟数据:

fbi <- rpois(100,12)
cens <- sample(0:1,100,replace=T) 
test <- data.frame(fbi,cens)

> head(test)
  fbi cens
1  18    0
2  14    0
3  17    1
4  11    1
5   9    0
6  10    1

1 个答案:

答案 0 :(得分:0)

使用您建议的虚拟数据以及我在答案中添加的数据,下面的行将绘制2个箱形图,总结fbi变量,包括所有案例,并仅使用非审查案例。

boxplot(test$fbi,test$fbi[test$cens==0],names=c("all cases","w/out censored"))

如果你想比较审查和未审查的案件,你可以这样做:

boxplot(fbi ~ cens,data=test,names=c("not censored","censored"))

修改

在回复下面的评论时,您正在寻找使用NADA库的以下代码吗?

library(NADA)
cenboxplot(test$fbi, as.logical(test$cen))

此处有cenboxplot在线功能的文档:http://rss.acs.unt.edu/Rdoc/library/NADA/html/cenboxplot.html