在R中通过boxplot绘制右删失数据的生命表四分位数

Question

我有一份回顾调查的数据。到调查时间为止没有经历过特定感兴趣事件的个体被归入审查观察类别，其余的则不被审查。如何绘制显示生命表四分位数的正确删除数据的箱线图，同时考虑到删除和未删除的观测结果？

(我感兴趣的变量'fbi‘是一个持续时间变量，因此对于未经审查的obs持续时间是可用的，对于经审查的情况，我已将持续时间'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

Answer 1

使用您建议的虚拟数据和我添加到答案中的虚拟数据，下面的线条将绘制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