diff --git a/MuPoMo_main_multipop/MuPoMo_main_multipop.R b/MuPoMo_main_multipop/MuPoMo_main_multipop.R
index 44ca96a..006c025 100644
--- a/MuPoMo_main_multipop/MuPoMo_main_multipop.R
+++ b/MuPoMo_main_multipop/MuPoMo_main_multipop.R
@@ -64,7 +64,7 @@ for (i in 1:(loop.31 - 1)) {
         sep = ""))))
     assign(nam7, temp3)
 }
-reference0.temp = rowMeans(merge31, na.rm = TRUE)
+reference0.temp = rowMeans(merge30, na.rm = TRUE)
 reference0      = ts(reference0.temp, start = Sweden$year[1], frequency = 1)
 
 ## plot the reference curve among all 31 smoothed curves
@@ -152,4 +152,4 @@ for (j in 1:iteration) {
         lines(eval(parse(text = paste("sm.kt.", names.31[i], ".female", sep = ""))), col = "grey")
     }
     lines(eval(parse(text = paste("reference", j - 1, sep = ""))), col = "blue", lwd = 3, lty = 5)
-}
\ No newline at end of file
+}