diff --git a/.Rproj.user/DFCA1C7C/sources/prop/421CC383 b/.Rproj.user/DFCA1C7C/sources/prop/421CC383
index 7ddb6a8..0f6b30e 100644
--- a/.Rproj.user/DFCA1C7C/sources/prop/421CC383
+++ b/.Rproj.user/DFCA1C7C/sources/prop/421CC383
@@ -1,5 +1,5 @@
 {
     "tempName": "Untitled1",
     "cursorPosition": "128,0",
-    "scrollLine": "119"
+    "scrollLine": "140"
 }
\ No newline at end of file
diff --git a/demographics.R b/demographics.R
index f849a32..2f6a3ea 100644
--- a/demographics.R
+++ b/demographics.R
@@ -1,5 +1,11 @@
+library("writexl")
+
 demographics <- read.csv("oxford_hgshsa_260820_demographics.csv", sep = ";")
 
+demographics$missing <- data.export$has.na
+
+write_xlsx(demographics,"hab_taxon_demographics_2020_final_exported.xlsx")
+
 data <- read.csv("oxford_hgshsa_260820_modded.csv", sep = ";")
 
 data <- na.omit(data)
diff --git a/hability_taxon_data_prep.R b/hability_taxon_data_prep.R
index d8102f9..498d8c4 100644
--- a/hability_taxon_data_prep.R
+++ b/hability_taxon_data_prep.R
@@ -1,18 +1,34 @@
 library("RTaxometrics")
 library("dplyr")
+library("writexl")
 
+#Read the data
 data <- read.csv("oxford_hgshsa_260820_modded.csv", sep = ";")
 
-data <- na.omit(data)
 
+#Prepare a table for exporting for OSF
+data.export <- data
+
+#Generate the xlsx file for OSF
+data.export$has.na <- 0
+data.export$has.na[rowSums(is.na(data.export)) > 0] <- 1
+write_xlsx(data.export,"hab_taxon_2020_final_exported.xlsx")
+
+#Get the filtered data
+data <- data.export[data.export$has.na == 0,]
+
+#Subset the dataset for involuntariness scale
 invol <- select(data, INV.1:INV.12)
 
+#Classify data with base rates .15 and .10
 invol15 <- ClassifyCases(invol, p = 0.15)
 
 invol10 <- ClassifyCases(invol, p = 0.10)
 
+#Remove problematic items
 invol_purified <- select(data, INV.3:INV.10)
 
+#Classify data with base rates .15 and .10
 invol_pur15 <- ClassifyCases(invol_purified, p = 0.15)
 
 invol_pur10 <- ClassifyCases(invol_purified, p = 0.10)
@@ -25,7 +41,7 @@ invol$challenge <- invol$INV.4 + invol$INV.5 + invol$INV.6 + invol$INV.8 + invol
 invol$cogper <- invol$INV.9 + invol$INV.11
 invol$amnesia <- invol$INV.12
 
-
+#Subset subscales columns and classify
 invol_subscales15 <- select(invol, ideomotor, challenge, cogper, amnesia)
 invol_subscales15 <- ClassifyCases(invol_subscales15, p = 0.15)
 
@@ -41,7 +57,7 @@ invol_purified$ideomotor <- invol_purified$INV.3 + invol_purified$INV.7
 invol_purified$challenge <- invol_purified$INV.4 + invol_purified$INV.5 + invol_purified$INV.6 + invol_purified$INV.8 + invol_purified$INV.10
 invol_purified$cogper <- invol_purified$INV.9
 
-
+#Subset subscales columns and classify
 invol_subscales_pur15 <- select(invol_purified, ideomotor, challenge, cogper)
 invol_subscales_pur15 <- ClassifyCases(invol_subscales_pur15, p = 0.15)
 
@@ -56,6 +72,7 @@ invol_subscales_pur10 <- ClassifyCases(invol_subscales_pur10, p = 0.10)
 #Create a H. suggestibility subset
 hab <- select(data, HGSHS.A1:HGSHS.A11, HGSHS.A12)
 
+#Classify data with base rates .15 and .10
 hab15 <- ClassifyCases(hab, p = 0.15)
 
 hab10 <- ClassifyCases(hab, p = 0.10)
@@ -66,6 +83,7 @@ hab10 <- ClassifyCases(hab, p = 0.10)
 #Create a H. suggestibility subset without "bad" items
 hab_purified <- select(data,  HGSHS.A3:HGSHS.A10)
 
+#Classify data with base rates .15 and .10
 hab_pur15 <- ClassifyCases(hab_purified, p = 0.15)
 
 hab_pur10 <- ClassifyCases(hab_purified, p = 0.10)
@@ -80,7 +98,7 @@ hab$challenge <- hab$HGSHS.A4 + hab$HGSHS.A5 + hab$HGSHS.A6 + hab$HGSHS.A8 + hab
 hab$cogper <- hab$HGSHS.A9 + hab$HGSHS.A11
 hab$amnesia <- hab$HGSHS.A12
 
