[Robast-commits] r1054 - branches/robast-1.1/pkg/ROptEst/R branches/robast-1.1/pkg/RobAStBase/R branches/robast-1.2/pkg/ROptEst/R branches/robast-1.2/pkg/RobAStBase/R pkg/ROptEst/R pkg/RobAStBase/R
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Tue Jul 24 15:05:38 CEST 2018
Author: ruckdeschel
Date: 2018-07-24 15:05:37 +0200 (Tue, 24 Jul 2018)
New Revision: 1054
Modified:
branches/robast-1.1/pkg/ROptEst/R/getInfLM.R
branches/robast-1.1/pkg/ROptEst/R/getInfRobIC_asAnscombe.R
branches/robast-1.1/pkg/ROptEst/R/getRiskIC.R
branches/robast-1.1/pkg/RobAStBase/R/getBiasIC.R
branches/robast-1.1/pkg/RobAStBase/R/getRiskIC.R
branches/robast-1.2/pkg/ROptEst/R/getInfLM.R
branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asAnscombe.R
branches/robast-1.2/pkg/ROptEst/R/getRiskIC.R
branches/robast-1.2/pkg/RobAStBase/R/getBiasIC.R
branches/robast-1.2/pkg/RobAStBase/R/getRiskIC.R
pkg/ROptEst/R/getInfLM.R
pkg/ROptEst/R/getInfRobIC_asAnscombe.R
pkg/ROptEst/R/getRiskIC.R
pkg/RobAStBase/R/getBiasIC.R
pkg/RobAStBase/R/getRiskIC.R
Log:
[RobAStBase], [ROptEst]
+ buglets in getRiskIC, getBiasIC -> lacked treatment of case that default of withCheck is not transmitted in method dispatch
+ in getLagrangeMultByIter and in "getInfRobIC", signature(L2deriv = "UnivariateDistribution", risk = "asAnscombe") catch NA return values
Modified: branches/robast-1.1/pkg/ROptEst/R/getInfLM.R
===================================================================
--- branches/robast-1.1/pkg/ROptEst/R/getInfLM.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.1/pkg/ROptEst/R/getInfLM.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -77,7 +77,7 @@
cat("current precision in IC algo:\t", prec, "\n")
print(round(c(A=prec[1],a=prec[2]),3))
}
- if(prec < tol) break
+ if(!is.na(prec)) if(prec < tol) break
if(iter > maxiter){
if(warnit)
cat("maximum iterations reached!\n",
Modified: branches/robast-1.1/pkg/ROptEst/R/getInfRobIC_asAnscombe.R
===================================================================
--- branches/robast-1.1/pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.1/pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -58,7 +58,7 @@
tol = toli, warn, noLow = noLow,
verbose = FALSE, checkBounds = FALSE)
trV <- upBerg$risk$trAsCov$value
- e.up <- FI/trV
+ if(!is.na(trV)) e.up <- FI/trV
upper <- upper * 3
}
upper <- upper / 3
Modified: branches/robast-1.1/pkg/ROptEst/R/getRiskIC.R
===================================================================
--- branches/robast-1.1/pkg/ROptEst/R/getRiskIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.1/pkg/ROptEst/R/getRiskIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -7,6 +7,7 @@
L2Fam = "missing"),
function(IC, risk, withCheck= TRUE){
L2Fam <- force(eval(IC at CallL2Fam))
+ if(missing(withCheck)) withCheck <- TRUE
getRiskIC(IC = IC, risk = risk, L2Fam = L2Fam, withCheck= withCheck)
})
@@ -15,6 +16,7 @@
neighbor = "missing",
L2Fam = "L2ParamFamily"),
function(IC, risk, L2Fam, withCheck= TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
Cov <- eval(IC at Risks[["asCov"]])
if(is.null(Cov)){
if(numberOfMaps(L2Fam at L2deriv)==1){ ## L2derivDim <- L2Fam at L2deriv
Modified: branches/robast-1.1/pkg/RobAStBase/R/getBiasIC.R
===================================================================
--- branches/robast-1.1/pkg/RobAStBase/R/getBiasIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.1/pkg/RobAStBase/R/getBiasIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -12,6 +12,8 @@
misF <- TRUE
L2Fam <- eval(IC at CallL2Fam)
}
+ if(missing(withCheck)) withCheck <- TRUE
+
D1 <- L2Fam at distribution
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(D1)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
Modified: branches/robast-1.1/pkg/RobAStBase/R/getRiskIC.R
===================================================================
--- branches/robast-1.1/pkg/RobAStBase/R/getRiskIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.1/pkg/RobAStBase/R/getRiskIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -16,9 +16,11 @@
risk = "asCov",
neighbor = "missing",
L2Fam = "missing"),
- function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE)
- getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
+ function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
+ return(getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
tol = tol, withCheck = withCheck))
+ })
setMethod("getRiskIC", signature(IC = "IC",
risk = "asCov",
@@ -28,6 +30,7 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
IC1 <- as(diag(dimension(IC at Curve)) %*% IC at Curve, "EuclRandVariable")
bias <- E(L2Fam, IC1)
@@ -46,8 +49,9 @@
neighbor = "missing",
L2Fam = "missing"),
function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
- tol = tol, withCheck = withCheck)
+ if(missing(withCheck)) withCheck <- TRUE
+ return(getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
+ tol = tol, withCheck = withCheck))
})
setMethod("getRiskIC", signature(IC = "IC",
@@ -58,6 +62,8 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
+
trCov <- getRiskIC(IC, risk = asCov(), L2Fam = L2Fam, withCheck = withCheck)$asCov
trCov$value <- sum(diag(as.matrix(trCov$value)))
@@ -73,7 +79,10 @@
neighbor = "UncondNeighborhood",
L2Fam = "missing"),
function(IC, risk, neighbor, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getBiasIC(IC = IC, neighbor = neighbor,
+
+ if(missing(withCheck)) withCheck <- TRUE
+
+ getBiasIC(IC = IC, neighbor = neighbor,
biastype = biastype(risk), normtype = normtype(risk), tol = tol,
withCheck = withCheck)
})
@@ -82,7 +91,8 @@
neighbor = "UncondNeighborhood",
L2Fam = "L2ParamFamily"),
function(IC, risk, neighbor, L2Fam, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getBiasIC(IC = IC, neighbor = neighbor, L2Fam = L2Fam,
+ if(missing(withCheck)) withCheck <- TRUE
+ getBiasIC(IC = IC, neighbor = neighbor, L2Fam = L2Fam,
biastype = biastype(risk), normtype = normtype(risk),
tol = tol, withCheck = withCheck)
})
@@ -94,6 +104,7 @@
neighbor = "UncondNeighborhood",
L2Fam = "missing"),
function(IC, risk, neighbor, tol = .Machine$double.eps^0.25, withCheck = TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
L2Fam <- eval(IC at CallL2Fam)
getRiskIC(IC = IC, risk = risk, neighbor = neighbor,
L2Fam = L2Fam, tol = tol, withCheck = withCheck)
@@ -107,6 +118,7 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
rad <- neighbor at radius
if(rad == Inf) return(Inf)
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfLM.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfLM.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfLM.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -77,7 +77,7 @@
cat("current precision in IC algo:\t", prec, "\n")
print(round(c(A=prec[1],a=prec[2]),3))
}
- if(prec < tol) break
+ if(!is.na(prec)) if(prec < tol) break
if(iter > maxiter){
if(warnit)
cat("maximum iterations reached!\n",
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asAnscombe.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -58,7 +58,7 @@
tol = toli, warn, noLow = noLow,
verbose = FALSE, checkBounds = FALSE)
trV <- upBerg$risk$trAsCov$value
- e.up <- FI/trV
+ if(!is.na(trV)) e.up <- FI/trV
upper <- upper * 3
}
upper <- upper / 3
Modified: branches/robast-1.2/pkg/ROptEst/R/getRiskIC.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getRiskIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.2/pkg/ROptEst/R/getRiskIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -7,6 +7,7 @@
L2Fam = "missing"),
function(IC, risk, withCheck= TRUE){
L2Fam <- force(eval(IC at CallL2Fam))
+ if(missing(withCheck)) withCheck <- TRUE
getRiskIC(IC = IC, risk = risk, L2Fam = L2Fam, withCheck= withCheck)
})
@@ -15,6 +16,7 @@
neighbor = "missing",
L2Fam = "L2ParamFamily"),
function(IC, risk, L2Fam, withCheck= TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
Cov <- eval(IC at Risks[["asCov"]])
if(is.null(Cov)){
if(numberOfMaps(L2Fam at L2deriv)==1){ ## L2derivDim <- L2Fam at L2deriv
Modified: branches/robast-1.2/pkg/RobAStBase/R/getBiasIC.R
===================================================================
--- branches/robast-1.2/pkg/RobAStBase/R/getBiasIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.2/pkg/RobAStBase/R/getBiasIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -12,6 +12,8 @@
misF <- TRUE
L2Fam <- eval(IC at CallL2Fam)
}
+ if(missing(withCheck)) withCheck <- TRUE
+
D1 <- L2Fam at distribution
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(D1)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
Modified: branches/robast-1.2/pkg/RobAStBase/R/getRiskIC.R
===================================================================
--- branches/robast-1.2/pkg/RobAStBase/R/getRiskIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ branches/robast-1.2/pkg/RobAStBase/R/getRiskIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -16,9 +16,11 @@
risk = "asCov",
neighbor = "missing",
L2Fam = "missing"),
- function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE)
- getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
+ function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
+ return(getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
tol = tol, withCheck = withCheck))
+ })
setMethod("getRiskIC", signature(IC = "IC",
risk = "asCov",
@@ -28,6 +30,7 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
IC1 <- as(diag(dimension(IC at Curve)) %*% IC at Curve, "EuclRandVariable")
bias <- E(L2Fam, IC1)
@@ -46,8 +49,9 @@
neighbor = "missing",
L2Fam = "missing"),
function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
- tol = tol, withCheck = withCheck)
+ if(missing(withCheck)) withCheck <- TRUE
+ return(getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
+ tol = tol, withCheck = withCheck))
})
setMethod("getRiskIC", signature(IC = "IC",
@@ -58,6 +62,8 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
+
trCov <- getRiskIC(IC, risk = asCov(), L2Fam = L2Fam, withCheck = withCheck)$asCov
trCov$value <- sum(diag(as.matrix(trCov$value)))
@@ -73,7 +79,10 @@
neighbor = "UncondNeighborhood",
L2Fam = "missing"),
function(IC, risk, neighbor, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getBiasIC(IC = IC, neighbor = neighbor,
+
+ if(missing(withCheck)) withCheck <- TRUE
+
+ getBiasIC(IC = IC, neighbor = neighbor,
biastype = biastype(risk), normtype = normtype(risk), tol = tol,
withCheck = withCheck)
})
@@ -82,7 +91,8 @@
neighbor = "UncondNeighborhood",
L2Fam = "L2ParamFamily"),
function(IC, risk, neighbor, L2Fam, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getBiasIC(IC = IC, neighbor = neighbor, L2Fam = L2Fam,
+ if(missing(withCheck)) withCheck <- TRUE
+ getBiasIC(IC = IC, neighbor = neighbor, L2Fam = L2Fam,
biastype = biastype(risk), normtype = normtype(risk),
tol = tol, withCheck = withCheck)
})
@@ -94,6 +104,7 @@
neighbor = "UncondNeighborhood",
L2Fam = "missing"),
function(IC, risk, neighbor, tol = .Machine$double.eps^0.25, withCheck = TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
L2Fam <- eval(IC at CallL2Fam)
getRiskIC(IC = IC, risk = risk, neighbor = neighbor,
L2Fam = L2Fam, tol = tol, withCheck = withCheck)
@@ -107,6 +118,7 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
rad <- neighbor at radius
if(rad == Inf) return(Inf)
Modified: pkg/ROptEst/R/getInfLM.R
===================================================================
--- pkg/ROptEst/R/getInfLM.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ pkg/ROptEst/R/getInfLM.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -77,7 +77,7 @@
cat("current precision in IC algo:\t", prec, "\n")
print(round(c(A=prec[1],a=prec[2]),3))
}
- if(prec < tol) break
+ if(!is.na(prec)) if(prec < tol) break
if(iter > maxiter){
if(warnit)
cat("maximum iterations reached!\n",
Modified: pkg/ROptEst/R/getInfRobIC_asAnscombe.R
===================================================================
--- pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -58,7 +58,7 @@
tol = toli, warn, noLow = noLow,
verbose = FALSE, checkBounds = FALSE)
trV <- upBerg$risk$trAsCov$value
- e.up <- FI/trV
+ if(!is.na(trV)) e.up <- FI/trV
upper <- upper * 3
}
upper <- upper / 3
Modified: pkg/ROptEst/R/getRiskIC.R
===================================================================
--- pkg/ROptEst/R/getRiskIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ pkg/ROptEst/R/getRiskIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -7,6 +7,7 @@
L2Fam = "missing"),
function(IC, risk, withCheck= TRUE){
L2Fam <- force(eval(IC at CallL2Fam))
+ if(missing(withCheck)) withCheck <- TRUE
getRiskIC(IC = IC, risk = risk, L2Fam = L2Fam, withCheck= withCheck)
})
@@ -15,6 +16,7 @@
neighbor = "missing",
L2Fam = "L2ParamFamily"),
function(IC, risk, L2Fam, withCheck= TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
Cov <- eval(IC at Risks[["asCov"]])
if(is.null(Cov)){
if(numberOfMaps(L2Fam at L2deriv)==1){ ## L2derivDim <- L2Fam at L2deriv
Modified: pkg/RobAStBase/R/getBiasIC.R
===================================================================
--- pkg/RobAStBase/R/getBiasIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ pkg/RobAStBase/R/getBiasIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -12,6 +12,8 @@
misF <- TRUE
L2Fam <- eval(IC at CallL2Fam)
}
+ if(missing(withCheck)) withCheck <- TRUE
+
D1 <- L2Fam at distribution
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(D1)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
Modified: pkg/RobAStBase/R/getRiskIC.R
===================================================================
--- pkg/RobAStBase/R/getRiskIC.R 2018-07-24 12:36:46 UTC (rev 1053)
+++ pkg/RobAStBase/R/getRiskIC.R 2018-07-24 13:05:37 UTC (rev 1054)
@@ -16,9 +16,11 @@
risk = "asCov",
neighbor = "missing",
L2Fam = "missing"),
- function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE)
- getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
+ function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
+ return(getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
tol = tol, withCheck = withCheck))
+ })
setMethod("getRiskIC", signature(IC = "IC",
risk = "asCov",
@@ -28,6 +30,7 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
IC1 <- as(diag(dimension(IC at Curve)) %*% IC at Curve, "EuclRandVariable")
bias <- E(L2Fam, IC1)
@@ -46,8 +49,9 @@
neighbor = "missing",
L2Fam = "missing"),
function(IC, risk, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
- tol = tol, withCheck = withCheck)
+ if(missing(withCheck)) withCheck <- TRUE
+ return(getRiskIC(IC = IC, risk = risk, L2Fam = eval(IC at CallL2Fam),
+ tol = tol, withCheck = withCheck))
})
setMethod("getRiskIC", signature(IC = "IC",
@@ -58,6 +62,8 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
+
trCov <- getRiskIC(IC, risk = asCov(), L2Fam = L2Fam, withCheck = withCheck)$asCov
trCov$value <- sum(diag(as.matrix(trCov$value)))
@@ -73,7 +79,10 @@
neighbor = "UncondNeighborhood",
L2Fam = "missing"),
function(IC, risk, neighbor, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getBiasIC(IC = IC, neighbor = neighbor,
+
+ if(missing(withCheck)) withCheck <- TRUE
+
+ getBiasIC(IC = IC, neighbor = neighbor,
biastype = biastype(risk), normtype = normtype(risk), tol = tol,
withCheck = withCheck)
})
@@ -82,7 +91,8 @@
neighbor = "UncondNeighborhood",
L2Fam = "L2ParamFamily"),
function(IC, risk, neighbor, L2Fam, tol = .Machine$double.eps^0.25, withCheck = TRUE){
- getBiasIC(IC = IC, neighbor = neighbor, L2Fam = L2Fam,
+ if(missing(withCheck)) withCheck <- TRUE
+ getBiasIC(IC = IC, neighbor = neighbor, L2Fam = L2Fam,
biastype = biastype(risk), normtype = normtype(risk),
tol = tol, withCheck = withCheck)
})
@@ -94,6 +104,7 @@
neighbor = "UncondNeighborhood",
L2Fam = "missing"),
function(IC, risk, neighbor, tol = .Machine$double.eps^0.25, withCheck = TRUE){
+ if(missing(withCheck)) withCheck <- TRUE
L2Fam <- eval(IC at CallL2Fam)
getRiskIC(IC = IC, risk = risk, neighbor = neighbor,
L2Fam = L2Fam, tol = tol, withCheck = withCheck)
@@ -107,6 +118,7 @@
if(dimension(Domain(IC at Curve[[1]])) != dimension(img(L2Fam at distribution)))
stop("dimension of 'Domain' of 'Curve' != dimension of 'img' of 'distribution' of 'L2Fam'")
+ if(missing(withCheck)) withCheck <- TRUE
rad <- neighbor at radius
if(rad == Inf) return(Inf)
More information about the Robast-commits
mailing list