我有来自回顾性调查的数据。没有经历特定事件的个人在调查时间之前被放入审查观察类别，其余未经审查。如何绘制这个正确的审查数据的boxplot，显示生命表四分位数同时考虑审查和未审查的观察？我的利益变量'fbi'是一个持续时间变量，因此对于未经审查的obs持续时间是可用的并且对于审查，我已经用原始到调查日期之间的时间间隔和另一个二分变量“cens'取代了持续时间'fbi' “有没有认识到审查和未经审查的情况下）要通过boxplot绘制生命表四分位数R

的数据可以与模拟：

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

2013-02-02 Shruti

你能提供的可重复的一些示例数据？ - http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example – thelatemail

对不起，我在这里真的很新..我没有得到如何正确的代码在这里的文本框。 ....... 请帮助我！ – Shruti

但我可以告诉你数据模拟这个：> fbi <-rpois（100,12）> cens <-sample（0：1,100，replace = T）> data <-data.frame（a，b），这里fbi，与cens = 0相对应的观察值是那些持续时间fbi由于它们仍在继续而没有被报告的，所以我已经用从起始到调查日期的持续时间代替了对应于零的fbi。（!!对不起，严重的代码） – Shruti

使用你的建议将假数据，并且我加入到了答案，下面的线将绘制2个箱线图总结包括所有情况的fbi变量，并仅使用非审查d案件。

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

来源

2013-02-02 06:13:25 thelatemail

非常感谢您的回复！ – Shruti

非常感谢您的回复！但我正在寻找经过审查的箱型图......即显示生命表四分位数的箱型图。我听说过有关“NADA”的软件包，但不知道它是否适用于正确的审查数据。 – Shruti

还有一件事！你能否请你逐步解释我是如何在文本框中编写代码的？如果我想在这里的评论框中复制R脚本中的代码，我该怎么办？ – Shruti

要通过boxplot绘制生命表四分位数R

回答

编辑

相关问题