Stats anx plots

library(corrplot)
library(GGally)
library(readr)
library(car)
library(plotly)
library(ggplot2)
library(RColorBrewer)
library(ggpubr)
library(colortools)
library(cowplot)

setwd("~/Google Drive/PhD/Projects/Stats anxiety/data")

#Read data in

stats.anx.sona.csv <- read_csv('stats.anx.merged.cleaned.may11.v5.csv')

stats.anx.sona.csv$stats.anx <- rowMeans(subset(stats.anx.sona.csv,select = c(stat.anx.1:stat.anx.16)), na.rm=TRUE)

# Correlation plot
# variables for correlations
stats.anx<- stats.anx.sona.csv[,c("gender.category", "program.type", "faculty", "self.efficacy","asi",  "frost.com",    "frost.da", "stat.anx.tc",  "stat.anx.i",
"stat.anx.ah",  "stat.anx.ws", "stat.anx.fst",  "stat.anx.sc",  "math.anx", "stats.anx")]

#renaming for plot
colnames(stats.anx)[1] <- "Gender"
colnames(stats.anx)[2] <- "Program"
colnames(stats.anx)[3] <- "Faculty"
colnames(stats.anx)[4] <- "Self Efficacy"
colnames(stats.anx)[5] <- "Anxiety Sensitivity"
colnames(stats.anx)[6] <- "Concern Over Mistakes"
colnames(stats.anx)[7] <- "Doubts About Actions"
colnames(stats.anx)[8] <- "Statistics Test/Class Anxiety"
colnames(stats.anx)[9] <- "Statistics Interpretation Anxiety"
colnames(stats.anx)[10] <- "Asking for Statistics Help"
colnames(stats.anx)[11] <- "Worth of Statistics"
colnames(stats.anx)[12] <- "Fear of Statistics Teachers"
colnames(stats.anx)[13] <- "Computational Self-Concept"
colnames(stats.anx)[14] <- "Math Anxiety"
colnames(stats.anx)[15] <- "Statistics Anxiety"

#correlation Matrix
corrplot(cor(stats.anx, use="pairwise.complete.obs"), method="color",type="upper", 
                         addCoef.col = "black", tl.col="black", tl.srt=50, p.mat = cor1$p,sig.level = 0.05, 
                 insig = "blank", diag=FALSE, number.cex= 8/ncol(stats.anx))

#self-efficacy + statistics anxiety
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=stats.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and statistics anxiety") + 
  xlab("Self Efficacy") + 
  ylab("Statistics Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#self-efficacy + test/class anxiety
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=stat.anx.tc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and statistics test/class anxiety") + 
  xlab("Self Efficacy") + 
  ylab("Statistics Test/Class Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#self-efficacy + interpretation anxiety
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=stat.anx.i)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and statistics interpretation anxiety") + 
  xlab("Self Efficacy") + 
  ylab("Statistics Interpretation Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#self-efficacy + fear of asking for help
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=stat.anx.ah)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and asking for statistics help anxiety") + 
  xlab("Self Efficacy") + 
  ylab("Asking for Statistics Help Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#self-efficacy + worth of statistics
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=stat.anx.ws)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and percieved worth of statistics") + 
  xlab("Self Efficacy") + 
  ylab("Percieved Worth of Statistics") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#self-efficacy + fear of statistics teacher
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=stat.anx.fst)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and fear of statistics teachers") + 
  xlab("Self Efficacy") + 
  ylab("Fear of Statistics Teachers") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#self-efficacy + computational self-concept
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=stat.anx.sc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and computational self-concept") + 
  xlab("Self Efficacy") + 
  ylab("Computational self-concept") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#self-efficacy + math anxiety
