[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