diff --git a/scripts/10_correction_batch_effects_unfiltered.Rmd b/scripts/10_correction_batch_effects_unfiltered.Rmd
index 31d08a5..ce10593 100644
--- a/scripts/10_correction_batch_effects_unfiltered.Rmd
+++ b/scripts/10_correction_batch_effects_unfiltered.Rmd
@@ -78,7 +78,7 @@ if (correction_variable_1 == "PLEASE FILL IN"){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -237,7 +237,7 @@ if (correction_variable_2 == "PLEASE FILL IN"){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -396,7 +396,7 @@ if (correction_variable_3 == "PLEASE FILL IN"){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
diff --git a/scripts/8_removal_of_outliers_detect_batch_effects.Rmd b/scripts/8_removal_of_outliers_detect_batch_effects.Rmd
index 73f59b4..81f2c90 100644
--- a/scripts/8_removal_of_outliers_detect_batch_effects.Rmd
+++ b/scripts/8_removal_of_outliers_detect_batch_effects.Rmd
@@ -117,7 +117,7 @@ variance_explained <- data.frame(PC = paste0("PC",1:ncol(PCs)),
 
 variance_explained %>%
   dplyr::slice(1:10) %>%
-  ggplot(aes(x = PC, y = var_explained, group = 1)) +
+  ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
   geom_point() +
   geom_line() +
   geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -209,7 +209,7 @@ saveRDS(PCA_object, paste0(user_choices$project_name, "/reports/PCA_object_filte
 pdf(paste0(user_choices$project_name, "/reports/PCA_screeplot_filtered.pdf"))
 variance_explained %>%
   dplyr::slice(1:10) %>%
-  ggplot(aes(x = PC, y = var_explained, group = 1)) +
+  ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
   geom_point() +
   geom_line() +
   geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -252,7 +252,7 @@ if(!identical(extreme_outliers_filtered %>% pull(personid),character(0))){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -344,7 +344,7 @@ if(!identical(extreme_outliers_filtered %>% pull(personid),character(0))){
   pdf(paste0(user_choices$project_name, "/reports/PCA_screeplot_filtered_no_outliers.pdf"))
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -463,7 +463,7 @@ variance_explained <- data.frame(PC = paste0("PC",1:ncol(PCs)),
 
 variance_explained %>%
   dplyr::slice(1:10) %>%
-  ggplot(aes(x = PC, y = var_explained, group = 1)) +
+  ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
   geom_point() +
   geom_line() +
   geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -555,7 +555,7 @@ saveRDS(PCA_object, paste0(user_choices$project_name, "/reports/PCA_object_unfil
 pdf(paste0(user_choices$project_name, "/reports/PCA_screeplot_unfiltered.pdf"))
 variance_explained %>%
   dplyr::slice(1:10) %>%
-  ggplot(aes(x = PC, y = var_explained, group = 1)) +
+  ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
   geom_point() +
   geom_line() +
   geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -599,7 +599,7 @@ if(!identical(extreme_outliers_unfiltered %>% pull(personid),character(0))){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -691,7 +691,7 @@ if(!identical(extreme_outliers_unfiltered %>% pull(personid),character(0))){
   pdf(paste0(user_choices$project_name, "/reports/PCA_screeplot_unfiltered_no_outliers.pdf"))
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
diff --git a/scripts/9_correction_batch_effects_filtered.Rmd b/scripts/9_correction_batch_effects_filtered.Rmd
index 6091378..28c4525 100644
--- a/scripts/9_correction_batch_effects_filtered.Rmd
+++ b/scripts/9_correction_batch_effects_filtered.Rmd
@@ -88,7 +88,7 @@ if (correction_variable_1 == "PLEASE FILL IN"){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -247,7 +247,7 @@ if (correction_variable_2 == "PLEASE FILL IN"){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +
@@ -406,7 +406,7 @@ if (correction_variable_3 == "PLEASE FILL IN"){
   
   variance_explained %>%
     dplyr::slice(1:10) %>%
-    ggplot(aes(x = PC, y = var_explained, group = 1)) +
+    ggplot(aes(x = factor(PC, level = PC), y = var_explained, group = 1)) +
     geom_point() +
     geom_line() +
     geom_text(aes(label = cumm_var_explained), vjust = 1.5) +