``` r
# Calculate DESeq object
Peve_counts_sRNA_rowNames <- Peve_counts_sRNA 
#%>% column_to_rownames(var = "Name")

dds_Peve_sRNA <- DESeqDataSetFromMatrix(countData = Peve_counts_sRNA_rowNames,
                              colData = Peve_metadata_sRNA,
                              design = ~ 1) 

# Run differential expression analysis 
# (Note that this DESeq() function runs all necessary steps, including data normalization, 
# estimating size factors, estimating dispersions, gene-wise dispersion estimates, mean-dispersion 
# relationship, final dispersion estimates, fitting model, and testing)
# Using design = ~1 because we don't have treatment groups

dds_Peve_sRNA <- DESeq(dds_Peve_sRNA)
```