I'm developing an R-based workflow using crlmm, which uses a normalization method different from what Illumina's proprietary BeadStudio software performs. I'm not sure exactly how it differs, and is it better or worse for my purposes? Specifically, I am starting with paired Grn/Red IDAT files for tumor-normal samples per patient. After probe-level copy-number estimation using crlmm, I'll do segmentation of the probes (using CBS), and identify somatic CNVs (subtract normal CNAs from the tumor). Those somatic CNVs will be fed into GISTIC to find recurrently/focally altered genes.
Any insight into my proposed workflow would be massively helpful. I'm sure that there might be a better solution. I would like to avoid using any proprietary software since this needs to be automated from start to end.