0

遡及調査のデータがあります。調査時間までに関心のある特定のイベントを経験していない個人は、打ち切られた観察カテゴリに入れられ、残りは打ち切られていません。打ち切られた観測と打ち切られていない観測の両方を考慮して、生命表の四分位数を示すこの正しい打ち切られたデータの箱ひげ図を描く方法は?

(私の興味のある変数「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
4

1 に答える 1

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

于 2013-02-02T06:13:25.550 に答える