Estimate the hyper-parameter by calling the gls function. The errors are allowed to be correlated and/or have unequal variances.
EstHyper(y, D, init.val)a vector of estimated hyper-parameter values plus log likelihood.
a vector of z-values.
a distance matrix defined based on the prior structure. Diagonal have to be zeros.
initial values for the transformed hyper-parameter. Default is 0.
Jun Chen
Xiao, J., Cao, H., Chen, J. (2017). False discovery rate control incorporating phylogenetic tree increases detection power in microbiome-wide multiple testing. Bioinformatics, 33(18), 2873-2881.