Hi,
I'm using deseq2 for a RNAseq experiment with 2 genotypes and 5 conditions (2 replicates for each combination of genotype and condition). When I use the deseq command to perform a comparison, I get different results if all experiments are in the DESeqDataSet than when only the 2 compared experiments are in the DESeqDataSet, what is expected as the estimation steps for the model parameters look at all the samples in the DESeqDataSet. Then, if I try to apply lfcShrink with the DESeqDataSet with all experiments, it takes forever. I've waited for 24 hours and still no results. As I have many comparisons, this is not an option. Same thing with rlog. Is it normal? How should I proceed?
Thanks,
Thierry
I haven't tried on another machine, but if it works on yours, it means that mine has a problem. I fully reinstalled R and remove all libraries, reinstalled DESeq2, still the same... I suspect something more vicious, maybe related to GCC compiler or C++ libraries. I'll have a look at it next week, here everything is closing with Irma coming. I'll let you know when I find a solution, if I find one (hope so :).
Sorry, I didn't have the opportunity to work on it. Everything shut down here because of Irma.
I worked on a machine with Ubuntu 14.4. I tried to reinstall all packages and R, but nothing changed. I have the feeling that it's related to GCC as it says when I install DESeq2 that openMP is not working while I have gcc-4.8 and gcc-4.9 (default) installed on my computer, both enabling the use of openMP. I tried to find a way to tell R which gcc to use, but it's really not clear how it works and I couldn't find a way to do it properly. I have a virtualbox with Windows 10 so I installed R on it and I'm now able to run lfcShrink normally. It's great in the sense that I'm not stuck and I can analyze the data I have, but it's frustrating not to be able to understand why it's not working with Ubuntu. Do you have any idea?
Sorry, no ideas...
No worries, I found a way to make it work. I'll try it again when I upgrade to Ubuntu 16, never knows.
I upgraded to Ubuntu 16 and ... it works!!!!!!!!!
I can't tell why it works now, but it might be helpful for others. Michael, thanks a lot for your comments!