Using greml I get snp heritability for binary trait of 0.029642 with standard error of 0.006261.

I know I can convert SNP heritability to liability scale like this:

#K=pop prevalence
#P=proportion of cases in study
#hsq=Heritability estimate (on observed scale)
#bigT = liability threshold
#tau = density of gaussian
K=my pop prev.
P=my proportion of cases
h2=0.029642
zv <- dnorm(qnorm(K))

h2_liab <- h2 * K^2 * ( 1 - K)^2 / P / (1-P) / zv^2

But how do I get standard error on the liability scale or is it just the same as the SNP heritability standard error. Thanks!