add underlying data

00-scripts/step5/figures/figure2.R

@@ -216,8 +216,6 @@ legMatCol=matrix(
 
 ##########################################################################
 ##write tables
-data.frame(gd=Mr$GD_mean,richness=Mr$richness)
-
 write.table(data.frame(gd=Mr$GD_mean,richness=Mr$richness)
 ,"14-tables/figure_2a_underlying_data.tsv",col.names=T,row.names=F,sep="\t")
 write.table(data.frame(gd=Fr$GD_mean,richness=Fr$richness)

00-scripts/step5/figures/figure3.R

@@ -197,6 +197,9 @@ p_hp_F=ggplot(data=hp_df_F, aes(x=nom_var, y=hier_part)) +
         panel.border=element_rect(fill=NA,size=1),
         panel.background = element_blank())+
   ggtitle("(d)")
+## write underlying data tables
+write.table(hp_df_F,"14-tables/figure_3d_underlying_data.tsv",col.names=T,row.names=F,sep="\t")
+
 
 
 ##########################################################################
@@ -352,6 +355,9 @@ p_hp_M=ggplot(data=hp_df_M, aes(x=nom_var, y=hier_part)) +
         panel.border=element_rect(fill=NA,size=1),
         panel.background = element_blank())+
   ggtitle("(b)")
+## write underlying data tables
+write.table(hp_df_M,"14-tables/figure_3b_underlying_data.tsv",col.names=T,row.names=F,sep="\t")
+
 
 
 ##########################################################################

14-tables/figure_3b_underlying_data.tsv

+"nom_var"	"hier_part"
+"Sea surface
+temperature"	75.2149245547393
+"Regions"	9.21859678770979
+"Autocor"	3.4099407731038
+"Number of species"	9.02859567447895
+"Number of sequences"	3.12794220996813

14-tables/figure_3d_underlying_data.tsv

+"nom_var"	"hier_part"
+"Number of species"	2.11284354824658
+"Number of sequences"	1.00988243793768
+"Mean regional
+temperature"	10.4589871942955
+"Regions"	36.1126989042008
+"Average slope"	23.6091149269876
+"Autocor"	26.6964729883319