-
+#Subset subscales columns and classify
 hab_subsc15 <- select(hab, ideomotor, challenge, cogper, amnesia)
 hab_subsc15 <- ClassifyCases(hab_subsc15, p = 0.15)
 
@@ -96,7 +114,7 @@ hab_purified$ideomotor <- hab_purified$HGSHS.A3 + hab_purified$HGSHS.A7
 hab_purified$challenge <- hab_purified$HGSHS.A4 + hab_purified$HGSHS.A5 + hab_purified$HGSHS.A6 + hab_purified$HGSHS.A8 + hab_purified$HGSHS.A10
 hab_purified$cogper <- hab_purified$HGSHS.A9
 
-
+#Subset subscales columns and classify
 hab_subsc_pur15 <- select(hab, ideomotor, challenge, cogper)
 hab_subsc_pur15 <- ClassifyCases(hab_subsc_pur15, p = 0.15)
 
@@ -105,11 +123,7 @@ hab_subsc_pur10 <- select(hab, ideomotor, challenge, cogper)
 hab_subsc_pur10 <- ClassifyCases(hab_subsc_pur10, p = 0.10)
 
 
-# Performing the analysis
-
-## Involuntariness data: separate items
-
-### P = 0.15
+# Perform the analysis
 
 cat("invol15")
 RunTaxometrics(invol15)
diff --git a/plots.R b/plots.R
index 17e70a9..5b8a5ce 100644
--- a/plots.R
+++ b/plots.R
@@ -1,6 +1,10 @@
 library("ggplot2")
 library(gridExtra)
 library(grid)
+library(ggpubr)
+library(extrafont)
+font_import()
+loadfonts(device="win")       #Register fonts for Windows bitmap output
 
 
 # Involuntariness subscales; p = .15 --------------------------------------
