this is newversion. Updated
1 dispersion based gamma regression
2 estimating sigma^2 from DispSeq