Hi, everyone, sorry if it s a stupid question.
My dataset is miRNA's, we are using DESeq to check a case versus control in two different condition.
Since my samples are all paired between conditions I want to estimate dispertion using the fact that I have two sample per subjects.
I am not sure if what I am doing is good, could you make me understand my own command ?
Thank you !!
My dataset is miRNA's, we are using DESeq to check a case versus control in two different condition.
Since my samples are all paired between conditions I want to estimate dispertion using the fact that I have two sample per subjects.
Code:
cds=estimateDispersions( cds, method='pooled-CR', sharingMode="maximum", fitType="local", modelFormula = count ~ miRNAdesign$subject ) cds2=cds[,is_condition1] exact test and blablabla ...
Thank you !!