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 exercicio faz parte do curso
ChIP-seq com Bioconductor em R
exercicio interativo prático
Tente este exercicio completando este código de exemplo.
# Create Venn diagram
___(ar_diff, mask=1:4)