У меня проблемы с использованием coxph(). Я имею две категориальные переменные: "tecnologia" и "pais", и я хочу оценить возможный эффект взаимодействия "pais" на "tecnologia". "Tecnologia" - это переменный фактор с двумя уровнями: gps и convencional. И "pais" как 2 уровня: PT и ES. Я не знаю, почему это предупреждение продолжает появляться. Здесь код и вывод:
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
Я открываю еще один вопрос по этому вопросу, хотя несколько месяцев назад я сделал подобное, потому что снова столкнулся с той же проблемой, с другими данными. И на этот раз я уверен, что это не проблема, связанная с данными.
Может кто-нибудь мне помочь? спасибо
ОБНОВЛЕНИЕ: проблема не кажется идеальной классификацией
> 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