[Robast-commits] r1176 - branches/robast-1.2/pkg/RobExtremes/R

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

Author: ruckdeschel
Date: 2019-03-01 11:12:51 +0100 (Fri, 01 Mar 2019)
New Revision: 1176

Modified:
   branches/robast-1.2/pkg/RobExtremes/R/getStartICPareto.R
   branches/robast-1.2/pkg/RobExtremes/R/internal-getpsi.R
Log:
[RobExtremes] branch 1.2 / asBias was (not) yet set -- now the following code works:

library(RobExtremes)
model <- GEVFamilyMuUnknown()
library(robustbase)
data(los)
system.time(RMX.fast <- roptest(los, risk = RMXRRisk(samplesize = 207),
                                L2Fam = model, eps.lower = 0,
                                eps.upper = 0.1, steps = 3))
RMX.fast
confint(RMX.fast)
## Fehler ...
confint(RMX.fast, symmetricBias())

Modified: branches/robast-1.2/pkg/RobExtremes/R/getStartICPareto.R
===================================================================
--- branches/robast-1.2/pkg/RobExtremes/R/getStartICPareto.R	2019-02-27 16:00:34 UTC (rev 1175)
+++ branches/robast-1.2/pkg/RobExtremes/R/getStartICPareto.R	2019-03-01 10:12:51 UTC (rev 1176)
@@ -63,9 +63,13 @@
    }else weight(w) <- minbiasweight(w, neighbor = nb, biastype = biast,
                           normW = normt)
 
+   Risk <- list(asBias = list(value = b, biastype = biast,
+                                       normtype = normt,
+                                       neighbortype = class(nb)))
+
    res <- list(a = a, A = A, b = b, d = 0*a,
                normtype = normt, biastype = biast, w = w,
-               info = c("optIC", ICT), risk = list(),
+               info = c("optIC", ICT), risk = Risk,
                modifyIC = NULL)
 
 

Modified: branches/robast-1.2/pkg/RobExtremes/R/internal-getpsi.R
===================================================================
--- branches/robast-1.2/pkg/RobExtremes/R/internal-getpsi.R	2019-02-27 16:00:34 UTC (rev 1175)
+++ branches/robast-1.2/pkg/RobExtremes/R/internal-getpsi.R	2019-03-01 10:12:51 UTC (rev 1176)
@@ -53,10 +53,13 @@
                           normW = normt)
    }else weight(w) <- minbiasweight(w, neighbor = nb, biastype = biast,
                           normW = normt)
+   Risk <- list(asBias = list(value = b, biastype = biast,
+                                       normtype = normt,
+                                       neighbortype = class(nb)))
 
    res <- list(a = a, A = A, b = b, d = 0*a,
                normtype = normt, biastype = biast, w = w,
-               info = c("optIC", ICT), risk = list(),
+               info = c("optIC", ICT), risk = Risk,
                modifyIC = NULL)
 
 
@@ -120,9 +123,13 @@
    }else weight(w) <- minbiasweight(w, neighbor = nb, biastype = biast,
                           normW = normt)
 
+   Risk <- list(asBias = list(value = b, biastype = biast,
+                                       normtype = normt,
+                                       neighbortype = class(nb)))
+
    res <- list(a = a, A = A, b = b, d = 0*a,
                normtype = normt, biastype = biast, w = w,
-               info = c("optIC", ICT), risk = list(),
+               info = c("optIC", ICT), risk = Risk,
                modifyIC = NULL)