I am analyzing paired samples (RNA-seq). I am facing problem with edgeR at
y <- estimateGLMTrendedDisp(d, design) step where it is saying error
message "Error in t(y) + prior.count.scaled : non-conformable arrays"

although it was ok with the step y <- estimateGLMCommonDisp(d, design)

Apparently in the assignment step something is happening [design <-
model.matrix (~Patient+Tissue)] and when I am trying to see whats there in
store at Design, its showing the matrix leaving the one of the sample
column absent.

Is that exactly what I should expect? Or this is the reason behind the
error I am facing. Any suggestion is welcome.

Thanks



Navonil De Sarkar
Senior Research Fellow
Human Genetics Unit
Indian Statistical Institute
203 B.T. Road Kolkata 700108
Mob:+91 9088153412
Skype. navonil.de.sarkar
E.mail ID. navonil.de@gmail.com

"I am all in favor of keeping dangerous weapons out of hands of Dumboes
Lets start with" Key board

	[[alternative HTML version deleted]]

