suppressPackageStartupMessages({
library(dplyr)
library(ggformula)
library(mosaic)
library(supernova)
library(lsr)})
EAMMIdata <- read.csv(file = "https://bit.ly/2NGiyrS", header = TRUE)
head(EAMMIdata)
nrow(EAMMIdata)
EAMMIdata <- filter(EAMMIdata, marriage1_2 < 100)
EAMMIdata <- filter(EAMMIdata, income != 1)
nrow(EAMMIdata)
gf_histogram(~income, data = EAMMIdata)%>%
gf_vline(xintercept = 4.504, color = "blue" )
Warning message:
“geom_vline(): Ignoring `mapping` because `xintercept` was provided.”
gf_boxplot(income~1, data = EAMMIdata)
empty.model <- lm(income ~ NULL, data = EAMMIdata)
empty.model
anova(empty.model)
supernova(empty.model)
income.stats <- favstats(~income, data = EAMMIdata)
SDoM <- do(10000) * mean(rnorm(1, mean = income.stats$mean, sd = income.stats$sd))
gf_histogram(~mean, data = SDoM, fill = "blue")
favstats(~income, data = EAMMIdata)
tally(~mean >= 4.504, data = SDoM, format = "proportion")
confint(empty.model, level = 0.95)
$$Y_i = b_0 + b_1 X_i + e_i$$
ERROR: Error in parse(text = x, srcfile = src): <text>:1:1: unexpected '$'
1: $
^
EAMMIdata$Parenting2Group <- ntile(EAMMIdata$marriage1_2, 2)
EAMMIdata$Parenting2Group <- factor(EAMMIdata$Parenting2Group, levels = c(1,2), labels = c("low", "high"))
gf_histogram(~income, data = EAMMIdata) %>%
gf_facet_grid(Parenting2Group~.)
gf_boxplot(Parenting2Group~income, data = EAMMIdata)
parentinggroupmodel <- lm(income~Parenting2Group, data = EAMMIdata)
parentinggroupmodel
anova(parentinggroupmodel)
superanova(parentinggroupmodel)
income.stats <- favstats(income ~Parenting2Group, data = EAMMIdata)
simincome <-rnorm(10000, mean = income.stats$mean, sd = income.stats$sd)
simpop <- data.frame(simincome)
gf_histogram(~simincome, data = simpop, fill = "red")
summary(parentinggroupmodel)
confint(parentinggroupmodel, level = 0.95)
SDoPRE <- do(1000) * PRE(income ~ shuffle(Parenting2Group), data = EAMMIdata)
samplePRE <- PRE(income ~ Parenting2Group, data = EAMMIdata)
tally(~PRE>samplePRE, data = SDoPRE)
xpf(samplePRE, df1 = 1, df2 = 2183)
gf_jitter(income~marriage1_2, data = EAMMIdata,width = .1, height = 0, alpha = 0.5)%>%
gf_lm(color = "coral")
gf_histogram(~income, data = EAMMIdata) %>%
gf_facet_grid(marriage1_2~.)
parenting.model <- lm(income~marriage1_2, data = EAMMIdata)
parenting.model
anova(parenting.model)
supernova(parenting.model)
income.stats <- favstats(income ~marriage1_2, data = EAMMIdata)
simincome <-rnorm(10000, mean = income.stats$mean, sd = income.stats$sd)
simpop <- data.frame(simincome)
gf_histogram(~simincome, data = simpop, fill = "darkorange2") %>%
gf_labs(title = "Simulated Population Distribution")
Warning message in rnorm(10000, mean = income.stats$mean, sd = income.stats$sd):
“NAs produced”
Warning message:
“Removed 3638 rows containing non-finite values (stat_bin).”
summary(parenting.model)
confint(parenting.model, level = .95)
#alpha level = .05
SDoPRE <- do(1000) * PRE(income ~ shuffle(marriage1_2), data = EAMMIdata)
samplePRE <- PRE(income ~ marriage1_2, data = EAMMIdata)
tally(~PRE>samplePRE, data = SDoPRE)
xpf(samplePRE, df1 = 1, df2 = 3079)