Чтобы нарисовать квартили таблицы смертности через коробчатую диаграмму для цензурированных справа данных в R

У меня есть данные ретроспективного опроса. Лица, которые не пережили конкретное событие, представляющее интерес до момента проведения исследования, помещаются в категорию цензурированного наблюдения, а остальные не подвергаются цензуре. Как нарисовать коробчатую диаграмму для этих подвергнутых цензуре справа данных, которые показывают квартили таблицы смертности с учетом как цензурированных, так и нецензурных наблюдений?

(Моя интересующая переменная '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

r boxplot survival-analysis
person Shruti    schedule 02.02.2013    source источник
comment
Не могли бы вы предоставить воспроизводимые образцы данных? - stackoverflow.com / questions / 5963269 /   -  person thelatemail    schedule 02.02.2013
comment
Извините, я здесь новичок ... Я не понимаю, как правильно вводить код здесь, в текстовом поле ........ пожалуйста, помогите мне!   -  person Shruti    schedule 02.02.2013
comment
но я могу сказать вам, что данные имитируют это: ›fbi‹ -rpois (100,12) ›cens‹ -sample (0: 1,100, replace = T) ›data‹ -data.frame (a, b), здесь в fbi , наблюдения, соответствующие cens = 0, - это те, для которых длительность fbi не сообщалась, поскольку они все еще продолжались, поэтому я заменил fbi, соответствующий нулю, на продолжительность от начала до даты исследования. (!! извините за плохо написанный код)   -  person Shruti    schedule 02.02.2013
comment
Извините!! это data.frame (fbi, cens)   -  person Shruti    schedule 02.02.2013
comment
Добро пожаловать в SO. Я отредактировал вопрос, включив в него данные, которые вы предоставили в комментариях выше. В следующий раз просто добавьте такую ​​(очень полезную!) Информацию прямо в вопрос, чтобы потенциальные респонденты знали, как выглядят ваши данные.   -  person thelatemail    schedule 02.02.2013

Ответы (1)


arrow_upward
0
arrow_downward

Используя фиктивные данные, которые вы предложили и которые я добавил к ответу, в строке ниже будут построены 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.

person thelatemail    schedule 02.02.2013
comment
Большое спасибо за ответ! но я ищу цензурированные коробчатые диаграммы ... то есть ящичные диаграммы, которые показывают квартили таблицы смертности. Я слышал о пакете NADA, но не знаю, будет ли он работать здесь или нет для правильных подвергнутых цензуре данных. - person Shruti; 02.02.2013
comment
еще кое-что! не могли бы вы объяснить мне пошагово, как вы писали код здесь, в текстовом поле? Что мне делать, если я хочу скопировать свой код из сценария R в поле для комментариев здесь? - person Shruti; 02.02.2013
comment
@Shruti - хорошо, я не очень знаком с такими вещами, но отредактировал свой ответ, включив функцию блочной диаграммы с использованием библиотеки NADA. - person thelatemail; 02.02.2013
comment
@Shruti - если вы вставите свой код в текстовое поле, затем выделите код и нажмите кнопку, которая выглядит как {}, вы получите отформатированный код. - person thelatemail; 02.02.2013
comment
Спасибо! но я не могу найти кнопку {}. - person Shruti; 02.02.2013