Hi all,

after I got a dataset analyzed with this archaic ;-) NFCORE-atacseq ( https://nf-co.re/atacseq/2.1.2 ) , I've been asked to look for differences according to the sex of individuals using diffbind ( https://bioconductor.org/packages/release/bioc/html/DiffBind.html ) and I must say, I've no idea of what I'm doing.

So, I've got this samplesheet containing the sample names (Tissue+sex) , the path to the bams and the peaks, I'm not really sure about those "Tissue Condition QC Factor Replicate".

I use the R script below to run diffbind inherited from a colleague who left the lab. Everytime I try to run it, the process is killed even with the following slurm directives:

#SBATCH --job-name=diffbind
#SBATCH --nodes=1
#SBATCH --ntasks=1
#SBATCH --cpus-per-task=1
#SBATCH --mem-per-cpu=50G

so my questions are:

1) Is there anything I should change to get the diffbind by sex ?

2) how can I avoid this dba.count error ? I got the same error with bParallel=TRUE

Thanks