same as before

R/gevdfit.R

@@ -95,6 +95,7 @@ gev.d.fit<-
     if(!is.null(init.vals)){init.vals <- data.frame(mu = init.vals[1], sigma = init.vals[2], xi = init.vals[3]
                                                     ,theta = init.vals[4], eta = init.vals[5])}
     if(any(is.null(c(muinit,siginit,shinit,etainit)))& is.null(init.vals)){
+      message('Initial values are calculated.')
       init.vals <- gev.d.init(xdat,ds,ifelse(is.null(thetainit),0,thetainit[1]))
     }
     # TODO: transform initial values with link function
@@ -283,7 +284,7 @@ gev.d.init <- function(xdat,ds,thetainit){
   # mean of shape parameters 
   shinit <- mean(mles[,3],na.rm = TRUE)
   
-  return(list(mu=muinit,sigma=siginit,xi=shinit,eta=etainit))
+  return(list(mu=muinit,sigma=siginit,xi=shinit,theta=thetainit,eta=etainit))
 }