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

Question

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

Answer 1

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

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

Answer 2

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

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

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