My question is more naive. I am using ChIPseeker to compare multiple Peak files. How to read them all and store in "files" object.
files <- getSampleFiles()
Let say I have two peak files
How to read and calculate tagMatrix for my multiple ChIP-seq files.. I have sam, bam, bed format aligned files as well. How can I come to tagMatrixList from the files that I have to come to this point:
tagMatrixList <- lapply(files, getTagMatrix, windows=promoter)
Could you please show me the beginning codes/steps.