Dear all,

I'm currently using edgeR to detect the differentially expressed genes from a RNAseq datasets, and I'm also using the gof() function to test for potential outliers. I have two questions regarding the outlier detection, and would like to have your suggestions. 

1) How the outlier is defined? Is it the gene that have a deviance larger than a threshold? How is the deviance contained in the glmfit data calculated?

2) In gof() function, it assumes the deviance should follow a chi-squared distribution. But what is the statistic basis for this assumption? 

Thanks!

Yuan
	[[alternative HTML version deleted]]

