J'analyse une population de patients que j'ai séparée en 2 groupes (obèse et non obèses). Je souhaiterais apparier des patients obèses avec des non obèses ayant certaines caractéristiques en commun (age, score de gravité...)
J'utilise la fonction "Match" et j'arrive à apparier les patients en choisissant une différence maximale commune pour toutes les variables mais je n'arrive pas à paramétrer la distance différente pour chacune des variables d'intérêt:
Ce script fonctionne
Code : Tout sélectionner
Appariement<-Match(Y=NULL, Tr=df$obese, X=cbind(df$Age,df$score),
distance.tolerance =c(5),
caliper=1,replace=F,M=1,exact=T)J'obtiens bien un appariement où les patients non obèses et les patients obèses ont une différence d'âge maximale de 5 ans et une différence de leur score de gravité maximale de 5 (32 patients appariés sur 35 potentiels)
Par contre si je veux une différence d'âge max de 5 ans et une différence de score de 3 ans, avec le script suivant:
Code : Tout sélectionner
Appariement<-Match(Y=NULL, Tr=df$obese, X=cbind(df$Age,df$score),
distance.tolerance =c(5,3),
caliper=1,replace=F,M=1,exact=T)
J'obtiens la réponse suivante:
Warning messages:
1: In if (distance.tolerance < 0) { :
the condition has length > 1 and only the first element will be used
2: In ecaliper[i] <- distance.tolerance :
number of items to replace is not a multiple of replacement length
3: In ecaliper[i] <- distance.tolerance :
number of items to replace is not a multiple of replacement length
Après m'être renseigné et avoir lu la notice du package, je crois que "distance.tolerance" ne peut contenir qu'un chiffre (précision des mesures) et que l'ondevrait paramétrer le caliper. J'ai essayé mais n'obtiens pas de résulat satisfaisant (appariement de seulement 3 patients même si je règle le caliper pour avoir la différence maximale de "5" comme pour la procédure ayant fonctionné plus haut).
Est-ce que quelqu'un pourrait m'indiquer un moyen simple d'apparier mes patients en paramétrant chacune des variables d'intérêt?
Merci beaucoup de votre aide.
Tài Pham