use mean (instead of median) for etainit

R/gevdfit.R

@@ -278,7 +278,7 @@ gev.d.init <- function(xdat,ds,link){
   # sig0 <- exp Intercept
   siginit <- link$siglink$linkfun(exp(lmsig$coefficients[[1]]))
   # eta <- mean of negativ slopes 
-  etainit <- link$etalink$linkfun(median(c(-lmsig$coefficients[[2]],-lmmu$coefficients[[2]])))
+  etainit <- link$etalink$linkfun(mean(c(-lmsig$coefficients[[2]],-lmmu$coefficients[[2]])))
   # mean of mu_d/sig_d 
   # could try:
   # mu0/sig0 = exp(lmmu$coefficients[[1]])/exp(lmsig$coefficients[[1]])