Quais picos são diferentes?
Vamos analisar o quanto as amostras diferem em termos de chamadas de pico. Observando os picos mesclados com os quais você trabalhou antes, quantos picos são compartilhados entre diferentes amostras?
Um gráfico que pode responder a essa pergunta é um diagrama de Venn. Eles podem ser criados diretamente a partir do objeto do DBA usando a função dba.plotVenn(). Uma alternativa é usar um gráfico UpSet, que você pode criar com o pacote UpSetR. Ele substitui as elipses de um diagrama de Venn por uma série de gráficos de barras, indicando o tamanho de vários subconjuntos. A função upset() exige que você forneça um data frame com uma coluna para cada conjunto, isto é, amostra, indicando a pertinência de cada pico ao conjunto. Os dados necessários estão no campo called do objeto ar_diff.
Este exercício faz parte do curso
ChIP-seq com Bioconductor em R
Exercício interativo prático
Experimente este exercício completando este código de exemplo.
# Create Venn diagram
___(ar_diff, mask=1:4)