## <FIXME>
## maybe need to refine below in case of non-Brownian corStruct
if (!missing(corStruct)) phy <- attr(corStruct, "tree")
- dfP <- sum(phy$edge.length)*N / sum(diag(R))
+ dfP <- sum(phy$edge.length)*N / sum(diag(vcv(phy))) # need the variances
## </FIXME>
## compute QIC: