IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

R Discussion :

comment traiter des données binaires avec glmer


Sujet :

R

Mode arborescent

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre confirmé
    Femme Profil pro
    Étudiant
    Inscrit en
    Avril 2017
    Messages
    80
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France, Aisne (Picardie)

    Informations professionnelles :
    Activité : Étudiant

    Informations forums :
    Inscription : Avril 2017
    Messages : 80
    Par défaut comment traiter des données binaires avec glmer
    Bonjour,

    j'ai des données non normales alors j'ai transformé ces données en données binaires (0 = pas de sucre et 1 = du sucre), j'aimerais les analyser, sauf que j'ai des messages d'erreur que j'utilise glm (sans tenir compte des répétitions) ou bien glmer avec ou sans interaction et en tenant compte des répétitions.

    ma variable mesurée est : sug (binary data: 0 pour pas de sucre mesuré et 1 si du sucre à été mesuré)
    mes effets fixes: temp (2 Temperatures différentes) et var (6 variétés différentes testées) + l'interaction des deux
    mes réplicats = rep (le sucre à été mesuré sur 5 réplicats biologiques)

    J'ai donc un split plot design et j'aimerais voir les différences variétales, entre températures et pour les interactions au niveau du sucres.

    Pouvez-vous m'aider svp? Voici mon script ci-dessous et le fichier "fille" en pièce-jointe.
    Normalement je dois obtenir un effet de l'interaction temp * var car je sais que j'ai 3 variétés qui ne sucrent pas pour temp = 8 mais qui sucrent pour temp = 4 et pour les 3 autres variétés elles ne sucrent pas pour les 2 températures.

    Merci par avance.
    Margot

    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
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
    29
    30
    31
    32
    33
    34
    35
    36
    37
    38
    39
    file<-read.csv2("file.csv", na.strings = "NA",header = TRUE)
    head(file)
    str(file)
    file$sug<-as.numeric(file$sug)
    file$temp<-as.factor(file$temp)
    library(lme4)
    library(car)
    ##Measured variable : sug (binary data: 0 means no sugar and 1 means sugar)
    ##fixed effect to observed: temp, var and interaction temp, var
    ##replicates = rep (5 biological replicates per sug analysis)
     
    #split plot design?
       # 2 temperatures (temp)
       # 6 varieties (var)
       #observation of the sugar content (sug) on 5 biological replicates
     
    #model without random factor (rep) / glm function
    model_global = glm(sug ~ temp * var , family="binomial", data =file)
    Anova(model_global)
        #warning:
            # Warning messages:
                   #   1: glm.fit: fitted probabilities numerically 0 or 1 occurred
                    # 2: glm.fit: fitted probabilities numerically 0 or 1 occurred
                   # 3: glm.fit: fitted probabilities numerically 0 or 1 occurred
     
    #model with random factor (rep) / glmer fuction
    model_global_int = glmer(sug ~ temp * var + (1 | rep), family=binomial, data =file)
    Anova(model_global_int)
     
         #Warning messages:
           #   1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
              #                     unable to evaluate scaled gradient
                #                   2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv,  :
                                       # Hessian is numerically singular: parameters are not uniquely determined
    model_global_plus = glmer(sug ~ temp + var + (1 | rep), family=binomial, data =file)
    Anova(model_global_plus)
     
          #Error in pwrssUpdate(pp, resp, tol = tolPwrss, GQmat = GQmat, compDev = compDev,  :
                         #  (maxstephalfit) PIRLS step-halvings failed to reduce deviance in pwrssUpdate
    Fichiers attachés Fichiers attachés

Discussions similaires

  1. [XQuery] traiter des données XML avec xquery
    Par fatjoe dans le forum Format d'échange (XML, JSON...)
    Réponses: 6
    Dernier message: 29/03/2010, 02h17
  2. Réponses: 2
    Dernier message: 09/01/2009, 17h33
  3. Réponses: 9
    Dernier message: 04/06/2008, 11h38
  4. Comment récupérer des données enregistrées avec Rapid File ?
    Par didmarj2a dans le forum Bases de données
    Réponses: 2
    Dernier message: 21/03/2008, 09h49
  5. Comment traiter les données reçues avec WinSock ?
    Par Shredder dans le forum Web & réseau
    Réponses: 3
    Dernier message: 01/10/2007, 13h54

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo