Counts normaliseren met DESeq2
We hebben het DESeq2-object aangemaakt en willen nu kwaliteitscontrole uitvoeren op onze monsters. Daarvoor moeten we de genormaliseerde counts genereren (genormaliseerd voor bibliotheekgrootte — het totale aantal gen-tellingen per monster — met inachtneming van de bibliotheeksamenstelling). Gebruik het DESeq2-object om de genormaliseerde count-matrix te verkrijgen.
Deze oefening maakt deel uit van de cursus
RNA-Seq met Bioconductor in R
Oefeninstructies
Schat de size factors voor de smoc2-countdata met de functie
estimateSizeFactors()en sla ze terug op in het objectdds_smoc2Extraheer de genormaliseerde countwaarden uit
dds_smoc2en sla ze op alssmoc2_normalized_countsmet de functiecounts().
Praktische interactieve oefening
Probeer deze oefening eens door deze voorbeeldcode in te vullen.
# Determine the size factors to use for normalization
dds_smoc2 <- ___(___)
# Extract the normalized counts
smoc2_normalized_counts <- ___(___, ___)