coxph () X matriz considerada singular;
Estoy teniendo algunos problemas para usar coxph (). Tengo dos variables categóricas: "tecnologia" y "pais", y quiero evaluar el posible efecto de interacción de "pais" en "tecnologia". "Tecnologia" es un factor variable con 2 niveles: gps y convencional. Y "pais" como 2 niveles: PT y ES. No tengo idea de por qué esta advertencia sigue apareciendo. Aquí está el código y la salida:
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
Estoy abriendo otra pregunta sobre este tema, aunque hice una similar hace algunos meses, porque estoy enfrentando el mismo problema otra vez, con otros datos. Y esta vez estoy seguro de que no es un problema relacionado con los datos.
¿Puede alguien ayudarme? Gracias
ACTUALIZAR: El problema no parece ser una clasificación perfecta.
> 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