J'aimerais faire une régression logistique conditionnelle sur un jeu de données appariées.
J'utilise donc la fonction clogit() du package survival, comme ceci (analyse univariée):
Code : Tout sélectionner
clogit(STATUT~ETUDES+strata(ID),data = descT)
Et ça me renvoie ça:
Code : Tout sélectionner
Error in coxph(formula = Surv(rep(1, 1498L), STATUT) ~ ETUDES + strata(ID), :
Cox model doesn't support "mright" survival data
L'appariement est en 2 pour 1 (sur l'âge te le sexe), sauf pour 2 "groupes" où il est en 1 pour 1. J'ai refait la même chose en supprimant ces 2 groupes, et ça marque la même chose, alors a priori le problème n'est pas là.
Ma variable explicative ETUDES est une variable factorielle avec 4 groupes (codés 1 à 4).
Quel est d'après vous le truc qui coince?
Voici une partie de ma bdd:
Code : Tout sélectionner
head(descT)
# A tibble: 6 x 8
ID STATUT AGE SEXE ETUDES TABAC TABAC_PA
<int> <fctr> <dbl> <fctr> <fctr> <fctr> <fctr>
1 1400002 cas 43.92608 1 4 0 0
2 1400002 témoin 45.41821 1 2 2 3
3 1400002 témoin 46.20945 1 4 0 0
4 1400003 cas 49.85626 1 2 1 3
5 1400003 témoin 51.59206 1 4 2 3
6 1400003 témoin 47.94524 1 4 1 3
# ... with 1 more variables: ALCOOL <fctr>
>