我有來自回顧性調查的數據。沒有經歷特定事件的個人在調查時間之前被放入審查觀察類別，其餘未經審查。如何繪製這個正確的審查數據的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

回答

編輯

相關問題