suppressPackageStartupMessages({
library(dplyr)
library(ggformula)
library(mosaic)
library(supernova)
library(lsr)})
EAMMIdata <- read.csv(file = "https://bit.ly/2NGiyrS", header = TRUE)
head(EAMMIdata)
gf_histogram(~ Stress, data = EAMMIdata)
Warning message:
“Removed 13 rows containing non-finite values (stat_bin).”
gf_boxplot(~Stress, data = EAMMIdata)
Warning message:
“Removed 13 rows containing non-finite values (stat_boxplot).”
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))
Warning message:
“geom_vline(): Ignoring `mapping` because `xintercept` was provided.”
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)
Warning message:
“Removed 1039 rows containing missing values (geom_point).”
gf_point(Stress~age, data = EAMMIdata)
Warning message:
“Removed 1039 rows containing missing values (geom_point).”
Condition.model <- lm(Stress~age, data = EAMMIdata)
Condition.model
supernova(Condition.model)
Note: 1050 cases removed due to missing value(s).
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))
Warning message:
“geom_vline(): Ignoring `mapping` because `xintercept` was provided.”
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)
Note: 1050 cases removed due to missing value(s).
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")
Warning message:
“Removed 13 rows containing non-finite values (stat_boxplot).”
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~.)
Warning message:
“Removed 13 rows containing non-finite values (stat_bin).”
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)
Note: 1050 cases removed due to missing value(s).
SDob1<- do(1000)* b1(Stress~age, data= resample(EAMMIdata, 3134))
gf_histogram(~b1, data = SDob1)%>%
gf_vline(xintercept = b1 (Stress~age, data = EAMMIdata))
Warning message:
“geom_vline(): Ignoring `mapping` because `xintercept` was provided.”
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)
Note: 1050 cases removed due to missing value(s).