Подмножество категориального ковариата регрессии в R - PullRequest
Новая
       15

Подмножество категориального ковариата регрессии в R

1 голос
/ 27 мая 2019

Я пытаюсь установить подкласс категориального ковариата под названием «Public», который имеет два уровня (0 и 1) в моей регрессионной модели.Я хотел бы взглянуть на то, как этот ковариат влияет на среднее значение моих кривых регрессии и значения дисперсии путем поднабора, когда модель становится равной «1».К сожалению, я получаю сообщение об ошибке: 

Error in `contrasts<-`(`*tmp*`, value = contr.funs[1 + isOF[nn]]) : 
  contrasts can be applied only to factors with 2 or more levels

Мои данные находятся в фрейме данных, озаглавленном "erad.surv", и выглядят так: 

    Surv(Erad_Time, Event)        Dist Public (weights)
1                       12 -0.51385627      1         1
2                      39+  0.65870948      1         1
3                      39+  0.47563451      0         1
4                      39+  0.73255621      1         1
5                      39+  0.32500333      1         1
6                      39+ -0.60323102      1         1
7                       31  0.60391220      1         1

Наконец, моя модель регрессиитот, что приведен ниже.Я также попробовал, казалось бы, любую другую комбинацию команды «subset», но безрезультатно. 

model.Public <- flexsurvreg(formula = Surv(Erad_Time, Event) ~ Dist + Public, data = erad.surv, subset = (Public == "1"), dist=ig)

Буду признателен за любые мысли или помощь!

...