Programming the statistical procedures from SAS

glimmix: dispersion estimates for small samples

Reply
Occasional Contributor
Posts: 6

glimmix: dispersion estimates for small samples

I am interested in the topic of dispersion estimates for count data coming from RNAseq data. The glimmix procedure provides the option empirical=mbn to correct for small sample size,and it appears to have good statistical properties in terms of type I error rates, however there is also the possibility to adjust the denumerator degrees of freedom for the F-tests with the KenwardRoger method. Does anybody have an idea what the best choice is? In our experiments, we typically only have 3 biological replicates for each treatment.

Respected Advisor
Posts: 2,655

Re: glimmix: dispersion estimates for small samples

I would go with the empirical= option for this, unless you have a lot of missing data.  If you are on SAS/STAT 12.1 or higher, consider ddfm=KR2, the improved Kenward-Rogers adjustment.

Steve Denham

Ask a Question
Discussion stats
  • 1 reply
  • 197 views
  • 0 likes
  • 2 in conversation