EAMMIdata <- read.csv(file = "https://bit.ly/2NGiyrS", header = TRUE) head(EAMMIdata)

gf_histogram(~ Stress, data = EAMMIdata)

gf_boxplot(~Stress, data = EAMMIdata)

Stress.model <- lm(Stress ~ NULL, data = EAMMIdata) Stress.model

anova(Stress.model)

SDob0.model <- do(1000)*b0(Stress~NULL, data = resample(EAMMIdata, 3134)) gf_histogram(~b0, data = SDob0.model)%>% gf_vline(xintercept = b0(Stress ~ NULL, data = EAMMIdata))

SDob0.model <- do(1000)*b0(Stress ~ NULL, data = resample(EAMMIdata, 3134)) favstats(~b0, data = SDob0.model)

SDob1.model <- do(1000)*b0(Stress~NULL, data = resample(EAMMIdata,3134)) confint(SDob1.model)

gf_jitter(Stress~age, data= EAMMIdata)

gf_point(Stress~age, data = EAMMIdata)

Condition.model <- lm(Stress~age, data = EAMMIdata) Condition.model

supernova(Condition.model)

SDob1.model <- do(1000)*b1(Stress~age, data = resample(EAMMIdata,3134)) gf_histogram(~b1, data = SDob1.model)%>% gf_vline(xintercept = b1 (Stress~age, data = EAMMIdata))

SDob1.model <- do(1000)*b1(Stress~age, data = resample(EAMMIdata,3134)) favstats(~b1, data=SDob1.model)

SDob1.model <- do (1000)* b1(Stress~age, data = resample(EAMMIdata,3134)) confint(SDob1.model)

Condition.model <- lm(Stress~age, data = EAMMIdata) supernova(Condition.model) empty2.model <- do(1000)*b1(Stress~resample(age), data = EAMMIdata) sampleb1 <- b1(Condition.model) gf_histogram(~b1, data = empty2.model, fill = ~b1>sampleb1)

EAMMIdata$age <- ntile(EAMMIdata$age, 2) EAMMIdata$age <- factor(EAMMIdata$age, levels = c(1,2), labels = c("young", "old")) gf_boxplot(Stress~age, data=EAMMIdata, fill = "blue")

EAMMIdata$age <- ntile(EAMMIdata$age, 2) EAMMIdata$age <- factor(EAMMIdata$age, levels = c(1,2), labels = c("young", "old")) gf_dhistogram(~Stress, data=EAMMIdata) %>% gf_facet_grid(age~.)

EAMMIdata$age <- ntile(EAMMIdata$age, 2) EAMMIdata$age <- factor(EAMMIdata$age, levels = c(1,2), labels = c("young", "old")) Stress.mod <-lm(Stress~age, data=EAMMIdata) Stress.mod

supernova(Stress.mod)

SDob1<- do(1000)* b1(Stress~age, data= resample(EAMMIdata, 3134)) gf_histogram(~b1, data = SDob1)%>% gf_vline(xintercept = b1 (Stress~age, data = EAMMIdata))

SDob1.model<- do(1000)*b1(Stress~age, data=resample(EAMMIdata,3134)) favstats (~b1, data = SDob1.model)

SDob1.model<- do(1000)*b1(Stress~age, data=resample(EAMMIdata,3134)) confint(SDob1.model)

StressLevel.model<-lm(Stress~age, data=EAMMIdata) supernova(StressLevel.model) StressLev2.model <- do(1000) *b1(resample(Stress) ~age, data=EAMMIdata) modelb1<-b1(StressLevel.model) gf_histogram(~b1, data=StressLev2.model, fill= ~b1>modelb1)