Commit e88d98b3 by Felix Fauer

### fixed fit. Now, diagnosis plots look good

parent 9d84b444
 ... @@ -298,7 +298,7 @@ gev.d.fit<- ... @@ -298,7 +298,7 @@ gev.d.fit<- mut.d <- mut*(sc0/((ds+theta)^eta )+tau) # new mut.d <- mut*(sc0/((ds+theta)^eta )+tau) # new #z\$data <- - log(as.vector((1 + xi * (xdat/sc.d-mut))^(-1/xi))) # old #z\$data <- - log(as.vector((1 + xi * (xdat/sc.d-mut))^(-1/xi))) # old z\$data <- - log(as.vector((1 + xi * ((xdat-mut)/sc.d))^(-1/xi))) # new z\$data <- - log(as.vector((1 + xi * ((xdat-mut.d)/sc.d))^(-1/xi))) # new z\$mle <- x\$par z\$mle <- x\$par test <- try( # catch error test <- try( # catch error z\$cov <- solve(x\$hessian) # invert hessian to get estimation on var-covar-matrix z\$cov <- solve(x\$hessian) # invert hessian to get estimation on var-covar-matrix ... ...
 ... @@ -66,7 +66,7 @@ plot(ann.max\$ds,ann.max\$xdat,log='xy',xlab = 'Duration [h]',ylab='Intensity [mm/ ... @@ -66,7 +66,7 @@ plot(ann.max\$ds,ann.max\$xdat,log='xy',xlab = 'Duration [h]',ylab='Intensity [mm/ * Step 2: fit d-GEV to annual maxima * Step 2: fit d-GEV to annual maxima ```{r fit} ```{r fit} fit <- gev.d.fit(xdat = ann.max\$xdat,ds = ann.max\$ds,sigma0link = make.link('log'), tau_zero = TRUE) fit <- gev.d.fit(xdat = ann.max\$xdat,ds = ann.max\$ds,sigma0link = make.link('log')) # checking the fit # checking the fit gev.d.diag(fit,pch=1,) gev.d.diag(fit,pch=1,) # parameter estimates # parameter estimates ... ...
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!