Coxph () Macierz X uważana za liczbę pojedynczą;
Mam problemy z używaniem koksów (). Mam dwie zmienne kategoryczne: „tecnologia” i „pais” i chcę ocenić możliwy wpływ interakcji „pais” na „tecnologia”. „Tecnologia” jest zmiennym czynnikiem z 2 poziomami: GPS i konwencjonalnym. I „pais” jako 2 poziomy: PT i ES. Nie mam pojęcia, dlaczego to ostrzeżenie ciągle pojawia się. Oto kod i wyjście:
cox_AC<-coxph(Surv(dados_temp$dias_seg,dados_temp$status)~tecnologia*pais,data=dados_temp)
Warning message:
In coxph(Surv(dados_temp$dias_seg, dados_temp$status) ~ tecnologia * :
X matrix deemed to be singular; variable 3
> cox_AC
Call:
coxph(formula = Surv(dados_temp$dias_seg, dados_temp$status) ~
tecnologia * pais, data = dados_temp)
coef exp(coef) se(coef) z p
tecnologiagps -0.152 0.859 0.400 -0.38 7e-01
paisPT 1.469 4.345 0.406 3.62 3e-04
tecnologiagps:paisPT NA NA 0.000 NA NA
Likelihood ratio test=23.8 on 2 df, p=6.82e-06 n= 127, number of events= 64
Otwieram kolejne pytanie na ten temat, chociaż zrobiłem podobne pytanie kilka miesięcy temu, ponieważ znów mam ten sam problem z innymi danymi. Tym razem jestem pewien, że nie jest to problem związany z danymi.
Czy ktoś może mi pomóc? Dziękuję Ci
AKTUALIZACJA: Problem nie wydaje się być doskonałą klasyfikacją
> xtabs(~status+tecnologia,data=dados)
tecnologia
status conv doppler gps
0 39 6 24
1 30 3 34
> xtabs(~status+pais,data=dados)
pais
status ES PT
0 71 8
1 49 28
> xtabs(~tecnologia+pais,data=dados)
pais
tecnologia ES PT
conv 69 0
doppler 1 8
gps 30 28