je souhaiterais intégrer dans une fonction une régression ordinale sur deux variables avec la fonction polr().
Voici le code :
Code : Tout sélectionner
RegLog<-function(Data.log,col_Y,col_variety,col_rep){
j2<-col_Y[1]-1
fit.polr<-vector("list",length = length(col_Y))# Initialisation liste modeles
summary.polr<-vector("list",length = length(col_Y))# Initialisation liste Resumes
for (j in col_Y){
fit.polr[[j-j2]]<-polr(Data.log[,j]~Data.log[,col_variety]+Data.log[,col_rep],data=Data.log)# modeles polr
summary.polr[[j-j2]]<-summary(fit.polr[[j-j2]])# resumes des modeles
}
}
col_Y est un vecteur regroupant les numéros de colonnes des variables d'intérêt. Lorsque je lance la fonction RegLog sur la fonction polr(), ça marche. Mais l'étape du summary non...
Fit.polr est une liste regroupant les deux modèles et summary.polr les résumés de chaque modèle.
L'erreur suivante est renvoyée :
Code : Tout sélectionner
Error in eval(expr, envir, enclos) : object 'Data.log' not found
Called from: eval(expr, envir, enclos)
Auriez-vous une explication et un moyen de remédier à ce problème?
Merci de votre aide,
Chloé