Open richfitz opened 10 years ago
Comment by mwpennell from Tuesday Dec 10, 2013 at 00:26 GMT
sorry i missed this one -- shouldn't have opened a new issue. i am not quite sure about why this is but the (n-(k-1))/n seems to be correct. i will change this.
Issue by richfitz from Wednesday Dec 04, 2013 at 05:39 GMT Originally opened as https://github.com/richfitz/modeladequacy/issues/56
I've taken a shot at this (see 79c2650), but I'm not sure that things are correct.
which all suggests that in
arbutus:::estimate.sigma2.gls()
we should use(n - (k - 1)) / n
as the scaling factor. This doesn't quite match my view of how these work, but I might be being dense.Ideas?
(this depends on very recent versions of both arbutus and diversitree)