Двухсторонняя цензурная модель в R (похожая на Zeligs Tobit)?

Question

Существует ли модель для зависимых переменных, которые подвергаются цензуре с обеих сторон? И если так, есть ли реализация в R? Я знаю только о моделях тобитов (например, в пакете Zelig), но они, очевидно, подвергаются цензуре только на левой стороне... Интересно, имеет ли смысл урезать обе стороны...

Answer 1

1. Есть разница между усечением и цензурой. Вы должны знать, в чем дело, прежде чем начинать моделирование. (в двух словах: цензура означает, что события могут быть обнаружены, но измерения не известны полностью (т. е. в вашем случае вы не знаете ни точное начало, ни точное окончание временного интервала, в котором субъекты подвергались риску для рассматриваемого события) Обрезание означает, что события могут наблюдаться только в том случае, если выполнено еще одно условие: популярным примером является выживание в доме престарелых, в котором могут принимать на работу только люди старше 65 лет - вхождение в исследуемую группу населения затем усекается в возрасте 65 лет.)

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

Answer 2

В очень реальном смысле большинство наблюдательных исследований популяций "людей на свободном расстоянии" подвергаются двойной цензуре... т.е. мы не наблюдаем людей в течение всей их жизни. Вот цитата к докторской диссертации, которая, кажется, хорошо излагает статистическую терминологию. Кроме того, некоторые из пакетов в R будут работать должным образом при настройке на интервальную цензуру или левую цензуру, включая пакеты на выживание, NADA, sand (с их веб-сайта DOE) и некоторые другие, по которым вы можете искать на веб-сайте Baron с помощью соответствующих стратегий поиска. по этой ссылке, которая настраивает эту страницу для получения как функций, так и записей r-help.

Изменить: добавив комментарии, чтобы прояснить, что речь идет об усечении, а не о цензуре.

Если вы хотите соответствовать усеченным дистрибутивам, посмотрите на пакет gamlss или создайте подходящую плотность для дважды усеченных дистрибутивов и используйте fitdistr в пакете MASS.