Analyseren van gegevensvarianten in modellen met R-dummies

Een variantieanalyse (ANOVA) is een zeer gebruikelijke techniek die wordt gebruikt met R om de gemiddelden te vergelijken tussen verschillende groepen gegevens. Om dit te illustreren, bekijkt u de dataset InsectSpray: >> str (InsectSprays). frame ': 72 obs. van 2 variabelen: $ count: num 10 7 20 14 14 12 10 23 17 20 … $ spray: factor met 6 niveaus "A", "B", "C", "D", …: 1 1 1 1 1 1 1 1 1 1 …

Deze dataset bevat de resultaten van een landbouwexperiment. Zes insecticiden werden getest op elk 12 velden en de onderzoekers telden het aantal vervelende bugs dat op elk veld bleef. Nu moeten de boeren weten of de insecticiden enig verschil maken, en zo ja, welke ze het best gebruiken. U beantwoordt deze vraag door de functie aov () te gebruiken om een ​​ANOVA uit te voeren.

Hoe het model te bouwen

Voor dit eenvoudige voorbeeld is het bouwen van het model een fluitje van een cent. U wilt in feite de middelen voor de variabele telling modelleren als een functie van de variabele spray. Je vertaalt dat naar R als volgt: >> AOVModel <- aov (count ~ spray, data = InsectSprays)

Je geeft twee argumenten door aan de functie aov () in deze regel code:

De formule telt ~ spray, die luidt als "tellen als een functie van spray"

• De argumentgegevens, waar u het gegevensframe opgeeft waarin de variabelen in de formule kunnen worden gevonden

  Elke modelleringsfunctie retourneert een modelobject met veel informatie over het ingerichte model. Zet dit modelobject altijd in een variabele. Op deze manier hoeft u het model niet opnieuw in te stellen wanneer u extra berekeningen moet uitvoeren.

• Hoe het modelobject te bekijken

Zoals met elk object, kunt u een modelobject bekijken door de naam ervan in de console te typen. Als u dat doet voor het objectmodel dat u hebt gemaakt, ziet u de volgende uitvoer:

>> AOVModel Oproep: aov (formula = count ~ spray, data = InsectSprays) Termen: spray Residuen Som van vierkanten 2668. 833 1015. 167 graden. of Freedom 5 66 Residuele standaardfout: 3. 921902 Geschatte effecten kunnen ongebalanceerd zijn

Dit zegt u niet zoveel, behalve het commando (of het

gesprek

) dat u hebt gebruikt om het model te bouwen en enige basisinformatie over het resultaat van de aanpassing.

In de uitvoer leest u ook dat de geschatte effecten mogelijk uit balans zijn. Dit is geen waarschuwing - het is een bericht dat is ingebouwd door de auteur van de functie aov (). Deze kan in twee situaties verschijnen: Je hebt niet hetzelfde aantal cases in elke groep. U hebt geen orthogonale contrasten ingesteld.

In dit geval is dit de tweede reden.