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

[noreply at r-forge.r-project.org]

Author: ruckdeschel
Date: 2014-04-03 13:07:10 +0200 (Thu, 03 Apr 2014)
New Revision: 739

Modified:
   branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R
Log:
...RobExtremes: now the sigma search range could include negative values -> fixed

Modified: branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R
===================================================================
--- branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R	2014-04-03 09:37:25 UTC (rev 738)
+++ branches/robast-1.0/pkg/RobExtremes/R/startEstGEV.R	2014-04-03 11:07:10 UTC (rev 739)
@@ -11,7 +11,7 @@
          mygev1 <- GEV(loc=0,scale=sig,shape=xi)
          CvMDist(x0,mygev1)
       }
-      intv <- if(xi<0) c(min(1e-5,xi/4*min(x)-mu), max(xi*min(x)-mu,2e-5)) else c(max(1e-5,xi*max(x)-mu),max(s0,2*xi*max(x)-mu))
+      intv <- if(xi<0) c(max(1e-5,xi/4*min(x)-mu), max(xi*min(x)-mu,2e-5)) else c(max(1e-5,xi*max(x)-mu),max(s0,2*xi*max(x)-mu))
       sigCvMMD1 <- optimize(funl, interval=intv)$minimum
       mygev <- GEVFamily(loc=0,scale=sigCvMMD1,shape=xi, withPos=withPos,
                          start0Est = fu, ..withWarningGEV=FALSE)