``` r
# extract normalized counts
# (normalization is automatically performed by deseq2)
Apul_counts_miRNA_norm <- counts(dds_Apul_miRNA, normalized=TRUE) %>% data.frame()

write.table(Apul_counts_miRNA_norm, file = "../output/03.10-D-Apul-sRNAseq-expression-DESeq2/Apul_counts_miRNA_normalized.txt", sep = "\t", row.names = TRUE, col.names = TRUE, quote = FALSE)


# variance stabilized data
vsd_Apul_miRNA <- varianceStabilizingTransformation(dds_Apul_miRNA, blind=TRUE)
wpn_vsd_Apul_miRNA <- getVarianceStabilizedData(dds_Apul_miRNA)
rv_wpn_Apul_miRNA <- rowVars(wpn_vsd_Apul_miRNA, useNames=TRUE)

Apul_counts_miRNA_vsd <- data.frame(wpn_vsd_Apul_miRNA)
write.table(Apul_counts_miRNA_vsd, file = "../output/03.10-D-Apul-sRNAseq-expression-DESeq2/Apul_counts_miRNA_variancestabilized.txt", sep = "\t", row.names = TRUE, col.names = TRUE,quote = FALSE)
```