# program to calculate the h2v and all the GCV
# assuming only one random effect in the model; more complicated equations are presented in Sae-Lim 2015 and in Mulder 2016 to calculate also the approximate se of these parameters
varav_exp=0.891085E-01
vare=
varp=1.0
vare_exp=vare/exp(0.5*varav_exp)
varav_add=vare_exp^2*exp(2*varav_exp)-vare^2
h2v=varav_add/(2*varp*varp+3*varav_add)
GCV_Ve=sqrt(varav_exp)
GCV_SDe=0.5*GCV_Ve
GCV_Vp=GCV_Ve*vare/varp
GCV_SDp=GCV_SDe*vare/varp