[Robast-commits] r740 - 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:21:15 CEST 2014


Author: ruckdeschel
Date: 2014-04-03 13:21:14 +0200 (Thu, 03 Apr 2014)
New Revision: 740

Modified:
   branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R
Log:
RobExtremes: now only let pass admissible starting estimators

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



More information about the Robast-commits mailing list