Home Forum Statistica con R risposta categoriale 4 categorie e confronto ANOVA

This topic contains 3 replies and has 2 voices.

Viewing 4 posts - 1 through 4 (of 4 total)
  • Author
    Posts
  • #5176

    silvia9
    Participant

    Ciao! devo fare un modello per una y categoriale con 4 modalità (e una sola esplicativa binomiale)
    poi effettuare un’anova per il confronto.
    Non riesco però a trovare il comando in R per far il modello, qual’è il più indicato?
    Grazie per l’aiuto!

    #5325

    pdeninis
    Participant

    Ciao a te.

    Non so se ho capito. Io faccio semplicemente così:

    model<-lm(vardip ~ factor, dataframe) #Model è il modello di regressione, i parametri 
                                         #passati la variabile dipendente ed il fattore
    summary(model)                       #Per avere i parametri della regressione
    summary(aov(model))                  #Per avere l'anova
    
    library(multcomp)
    library(sandwich)
    posthocAnova = glht(model, linfct = mcp(TREATC = "Tukey"), vcov = sandwich)
    summary(posthocAnova)                #Per i test post-hoc tra i livelli del fattore
    
    #Group mean
    aovAnova<-aov(vardip ~ factor)
    library(effects)
    EMMAnova <- effect("factor", aovAnova)
    summary(EMMAnova)

    Il risultato dell’ANOVA può essere diverso da quello dei software che usano l’SS Type III. R usa di default l’I.

    Se vuoi avere l’SS Type III

    #Ancova SS Type III
    modelAnova<-lm(vardip ~ factor, contrast=list(factor=contr.sum))
    library(car)
    Anova(modelAnova, Type="III")
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    #5337

    pdeninis
    Participant

    Il risultato dell’ANOVA può essere diverso da quello dei software che usano l’SS Type III. R usa di default l’I.

    La mia frase era per generalizzare a più tipi di modelli diversi.
    Di questa differenza puoi accorgertene solo se includi altri regressori non ortogonali, però. Con il solo fattore a 4 livelli non ci saranno differenze.

    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    #5339

    pdeninis
    Participant

    Avevo il dubbio di aver capito male… Mi scuso: non vuoi fare una semplice anova.

    Allora: y categoriale a 4 modalità e x binomiale?
    Leggi questo e spiega se pensi che faccia la caso tuo.

    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
    • This reply was modified 7 months, 3 weeks ago by  pdeninis.
Viewing 4 posts - 1 through 4 (of 4 total)

You must be logged in to reply to this topic.