ggplot(stats.anx.sona.csv, aes(x=self.efficacy, y=math.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between self efficacy and math anxiety") + 
  xlab("Self Efficacy") + 
  ylab("Math Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1))

#anxiety sensitivity + statistics anxiety
ggplot(stats.anx.sona.csv, aes(x=asi, y=stats.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and statistics anxiety") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Statistics Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#anxiety sensitivity + test/class anxiety
ggplot(stats.anx.sona.csv, aes(x=asi, y=stat.anx.tc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and statistics test/class anxiety") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Statistics Test/Class Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#anxiety sensitivity + interpretation anxiety
ggplot(stats.anx.sona.csv, aes(x=asi, y=stat.anx.i)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and statistics interpretation anxiety") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Statistics Interpretation Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#anxiety sensitivity + fear of asking for help
ggplot(stats.anx.sona.csv, aes(x=asi, y=stat.anx.ah)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and asking for statistics help anxiety") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Asking for Statistics Help Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#anxiety sensitivity + worth of statistics
ggplot(stats.anx.sona.csv, aes(x=asi, y=stat.anx.ws)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and percieved worth of statistics") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Percieved Worth of Statistics") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#anxiety sensitivity + fear of statistics teacher
ggplot(stats.anx.sona.csv, aes(x=asi, y=stat.anx.fst)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and fear of statistics teachers") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Fear of Statistics Teachers") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#anxiety sensitivity + computational self-concept
ggplot(stats.anx.sona.csv, aes(x=asi, y=stat.anx.sc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and computational self-concept") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Computational self-concept") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#anxiety sensitivity + math anxiety
ggplot(stats.anx.sona.csv, aes(x=asi, y=math.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between anxiety sensitivity and math anxiety") + 
  xlab("Anxiety Sensitivity") + 
  ylab("Math Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1))

#concern over mistakes + statistics anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=stats.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and statistics anxiety") + 
  xlab("Concern over Mistakes") + 
  ylab("Statistics Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#concern over mistakes + test/class anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=stat.anx.tc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and statistics test/class anxiety") + 
  xlab("Concern over Mistakes") + 
  ylab("Statistics Test/Class Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#concern over mistakes + interpretation anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=stat.anx.i)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and statistics interpretation anxiety") + 
  xlab("Concern over Mistakes") + 
  ylab("Statistics Interpretation Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#concern over mistakes + fear of asking for help
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=stat.anx.ah)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and asking for statistics help anxiety") + 
  xlab("Concern over Mistakes") + 
  ylab("Asking for Statistics Help Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#concern over mistakes + worth of statistics
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=stat.anx.ws)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and percieved worth of statistics") + 
  xlab("Concern over Mistakes") + 
  ylab("Percieved Worth of Statistics") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#concern over mistakes + fear of statistics teacher
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=stat.anx.fst)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and fear of statistics teachers") + 
  xlab("Concern over Mistakes") + 
  ylab("Fear of Statistics Teachers") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#concern over mistakes + computational self-concept
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=stat.anx.sc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and computational self-concept") + 
  xlab("Concern over Mistakes") + 
  ylab("Computational self-concept") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#concern over mistakes + math anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.com, y=math.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between concern over mistakes and math anxiety") + 
  xlab("Concern over Mistakes") + 
  ylab("Math Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1))

#doubts about actions + statistics anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=stats.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and statistics anxiety") + 
  xlab("Doubts About Actions") + 
  ylab("Statistics Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#doubts about actions + test/class anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=stat.anx.tc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and statistics test/class anxiety") + 
  xlab("Doubts About Actions") + 
  ylab("Statistics Test/Class Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#doubts about actions + interpretation anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=stat.anx.i)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and statistics interpretation anxiety") + 
  xlab("Doubts About Actions") + 
  ylab("Statistics Interpretation Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#doubts about actions + fear of asking for help
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=stat.anx.ah)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and asking for statistics help anxiety") + 
  xlab("Doubts About Actions") + 
  ylab("Asking for Statistics Help Anxiety") +
  theme_cowplot() +
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#doubts about actions + worth of statistics
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=stat.anx.ws)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and percieved worth of statistics") + 
  xlab("Doubts About Actions") + 
  ylab("Percieved Worth of Statistics") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#doubts about actions + fear of statistics teacher
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=stat.anx.fst)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and fear of statistics teachers") + 
  xlab("Doubts About Actions") + 
  ylab("Fear of Statistics Teachers") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#doubts about actions + computational self-concept
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=stat.anx.sc)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and computational self-concept") + 
  xlab("Doubts About Actions") + 
  ylab("Computational self-concept") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1)) 

#doubts about actions + math anxiety
ggplot(stats.anx.sona.csv, aes(x=frost.da, y=math.anx)) + 
  geom_point(size = 3, alpha = 0.3) +
  geom_smooth(method=lm, se=TRUE, fullrange=FALSE, level=0.95, color ="hotpink3", fill = "pink") +
  ggtitle("Association between doubts about actions and math anxiety") + 
  xlab("Doubts About Actions") + 
  ylab("Math Anxiety") +
  theme_cowplot() + 
    theme(plot.title = element_text(size = 13, face = "bold")) +
  stat_cor(p.accuracy = 0.0001, r.accuracy = 0.01) + 
  scale_x_continuous(limits = c(1,5), breaks = seq(0,5,1)) +
  scale_y_continuous(limits = c(1,6), breaks = seq(0,5,1))