@@ -818,8 +818,8 @@ ASV_physeq <- phyloseq(count_tab_phy, tax_tab_phy, sample_info_tab_phy)
818
818
richness_and_diversity_estimates_by_sample <- plot_richness(ASV_physeq , color = " groups" , measures = c(" Chao1" , " Shannon" ))
819
819
richness_and_diversity_estimates_by_group <- plot_richness(ASV_physeq , x = " groups" , color = " groups" , measures = c(" Chao1" , " Shannon" ))
820
820
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 )
823
823
```
824
824
825
825
** Parameter Definitions:**
@@ -861,10 +861,10 @@ relative_classes <- plot_bar(proportions_physeq, x = "groups", fill = "class")
861
861
samplewise_phyla <- plot_bar(proportions_physeq , fill = " phylum" )
862
862
samplewise_classes <- plot_bar(proportions_physeq , fill = " class" )
863
863
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 )
868
868
```
869
869
870
870
** Input Data:**
@@ -967,8 +967,8 @@ ordination_plot_u <- plot_ordination(vst_physeq, vst_pcoa, color = "groups") +
967
967
annotate(" text" , x = Inf , y = - Inf , label = paste(" R2:" , toString(round(r2_value , 3 ))), hjust = 1.1 , vjust = - 2 , size = 4 )+
968
968
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" )
969
969
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 )
972
972
973
973
```
974
974
@@ -1009,7 +1009,7 @@ Run the DESeq() function to normalize for sample read-depth and composition, tra
1009
1009
``` R
1010
1010
deseq_modeled <- DESeq(deseq_obj )
1011
1011
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 )
1013
1013
```
1014
1014
1015
1015
** Input Data:**
0 commit comments