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

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


# variance stabilized data
vsd_Apul_siRNA <- varianceStabilizingTransformation(dds_Apul_siRNA, blind=TRUE)
wpn_vsd_Apul_siRNA <- getVarianceStabilizedData(dds_Apul_siRNA)
rv_wpn_Apul_siRNA <- rowVars(wpn_vsd_Apul_siRNA, useNames=TRUE)

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