Vilka peaks skiljer sig åt?
Låt oss titta på hur mycket proverna skiljer sig åt i fråga om peak-anrop. Med utgångspunkt i de sammanslagna peaks du arbetat med tidigare – hur många peaks delas mellan olika prover?
Ett diagram som kan besvara den frågan är ett Venn-diagram. Det kan skapas direkt från DBA-objektet med funktionen dba.plotVenn(). Ett alternativ är att använda ett UpSet-diagram, vilket du kan skapa med paketet UpSetR. Det ersätter ellipserna i ett Venn-diagram med en serie stapeldiagram som visar storleken på olika delmängder. Funktionen upset() kräver att du anger en dataram med en kolumn för varje mängd, det vill säga varje prov, som anger mängdtillhörighet för varje peak. Den data som krävs finns i fältet called i objektet ar_diff.
Den här övningen är en del av kursen
ChIP-seq med Bioconductor i R
Interaktiv övning med praktiskt arbete
Testa den här övningen genom att slutföra den här exempelkoden.
# Create Venn diagram
___(ar_diff, mask=1:4)