@@ -6057,15 +6061,16 @@ mambac.plot.cat <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   ylim(mambac.y.range) +
   ggtitle(label = "Categorical comparison data",
           subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   scale_y_continuous(breaks=seq(round(mambac.y.range[1], digits = 1),
                                 round(mambac.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -6081,12 +6086,13 @@ mambac.plot.dim <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(mambac.y.range) +
   ggtitle(label = "Dimensional comparison data", subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -6108,14 +6114,15 @@ maxeig.plot.cat <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   scale_y_continuous(breaks=seq(round(maxeig.y.range[1], digits = 1),
                                 round(maxeig.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -6131,11 +6138,12 @@ maxeig.plot.dim <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -6162,14 +6170,15 @@ lmode.plot.cat <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   scale_y_continuous(breaks=seq(round(lmode.y.range[1], digits = 1),
                                 round(lmode.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -6191,11 +6200,12 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -6208,11 +6218,10 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
 
 #Grid
 
+grid.title <- textGrob(expression(bold("Involuntariness subscales ("~bolditalic("a priori")~"base rate: .15 )")), gp = gpar(fontsize = 14, fontfamily = "Arial"))
 grid.arrange(mambac.plot.cat, mambac.plot.dim, maxeig.plot.cat, maxeig.plot.dim, lmode.plot.cat, lmode.plot.dim,
              nrow = 3,
-             top=textGrob("Involuntariness subscales; p = .15",
-                          gp = gpar(fontsize = 14, fontface = 'bold')))
-
+             top=grid.title)
 
 
 
@@ -12271,15 +12280,16 @@ mambac.plot.cat <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   ylim(mambac.y.range) +
   ggtitle(label = "Categorical comparison data",
           subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   scale_y_continuous(breaks=seq(round(mambac.y.range[1], digits = 1),
                                 round(mambac.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -12295,12 +12305,13 @@ mambac.plot.dim <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(mambac.y.range) +
   ggtitle(label = "Dimensional comparison data", subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -12322,14 +12333,15 @@ maxeig.plot.cat <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   scale_y_continuous(breaks=seq(round(maxeig.y.range[1], digits = 1),
                                 round(maxeig.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -12345,11 +12357,12 @@ maxeig.plot.dim <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -12376,14 +12389,15 @@ lmode.plot.cat <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   scale_y_continuous(breaks=seq(round(lmode.y.range[1], digits = 1),
                                 round(lmode.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -12405,11 +12419,12 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -12421,10 +12436,10 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
   )
 
 #Grid
+grid.title <- textGrob(expression(bold("Involuntariness subscales ("~bolditalic("a priori")~"base rate: .10 )")), gp = gpar(fontsize = 14, fontfamily = "Arial"))
 grid.arrange(mambac.plot.cat, mambac.plot.dim, maxeig.plot.cat, maxeig.plot.dim, lmode.plot.cat, lmode.plot.dim,
              nrow = 3,
-             top=textGrob("Involuntariness subscales; p = .10",
-                          gp = gpar(fontsize = 14, fontface = 'bold')))
+             top=grid.title)
 
 
 
@@ -18496,15 +18511,16 @@ mambac.plot.cat <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   ylim(mambac.y.range) +
   ggtitle(label = "Categorical comparison data",
           subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   scale_y_continuous(breaks=seq(round(mambac.y.range[1], digits = 1),
                                 round(mambac.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -18520,12 +18536,13 @@ mambac.plot.dim <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(mambac.y.range) +
   ggtitle(label = "Dimensional comparison data", subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -18547,14 +18564,15 @@ maxeig.plot.cat <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   scale_y_continuous(breaks=seq(round(maxeig.y.range[1], digits = 1),
                                 round(maxeig.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -18570,11 +18588,12 @@ maxeig.plot.dim <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -18601,14 +18620,15 @@ lmode.plot.cat <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   scale_y_continuous(breaks=seq(round(lmode.y.range[1], digits = 1),
                                 round(lmode.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -18630,11 +18650,12 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -18646,13 +18667,13 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
   )
 
 #Grid
+grid.title <- textGrob(expression(bold("HGSHS:A subscales ("~bolditalic("a priori")~"base rate: .15 )")), gp = gpar(fontsize = 14, fontfamily = "Arial"))
 grid.arrange(mambac.plot.cat, mambac.plot.dim, maxeig.plot.cat, maxeig.plot.dim, lmode.plot.cat, lmode.plot.dim,
              nrow = 3,
-             top=textGrob("HGSHS:A subscales; p = .15",
-                          gp = gpar(fontsize = 14, fontface = 'bold')))
+             top=grid.title)
 
 
-# Suggestibility subscales; p = .15 --------------------------------------
+# Suggestibility subscales; p = .10 --------------------------------------
 # MAMBAC ------------------------------------------------------------------
 
 mambac.empcurve.x <- c(1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44,45,46,47,48,49,50)
@@ -24721,15 +24742,16 @@ mambac.plot.cat <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   ylim(mambac.y.range) +
   ggtitle(label = "Categorical comparison data",
           subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   scale_y_continuous(breaks=seq(round(mambac.y.range[1], digits = 1),
                                 round(mambac.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -24745,12 +24767,13 @@ mambac.plot.dim <- ggplot(mamb.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(mambac.y.range) +
   ggtitle(label = "Dimensional comparison data", subtitle = "MAMBAC") +
-  xlab("Cut") + ylab("Mean Difference") +
+  xlab("Cut") + ylab("Mean difference") +
   theme(
-    plot.title = element_text(size=16, hjust = 0.5),
-    plot.subtitle = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=16, hjust = 0.5, family = "Arial"),
+    plot.subtitle = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -24772,14 +24795,15 @@ maxeig.plot.cat <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   scale_y_continuous(breaks=seq(round(maxeig.y.range[1], digits = 1),
                                 round(maxeig.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -24795,11 +24819,12 @@ maxeig.plot.dim <- ggplot(maxe.emp, aes(x = xvalues, y = yvalues)) +
   geom_line(size = 1.5) +
   ylim(maxeig.y.range) +
   ggtitle(label = "MAXEIG") +
-  xlab("Score") + ylab("Eigen Value") +
+  xlab("Score") + ylab("Eigenvalue") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -24826,14 +24851,15 @@ lmode.plot.cat <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   scale_y_continuous(breaks=seq(round(lmode.y.range[1], digits = 1),
                                 round(lmode.y.range[2], digits = 1),
                                 by = 0.1)) +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     panel.grid.major = element_blank(),
     panel.grid.minor = element_blank(),
     panel.background = element_blank(),
@@ -24855,11 +24881,12 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
                    yend = lmode.vline2.cat.y[2]), size = 1) +
   ylim(lmode.y.range) +
   ggtitle(label = "L-Mode") +
-  xlab("Factor Score") + ylab("Density") +
+  xlab("Factor score") + ylab("Density") +
   theme(
-    plot.title = element_text(size=12, hjust = 0.5),
-    axis.title = element_text(size = 12),
-    axis.text = element_text(size = 12),
+    plot.title = element_text(size=12, hjust = 0.5, family = "Arial"),
+    axis.title = element_text(size = 12, family = "Arial"),
+    axis.text = element_text(colour = "black", size = 12, family = "Arial"),
+    axis.ticks = element_line(colour = "black"),
     axis.title.y = element_blank(),
     axis.text.y = element_blank(),
     axis.ticks.y = element_blank(),
@@ -24871,7 +24898,7 @@ lmode.plot.dim <- ggplot(lmo.emp, aes(x = xvalues, y = yvalues)) +
   )
 
 #Grid
+grid.title <- textGrob(expression(bold("HGSHS:A subscales ("~bolditalic("a priori")~"base rate: .10 )")), gp = gpar(fontsize = 14, fontfamily = "Arial"))
 grid.arrange(mambac.plot.cat, mambac.plot.dim, maxeig.plot.cat, maxeig.plot.dim, lmode.plot.cat, lmode.plot.dim,
              nrow = 3,
-             top=textGrob("HGSHS:A subscales; p = .10",
-             gp = gpar(fontsize = 14, fontface = 'bold')))
+             top=grid.title)