Discussion :

[J-r_R]

BOnjour,
J'ai une erreur lorsque j'exécute une boucle avec des régressions logistiques, alors que hors de la boucle, le script fonctionne parfaitement pour toutes les valeurs que devrait prendre le compteur 'x' dans la boucle (j'obtiens mes résultats de régression en appliquant successivement x<-1 ; x<-2 ; etc avant de lancer le script contenu dans la boucle).
Le problème étant que je suis amené à refaire cette manip plusieurs fois et qu'à la longue ça prend du temps de faire tourner une boucle à la main...

Mon code est le suivant :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
sink("MonFichier.txt")
  for (x in 1:length(L)) {
    A<-paste0("M_",L[x])
    B<-paste0("Ms_",L[x])
    eval(parse(text=paste0(A," <- glm(",L[x]," ~ ",paste(E, collapse="+"),",family='binomial',data=MaSource)")))
    eval(parse(text=paste0(B," <- glm(",L[x]," ~ ",paste(E, collapse="+"),"+",paste(S, collapse="+"),",family='binomial',data=MaSource)")))
 
    print(eval(parse(text=paste0("summary(",A,")"))))
    print("Odds Ratio")
    print(eval(parse(text=paste0("exp(coefficients(",A,"))"))))
    print(eval(parse(text=paste0("summary(",B,")"))))
    print("Odds Ratio")
    print(eval(parse(text=paste0("exp(coefficients(",B,"))"))))
  }
  sink()

Où L contient les noms des variables à expliquer successivement
E et S contiennent les variables explicatives

Le tout sous en utilisant RStudio version 1.4.1103

Si quelqu'un peut me donner une piste, je lui en saurai infiniment gré !

[J-r_R]

J'ai oublié de préciser : le message d'erreur est le suivant :
"Error in weights * y : argument non numérique pour un opérateur binaire"

[olivier.decourt]

Bonjour.
Aucun souci avec le code ci-dessous, inspiré du tien.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
data(iris)
 
E <- c("Sepal.Width")
S <- c("Sepal.Length")
 
L <- c("Petal.Width","Petal.Length")
 
for (x in seq_along(L)){
  A<-paste0("M_",L[x])
  B<-paste0("Ms_",L[x])
  eval(parse(text=paste0(A," <- lm(",L[x]," ~ ",paste(E, collapse="+"),",data=iris)")))
  eval(parse(text=paste0(B," <- lm(",L[x]," ~ ",paste(E, collapse="+"),"+",paste(S, collapse="+"),",data=iris)")))
 
  print(eval(parse(text=paste0("summary(",A,")"))))
  print(eval(parse(text=paste0("summary(",B,")"))))
}

Est-ce que le message d'erreur apparaît à toutes les itérations de la boucle ou seulement à certaines ?
Es-tu sûr que le texte fourni à eval+parse est bien le code que tu attends ? (A priori oui mais à vérifier)
As-tu bien vérifié les types de tes variables ? (en particulier celles que tu veux expliquer)
Y a-t-il des pondérations dans tes modèles ?

[J-r_R]

BOnjour,

Merci beaucoup Olivier pour ta suggestion.

Je m'aperçois que je n'ai pas assez insisté sur un point : le bug survient pour les régressions logistiques (glm... family='binomial')
J'ai fait le même script avec des régressions linéaires qui fonctionne parfaitement (pour des variables à expliquer différentes évidemment).

Pour répondre à tes questions : la boule ne s'exécute pas du tout et je n'ai que le message d'erreur comme retour (une seule fois).
Le code dans eval est bon : il fonctionne si je le sors de la boucle (et pour toutes les valeurs de x) !
pareil pour les types de variables
Et non : pas de pondération !!!

C'est là qu'est tout le mystère... (pour moi en tout cas :-s

[olivier.decourt]

Les variables à expliquer ont quoi comme min et max ?

[J-r_R]

Toutes les variables à expliquer prennent les valeurs 0/1

[J-r_R]

Hello!

Finalement j'ai pu contourner le problème en utilisant une autre fonction (puisque le script fonctionnait très bien avec les régressions linéaires, je suppose que le problème vient d'un bug de la fonction glm)
J'ai donc utilisé la fonction vglm du package VGAM. En outre elle me permet d'utiliser également des variables dépendantes ordinales

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
sink("MonFichier.txt")
  for (x in 1:length(L)) {
    A<-paste0("M_",L[x])
    B<-paste0("Ms_",L[x])
    eval(parse(text=paste0(A," <- vglm(",L[x]," ~ ",paste(E, collapse="+"),",family=acat(),data=MaSource)")))
    eval(parse(text=paste0(B," <- vglm(",L[x]," ~ ",paste(E, collapse="+"),"+",paste(S, collapse="+"),",family=acat(),data=MaSource)")))
 
    print(eval(parse(text=paste0("summary(",A,")"))))
    print("Odds Ratio")
    print(eval(parse(text=paste0("exp(coefficients(",A,"))"))))
    print(eval(parse(text=paste0("summary(",B,")"))))
    print("Odds Ratio")
    print(eval(parse(text=paste0("exp(coefficients(",B,"))"))))
  }
  sink()

Un grand merci à toi, Olivier pour le temps que tu as passé sur ma question ! Et si tu trouve une solution/explication ça m'intéresse toujours néanmoins...