View Single Post
Old 04-30-2012, 10:06 AM   #4
Simon Anders
Senior Member
Location: Heidelberg, Germany

Join Date: Feb 2010
Posts: 993

First of all: The command "featureData(ecs)$dispersion_CR_est" did not work because you are reading an outdated manual. The slot "dispersion_CR_est" is now called "dispBeforeSharing".

Use the command "openVignette()" to get documentation that fits your installed package versions.

In the current manual, you will see instructions to plot the estimated dispersions against the mean count value. Have a look, maybe your data is simply to noisy to get any results. (If you unsure how to interpret the plot, just post it here.)

Also have a look at the count values plot for few genes (with plotDEXSeq(ecs, "my_favourite_geneID", norCounts = TRUE)). Does it look like you should detect something?
Simon Anders is offline   Reply With Quote