#r #survival-analysis #exponential #reliability #weibull
#r #анализ выживаемости #экспоненциальный #надежность #вейбулл
Вопрос:
Я хочу сгенерировать 10%, 20% и 30% подвергнутые цензуре данные в R , используя следующие два распределения.
# generate the life-time observations
data1 = c(rweibull(15, 1,2), rweibull(35, 2,3))
# censoring time observations
data2 = rexp(50, c)
Но я хочу определить c по заданной пропорции цензуры, например, 10%, 20% и 30% . То есть мы можем вычислять по P(data1 >= data2) = 10% , P(data1 >= data2) = 20% и так далее. Мой вопрос в том, как найти c для заданной пропорции? Есть ли у нас какие-либо конкретные r функции? Спасибо!
Комментарии:
1. В статистических терминах вы запрашиваете аналитические решения для пересечения кумулятивных распределений смеси Вейбуллов и экспоненциального. Кажется крайне маловероятным, что может быть создано простое решение, но есть пакет, который позволяет выполнять высокоуровневые операции с дистрибутивами. Установите
distrпакет, а затем обработайте виньетку. Он маркирует себя как «Объектно-ориентированную реализацию распределений», и это отличный способ получить опыт работы со статистическими распределениями. Это похоже на мгновенный курс математической статистики верхнего уровня.2. Вы искали другие мнения? Я нашел интересное обсуждение экспоненциальных распределений с преобразованием мощности. Оказывается, они Weibull. Вы могли бы применить преобразование мощности к экспоненте, чтобы сделать ее квантиль 0.7 таким же, как у вашей случайной реализации несбалансированной смеси Вейбуллов. actuarialmodelingtopics.wordpress.com/2017/06/10 /…