diff --git a/02_Mesobromion_ExamineOutput.R b/02_Mesobromion_ExamineOutput.R
index b99cbff39df9ce244764c8110bd19a4b5355841c..52e9c37d7ca5d7c29a044983c72344a93793b288 100644
--- a/02_Mesobromion_ExamineOutput.R
+++ b/02_Mesobromion_ExamineOutput.R
@@ -399,7 +399,7 @@ corXY.ci.pa <- corXY.ci.pa %>%
   dplyr::select(Trait.comb, Test, n, ntraits, everything()) %>% 
   mutate(run="seq")
 
-rm( dataFiles1, corXY.all) 
+rm( dataFiles, corXY.all) 
 
 
 #### 2.2.2 One by one - Graph of r(XY)  ####
@@ -493,11 +493,11 @@ for(nn in 1:maxtraits){
   }
 }
 
-best.4traits <- corXY.ci %>% 
+best.4traits.pa <- corXY.ci.pa %>% 
   filter(as.character(Trait.comb)==best.progr[length(best.progr)]) %>% 
-  dplyr::select(trait1:trait8) %>% 
+  dplyr::select(trait1:trait4) %>% 
   mutate_all(~as.character(.))
-best.4traits <- as.character(best.4traits[1,1:4])
+best.4traits.pa <- as.character(best.4traits.pa[1,1:4])
 
 ### Create dataset with best combinations + all the one-way combinations
 mydata.best <- mydata %>%