[R] heatmap2 error key

2016-04-19 Thread Catalina Aguilar Hurtado
Hi I am trying to understand what happen with the heatmap.2 code that it used to work (last used in October 2015). I am able to get a heatmap but my colour key doesn't come up and instead I get an error in all my files the used to work. Does anyone has any idea of what could have changed? or how to

[R] How to make a proper use of blocking in limma using voom

2014-04-06 Thread Catalina Aguilar Hurtado
Hi all, I have a RNAseq data to analyse were I have a control and a one treatment for different individuals. I need to block the effects of the individual, but I am having several troubles to get the data that I need. I am using voom because my data is very heterogeneous and voom seams to do a goo

[R] How to make a proper use of blocking in limma using voom

2014-04-10 Thread Catalina Aguilar Hurtado
Hi Steve, Thanks for you reply. I have tried what you suggested and checked the limma guide sessions. Now I am getting the results from only one treatment which is great but the problem is all my adjusted Pvalues (adj.P.Val) are >0.05. Don't really know how to solve this, I know my data works be

[R] How to make a proper use of blocking in limma using voom

2014-04-12 Thread Catalina Aguilar Hurtado
ta works because I used DESeq and got around 200 differentially expressed genes FDR < 0.05. >From: Steve Lianoglou >Date: Mon, Apr 7, 2014 at 5:09 PM >Subject: Re: [R] How to make a proper use of blocking in limma using >voom >To: Catalina Aguilar Hurtado >Cc: r help

[R] DESeq vs DESeq2 different DEGs results

2014-05-05 Thread Catalina Aguilar Hurtado
Hi, I want to compare DESeq vs DESeq2 and I am getting different number of DEGs which I will expect to be normal. However, when I compare the 149 genes ID that I get with DESeq with the 869 from DESeq2 there are only ~10 genes that are in common which I don’t understand (using FDR <0.05 for bot

[R] DESeq2 nbinomLRT reduced formula

2014-02-17 Thread Catalina Aguilar Hurtado
Hi, I am trying to use DESeq2 but I having troubles determining the reduced formula for the nbinomLRT. I will like to compare the results with the DEGs I got by nbinomGLMTest using DESeq. With DESeq I did the following: co=read.table("C_LPS_PBS_1h_DS.txt",header=T) ind=c(1,2,4,5,1,3,4,5) trt=c(