@@ -818,8 +818,8 @@ ASV_physeq <- phyloseq(count_tab_phy, tax_tab_phy, sample_info_tab_phy)
818818richness_and_diversity_estimates_by_sample <- plot_richness(ASV_physeq , color = " groups" , measures = c(" Chao1" , " Shannon" ))
819819richness_and_diversity_estimates_by_group <- plot_richness(ASV_physeq , x = " groups" , color = " groups" , measures = c(" Chao1" , " Shannon" ))
820820
821- ggsave(paste0( " richness_and_diversity_estimates_by_sample_GLAmpSeq " , " .png" ) , plot = richness_and_diversity_estimates_by_sample )
822- ggsave(paste0( " richness_and_diversity_estimates_by_group_GLAmpSeq " , " .png" ) , plot = richness_and_diversity_estimates_by_group )
821+ ggsave(filename = " richness_and_diversity_estimates_by_sample_GLAmpSeq .png" , plot = richness_and_diversity_estimates_by_sample )
822+ ggsave(filename = " richness_and_diversity_estimates_by_group_GLAmpSeq .png" , plot = richness_and_diversity_estimates_by_group )
823823```
824824
825825** Parameter Definitions:**
@@ -861,10 +861,10 @@ relative_classes <- plot_bar(proportions_physeq, x = "groups", fill = "class")
861861samplewise_phyla <- plot_bar(proportions_physeq , fill = " phylum" )
862862samplewise_classes <- plot_bar(proportions_physeq , fill = " class" )
863863
864- ggsave(filename = " relative_phyla_GLAmpSeq" , " .png" , plot = relative_phyla )
865- ggsave(filename = " relative_classes_GLAmpSeq" , " .png" , plot = relative_classes )
866- ggsave(filename = " samplewise_relative_phyla_GLAmpSeq" , " .png" , plot = samplewise_phyla )
867- ggsave(filename = " samplewise_relative_classes_GLAmpSeq" , " .png" , plot = samplewise_classes )
864+ ggsave(filename = " relative_phyla_GLAmpSeq.png" , plot = relative_phyla )
865+ ggsave(filename = " relative_classes_GLAmpSeq.png" , plot = relative_classes )
866+ ggsave(filename = " samplewise_relative_phyla_GLAmpSeq.png" , plot = samplewise_phyla )
867+ ggsave(filename = " samplewise_relative_classes_GLAmpSeq.png" , plot = samplewise_classes )
868868```
869869
870870** Input Data:**
@@ -967,8 +967,8 @@ ordination_plot_u <- plot_ordination(vst_physeq, vst_pcoa, color = "groups") +
967967 annotate(" text" , x = Inf , y = - Inf , label = paste(" R2:" , toString(round(r2_value , 3 ))), hjust = 1.1 , vjust = - 2 , size = 4 )+
968968 annotate(" text" , x = Inf , y = - Inf , label = paste(" Pr(>F)" , toString(round(prf_value ,4 ))), hjust = 1.1 , vjust = - 0.5 , size = 4 )+ ggtitle(" PCoA" )
969969
970- ggsave(filename = paste0( beta_diversity_out_dir , output_prefix , " PCoA_w_labels_GLAmpSeq" , " .png" ) , plot = ordination_plot )
971- ggsave(filename = paste0( beta_diversity_out_dir , output_prefix , " PCoA_without_labels_GLAmpSeq" , " .png" ) , plot = ordination_plot_u )
970+ ggsave(filename = " PCoA_w_labels_GLAmpSeq.png" , plot = ordination_plot )
971+ ggsave(filename = " PCoA_without_labels_GLAmpSeq.png" , plot = ordination_plot_u )
972972
973973```
974974
@@ -1009,7 +1009,7 @@ Run the DESeq() function to normalize for sample read-depth and composition, tra
10091009``` R
10101010deseq_modeled <- DESeq(deseq_obj )
10111011
1012- write.table(counts(deseq_modeled , normalized = TRUE ), file = paste0( " normalized_counts_GLAmpSeq.tsv" ) , sep = " \t " , row.names = TRUE , quote = FALSE )
1012+ write.table(counts(deseq_modeled , normalized = TRUE ), file = " normalized_counts_GLAmpSeq.tsv" , sep = " \t " , row.names = TRUE , quote = FALSE )
10131013```
10141014
10151015** Input Data:**
0 commit comments