--- title: "Practicum 2: Stokstaartjes" output: html_document: code_download: yes highlight: tango number_sections: yes theme: cosmo toc: yes toc_float: yes pdf_document: toc: yes word_document: toc: yes --- Stokstaartjes, een mangoestensoort die voorkomt in zuidelijk Afrika, leven in groepen. Elke groep kent zijn dominant vrouwtje, die grotendeels het monopolie heeft op de voortplanting binnen de groep: ze is de moeder van meer dan tachtig procent van de jongen die hun eerste maand overleven. Hoe ze dit monopolie bewerktstelligt, is echter minder goed begrepen. Een mogelijke verklaring ligt in het feit dat het dominante vrouwtje ongeveer een maand voor ze bevalt de andere vrouwtjes terroriseert, wat voor hen vaak uitmondt in een tijdelijke verbanning uit de groep. In een onderzoek tracht men na te gaan of deze verbanning gerelateerd is met factoren die een nadelige invloed hebben op de voortplanting. De studiepopulatie waarop het onderzoek betrekking heeft, bedraagt individuen uit 15 willekeurig gekozen groepen stokstaartjes. In een eerste luik van het experiment werd het gewichtsverschil bekeken van 15 stokstaartjes die verbannen werden uit de groep. In een tweede luik van het experiment werden op gelijkaardige tijdstippen als het vertrek en de terugkomst van de verbannen vrouwtjes het gewichtsverschil gemeten van 15 andere stokstaartjes uit dezelfde regio, maar die in hun eigen groep konden blijven. Deze laatste groep van individuen worden aangeduid in de dataset met `verbannen == 0` . #Gewichtsverschil met en zonder verbanning Test op het 5% significantieniveau de onderstelling dat de gemiddelde gewichtsverandering bij vrouwtjes die verbannen waren uit de groep **verschilt van** vrouwtjes die in de groep konden blijven. Noteer tevens duidelijk uw conclusie. # R libraries inlezen ```{r} library(ggplot2) library(dplyr) library(tidyr) library(readr) ``` # Vraag 1 Dataset inlezen De dataset `stokstaartjesDat.txt` bestaat uit $60$ observaties van de volgende variabelen: - `gewicht` het gewicht van het stokstaartje - `tijdstip` tijdstip van de meting. `"voor"` indien gewogen werd voor de verbanning, anders `"na"` - `verbannen` $1$ indien verbannen, $0$ indien niet. ```{r} #In het tekstbestand zijn er drie kolomnamen gegeven op de eerste lijn; gewicht, tijdstip en verbannen. Er zijn echter vier waarden vanaf de tweede lijn, waarbij de eerste waarde een rij nummer is. Dit zal dus voor problemen zorgen bij het inlezen. Om dit probleem te verhelpen gebruiken we skip = 1 om de eerste rij met de kolomnamen over te slaan. Daarna geven we zelf kolomnamen door met col_names. Met col_types vermelden we specifiek welk type variabele iedere kolom is. Met col_select selecteren we dan de kolommen die we wel degelijk nodig hebben. stokstaartjes <- read_delim("https://raw.githubusercontent.com/statOmics/statistiekBasisCursusData/master/practicum3/stokstaartjesDat.txt", delim = "\ " , skip = 1, col_names = c("rij_nr","gewicht","tijdstip","verbannen"), col_types = list(col_integer(), col_double(), col_factor(), col_factor()), col_select = c("gewicht","tijdstip","verbannen")) stokstaartjes <- stokstaartjes %>% mutate( index = factor(c(rep(1:15,2),rep(16:30,2))) ) ```