[Robast-commits] r1325 - in branches/robast-1.3/pkg/ROptEst: . inst tests/Examples
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Sun Jan 12 01:51:36 CET 2025
Author: ruckdeschel
Date: 2025-01-12 01:51:36 +0100 (Sun, 12 Jan 2025)
New Revision: 1325
Modified:
branches/robast-1.3/pkg/ROptEst/DESCRIPTION
branches/robast-1.3/pkg/ROptEst/inst/NEWS
branches/robast-1.3/pkg/ROptEst/tests/Examples/ROptEst-Ex.Rout.save
Log:
[ROptEst] ported changes from v 1.3.5 to branch 1.3
Modified: branches/robast-1.3/pkg/ROptEst/DESCRIPTION
===================================================================
--- branches/robast-1.3/pkg/ROptEst/DESCRIPTION 2025-01-12 00:50:47 UTC (rev 1324)
+++ branches/robast-1.3/pkg/ROptEst/DESCRIPTION 2025-01-12 00:51:36 UTC (rev 1325)
@@ -1,5 +1,5 @@
Package: ROptEst
-Version: 1.3.5
+Version: 1.3.7
Date: 2024-08-29
Title: Optimally Robust Estimation
Description: R infrastructure for optimally robust estimation in general smoothly
@@ -8,7 +8,7 @@
Rieder, H., Kohl, M., and Ruckdeschel, P. (2008), <doi:10.1007/s10260-007-0047-7>.
Depends: R(>= 3.4), methods, distr(>= 2.8.0), distrEx(>= 2.8.0), distrMod(>= 2.8.1),
RandVar(>= 1.2.0), RobAStBase(>= 1.2.0)
-Imports: startupmsg, MASS, stats, graphics, utils, grDevices
+Imports: startupmsg(>= 1.0.0), MASS, stats, graphics, utils, grDevices
Suggests: RobLox
Authors at R: c(person("Matthias", "Kohl", role=c("cre", "cph"),
email="Matthias.Kohl at stamats.de", comment = c(ORCID = "0000-0001-9514-8910")),
Modified: branches/robast-1.3/pkg/ROptEst/inst/NEWS
===================================================================
--- branches/robast-1.3/pkg/ROptEst/inst/NEWS 2025-01-12 00:50:47 UTC (rev 1324)
+++ branches/robast-1.3/pkg/ROptEst/inst/NEWS 2025-01-12 00:51:36 UTC (rev 1325)
@@ -8,6 +8,13 @@
information)
#######################################
+version 1.3.5
+#######################################
+
+under the hood:
++ adapted reference output for new startupmsg
+
+#######################################
version 1.3.4
#######################################
Modified: branches/robast-1.3/pkg/ROptEst/tests/Examples/ROptEst-Ex.Rout.save
===================================================================
--- branches/robast-1.3/pkg/ROptEst/tests/Examples/ROptEst-Ex.Rout.save 2025-01-12 00:50:47 UTC (rev 1324)
+++ branches/robast-1.3/pkg/ROptEst/tests/Examples/ROptEst-Ex.Rout.save 2025-01-12 00:51:36 UTC (rev 1325)
@@ -1,7 +1,7 @@
-R Under development (unstable) (2019-03-28 r76288) -- "Unsuffered Consequences"
-Copyright (C) 2019 The R Foundation for Statistical Computing
-Platform: x86_64-w64-mingw32/x64 (64-bit)
+R Under development (unstable) (2025-01-10 r87562 ucrt) -- "Unsuffered Consequences"
+Copyright (C) 2025 The R Foundation for Statistical Computing
+Platform: x86_64-w64-mingw32/x64
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
@@ -24,32 +24,29 @@
> library('ROptEst')
Loading required package: distr
Loading required package: startupmsg
-:startupmsg> Utilities for Start-Up Messages (version 0.9.6)
+:startupmsg> *** Utilities for Start-Up Messages ***
:startupmsg>
-:startupmsg> For more information see ?"startupmsg",
-:startupmsg> NEWS("startupmsg")
-
+:startupmsg> Version information in start-up messages is
+:startupmsg> currently suppressed. To see such information on
+:startupmsg> startup as in versions of this pkg prior to this
+:startupmsg> versionr, set option "StartupBanner" to a value
+:startupmsg> different to {"off", "no-version", "no -
+:startupmsg> version"}, e.g., by options("StartupBanner" =
+:startupmsg> "complete") or by options("StartupBanner" = NULL)
+:startupmsg> or by options("StartupBanner" = "something else").
+:startupmsg>
+:startupmsg> Detailed information about which packages are
+:startupmsg> currently loaded or attached at which version
+:startupmsg> (regardless of whether these have start-up
+:startupmsg> messages managed by this package) can be obtained
+:startupmsg> by "sessionInfo()".
+:startupmsg>
+:startupmsg>
Loading required package: sfsmisc
-:distr> Object Oriented Implementation of Distributions (version
-:distr> 2.8.0)
+:distr> *** Object Oriented Implementation of Distributions ***
:distr>
-:distr> Attention: Arithmetics on distribution objects are
-:distr> understood as operations on corresponding random variables
-:distr> (r.v.s); see distrARITH().
:distr>
-:distr> Some functions from package 'stats' are intentionally masked
-:distr> ---see distrMASK().
-:distr>
-:distr> Note that global options are controlled by distroptions()
-:distr> ---c.f. ?"distroptions".
-:distr>
-:distr> For more information see ?"distr", NEWS("distr"), as well as
-:distr> http://distr.r-forge.r-project.org/
-:distr> Package "distrDoc" provides a vignette to this package as
-:distr> well as to several extension packages; try
-:distr> vignette("distr").
-
Attaching package: 'distr'
The following objects are masked from 'package:stats':
@@ -57,23 +54,10 @@
df, qqplot, sd
Loading required package: distrEx
-:distrEx> Extensions of Package 'distr' (version 2.8.0)
+:distrEx> *** Extensions of Package 'distr' ***
:distrEx>
-:distrEx> Note: Packages "e1071", "moments", "fBasics" should be
-:distrEx> attached /before/ package "distrEx". See
-:distrEx> distrExMASK().Note: Extreme value distribution
-:distrEx> functionality has been moved to
:distrEx>
-:distrEx> package "RobExtremes". See distrExMOVED().
-:distrEx>
-:distrEx> For more information see ?"distrEx", NEWS("distrEx"), as
-:distrEx> well as
-:distrEx> http://distr.r-forge.r-project.org/
-:distrEx> Package "distrDoc" provides a vignette to this package
-:distrEx> as well as to several related packages; try
-:distrEx> vignette("distr").
-
Attaching package: 'distrEx'
The following objects are masked from 'package:stats':
@@ -82,36 +66,16 @@
Loading required package: distrMod
Loading required package: RandVar
-:RandVar> Implementation of Random Variables (version 1.2.0)
+:RandVar> *** Implementation of Random Variables ***
:RandVar>
-:RandVar> For more information see ?"RandVar", NEWS("RandVar"), as
-:RandVar> well as
-:RandVar> http://robast.r-forge.r-project.org/
-:RandVar> This package also includes a vignette; try
-:RandVar> vignette("RandVar").
-
+:RandVar>
Loading required package: MASS
Loading required package: stats4
-:distrMod> Object Oriented Implementation of Probability Models
-:distrMod> (version 2.8.1)
+:distrMod> *** Object Oriented Implementation of Probability
+:distrMod> Models ***
:distrMod>
-:distrMod> Some functions from pkg's 'base' and 'stats' are
-:distrMod> intentionally masked ---see distrModMASK().
:distrMod>
-:distrMod> Note that global options are controlled by
-:distrMod> distrModoptions() ---c.f. ?"distrModoptions".
-:distrMod>
-:distrMod> For more information see ?"distrMod",
-:distrMod> NEWS("distrMod"), as well as
-:distrMod> http://distr.r-forge.r-project.org/
-:distrMod> There is a vignette to this package; try
-:distrMod> vignette("distrMod").
-:distrMod> Package "distrDoc" provides a vignette to the other
-:distrMod> distrXXX packages,
-:distrMod> as well as to several related packages; try
-:distrMod> vignette("distr").
-
Attaching package: 'distrMod'
The following object is masked from 'package:stats4':
@@ -129,21 +93,12 @@
Loading required package: RobAStBase
Loading required package: rrcov
Loading required package: robustbase
-Scalable Robust Estimators with High Breakdown Point (version 1.4-7)
+Scalable Robust Estimators with High Breakdown Point (version 1.7-6)
-:RobAStBase> Robust Asymptotic Statistics (version 1.2.1)
+:RobAStBase> *** Robust Asymptotic Statistics ***
:RobAStBase>
-:RobAStBase> Some functions from pkg's 'stats' and 'graphics'
-:RobAStBase> are intentionally masked ---see RobAStBaseMASK().
:RobAStBase>
-:RobAStBase> Note that global options are controlled by
-:RobAStBase> RobAStBaseoptions() ---c.f. ?"RobAStBaseoptions".
-:RobAStBase>
-:RobAStBase> For more information see ?"RobAStBase",
-:RobAStBase> NEWS("RobAStBase"), as well as
-:RobAStBase> http://robast.r-forge.r-project.org/
-
Attaching package: 'RobAStBase'
The following object is masked from 'package:graphics':
@@ -454,20 +409,14 @@
>
> ### ** Examples
>
-> N0 <- NormLocationScaleFamily(mean=0, sd=1)
-> N0.Rob1 <- InfRobModel(center = N0,
-+ neighbor = ContNeighborhood(radius = 0.5))
+> ## all (interesting) examples to this function need
+> ## more time than 5 seconds;
+> ## you can find them in
+> ## system.file("scripts", "examples_taking_longer.R",
+> ## package="ROptEst")
>
-> ## Don't run to reduce check time on CRAN
-> ## Not run:
-> ##D IC1 <- optIC(model = N0, risk = asCov())
-> ##D IC2 <- optIC(model = N0.Rob1, risk = asMSE())
-> ##D
-> ##D comparePlot(IC1,IC2, withMBR=TRUE)
-> ## End(Not run)
>
>
->
> cleanEx()
> nameEx("getL1normL2deriv")
> ### * getL1normL2deriv
@@ -564,6 +513,8 @@
> ## taken from script NormalScaleModel.R in folder scripts
> N0 <- NormScaleFamily(mean=0, sd=1)
> (N0.IC7 <- radiusMinimaxIC(L2Fam=N0, neighbor=nb, risk=ri, loRad=0, upRad=Inf))
+Warning in .local(L2Fam, neighbor, risk, ...) :
+ Had to modify radius bounds to [ 0.00241827117512196 2.82842712474619 ] after 2 iterations.
An object of class "ContIC"
### name: IC of contamination type
@@ -652,8 +603,9 @@
>
> myrisk <- asMSE()
> getRiskFctBV(myrisk)
-function(bias, var) return(bias^2+var)
-<environment: 0x0000000012f4e020>
+function (bias, var)
+return(bias^2 + var)
+<environment: 0x000001ad11b91578>
>
>
>
@@ -739,10 +691,10 @@
function (e1, e2, ...)
standardGeneric("CvMDist")
-<bytecode: 0x000000000c1740b8>
-<environment: 0x000000000c140200>
+<bytecode: 0x000001acf6e5c968>
+<environment: 0x000001acf6d81778>
Methods may be defined for arguments: e1, e2
-Use showMethods("CvMDist") for currently available ones.
+Use showMethods(CvMDist) for currently available ones.
> gennbCtrl()
$neighbor
@@ -955,71 +907,14 @@
> #############################
> ## 3. Normal (Gaussian) location and scale
> #############################
-> ## 24 determinations of copper in wholemeal flour
-> library(MASS)
-> data(chem)
-> plot(chem, main = "copper in wholemeal flour", pch = 20)
>
-> ## Family
-> NF <- NormLocationScaleFamily()
-> ## ML-estimate
-> MLest <- MLEstimator(chem, NF)
-> estimate(MLest)
- mean sd
-4.280417 5.185859
-> confint(MLest)
-A[n] asymptotic (LAN-based) confidence interval:
- 2.5 % 97.5 %
-mean 2.205679 6.355154
-sd 3.718798 6.652920
-Type of estimator: Maximum likelihood estimate
-samplesize: 24
-Call by which estimate was produced:
-MLEstimator(x = chem, ParamFamily = NF)
+> ## this example of a two dimensional parameter
+> ## to be estimated will need more time than
+> ## 5 seconds to run
+> ## you can find it in
+> ## system.file("scripts", "examples_taking_longer.R",
+> ## package="ROptEst")
>
-> ## Don't run to reduce check time on CRAN
-> ## Not run:
-> ##D ## compute optimally robust estimator (known contamination)
-> ##D ## takes some time -> you can use package RobLox for normal
-> ##D ## location and scale which is optimized for speed
-> ##D nb1 <- gennbCtrl(eps = 0.05)
-> ##D robEst <- robest(chem, NF, nbCtrl = nb1, steps = 3)
-> ##D estimate.call(robEst)
-> ##D attr(robEst,"timings")
-> ##D estimate(robest)
-> ##D
-> ##D confint(robest, symmetricBias())
-> ##D plot(pIC(robest))
-> ##D ## plot of relative and absolute information; cf. Kohl (2005)
-> ##D infoPlot(pIC(robest))
-> ##D
-> ##D tmp <- qqplot(chem, robest, cex.pch=1.5, exp.cex2.pch = -.25,
-> ##D exp.fadcol.pch = .55, withLab = TRUE, which.Order=1:4,
-> ##D exp.cex2.lbl = .12,exp.fadcol.lbl = .45,
-> ##D nosym.pCI = TRUE, adj.lbl=c(1.7,.2),
-> ##D exact.pCI = FALSE, log ="xy")
-> ##D
-> ##D ## finite-sample correction
-> ##D if(require(RobLox)){
-> ##D n <- length(chem)
-> ##D r <- 0.05*sqrt(n)
-> ##D r.fi <- finiteSampleCorrection(n = n, r = r)
-> ##D fsCor0 <- r.fi/r
-> ##D nb1 <- gennbCtrl(eps = 0.05)
-> ##D robest <- robest(chem, NF, nbCtrl = nb1, fsCor = fsCor0, steps = 3)
-> ##D estimate(robest)
-> ##D }
-> ##D
-> ##D ## compute optimally robust estimator (unknown contamination)
-> ##D ## takes some time -> use package RobLox!
-> ##D nb2 <- gennbCtrl(eps.lower = 0.05, eps.upper = 0.1)
-> ##D robest1 <- robest(chem, NF, nbCtrl = nb2, steps = 3)
-> ##D estimate(robest1)
-> ##D confint(robest1, symmetricBias())
-> ##D plot(pIC(robest1))
-> ##D ## plot of relative and absolute information; cf. Kohl (2005)
-> ##D infoPlot(pIC(robest1))
-> ## End(Not run)
>
>
>
@@ -1116,27 +1011,14 @@
> #############################
> ## 3. Normal (Gaussian) location and scale
> #############################
-> ## 24 determinations of copper in wholemeal flour
-> library(MASS)
-> data(chem)
-> plot(chem, main = "copper in wholemeal flour", pch = 20)
>
-> ## ML-estimate
-> MLest <- MLEstimator(chem, NormLocationScaleFamily())
-> estimate(MLest)
- mean sd
-4.280417 5.185859
-> confint(MLest)
-A[n] asymptotic (LAN-based) confidence interval:
- 2.5 % 97.5 %
-mean 2.205679 6.355154
-sd 3.718798 6.652920
-Type of estimator: Maximum likelihood estimate
-samplesize: 24
-Call by which estimate was produced:
-MLEstimator(x = chem, ParamFamily = NormLocationScaleFamily())
+> ## this example of a two dimensional parameter
+> ## to be estimated will need more time than
+> ## 5 seconds to run
+> ## you can find it in
+> ## system.file("scripts", "examples_taking_longer.R",
+> ## package="ROptEst")
>
-> ## Don't run to reduce check time on CRAN
>
>
>
@@ -1145,7 +1027,7 @@
> cleanEx()
> options(digits = 7L)
> base::cat("Time elapsed: ", proc.time() - base::get("ptime", pos = 'CheckExEnv'),"\n")
-Time elapsed: 10.58 0.46 11.05 NA NA
+Time elapsed: 8.93 0.28 9.2 NA NA
> grDevices::dev.off()
null device
1
More information about the Robast-commits
mailing list