[Robast-commits] r741 - branches/robast-1.0/pkg/RobExtremes/R

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Thu Apr 3 13:35:57 CEST 2014


Author: ruckdeschel
Date: 2014-04-03 13:35:57 +0200 (Thu, 03 Apr 2014)
New Revision: 741

Modified:
   branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R
Log:
argh: RobExtremes --- stills throws errors ... hopefully this is it

Modified: branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R
===================================================================
--- branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R	2014-04-03 11:21:14 UTC (rev 740)
+++ branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R	2014-04-03 11:35:57 UTC (rev 741)
@@ -18,9 +18,10 @@
       mde0 <- try(MDEstimator(x0, mygev, distance=CvMDist, startPar=c("scale"=sigCvMMD1,"shape"=xi)),silent=TRUE)
       if(!is(mde0,"try-error")){
           es <- estimate(mde0)
-          if(criterion(mde0)<crit0&(1+es[2]*min(x0)/es[1]>0)){
+          crit1 <- criterion(mde0)
+          if(crit1<crit0&(1+min(es[2]*x0/es[1])>0)){
              mdeb <- mde0
-             crit0 <- criterion(mde0)
+             crit0 <- crit1
           }
       }
   }



More information about the Robast-commits mailing list