[Robast-commits] r1122 - in branches/robast-1.2/pkg/ROptEst: R inst
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Sat Aug 11 16:34:07 CEST 2018
Author: ruckdeschel
Date: 2018-08-11 16:34:06 +0200 (Sat, 11 Aug 2018)
New Revision: 1122
Modified:
branches/robast-1.2/pkg/ROptEst/R/LowerCaseMultivariate.R
branches/robast-1.2/pkg/ROptEst/R/getComp.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/getInfRobIC_asBias.R
branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asCov.R
branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asGRisk.R
branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asHampel.R
branches/robast-1.2/pkg/ROptEst/R/getInfStand.R
branches/robast-1.2/pkg/ROptEst/R/getInfV.R
branches/robast-1.2/pkg/ROptEst/R/getMaxIneff.R
branches/robast-1.2/pkg/ROptEst/R/getRadius.R
branches/robast-1.2/pkg/ROptEst/R/getStartIClcsc.R
branches/robast-1.2/pkg/ROptEst/R/leastFavorableRadius.R
branches/robast-1.2/pkg/ROptEst/R/optRisk.R
branches/robast-1.2/pkg/ROptEst/R/radiusMinimaxIC.R
branches/robast-1.2/pkg/ROptEst/R/roptest.new.R
branches/robast-1.2/pkg/ROptEst/R/updateNorm.R
branches/robast-1.2/pkg/ROptEst/inst/NEWS
Log:
[ROptEst] branch 1.2
+ now specified that we want to use distr::solve
+ internal function .getComp, determining by symmetry slots which entries in LMs a and A
have to be computed, now fills the lower triangle of A with FALSE (was not used so far,
but can be used in a faster computation method for checkIC makeIC to determine whether
it is cleverer to integrate in k or in p space)
+ begun with speedup code for checkIC makeIC
Modified: branches/robast-1.2/pkg/ROptEst/R/LowerCaseMultivariate.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/LowerCaseMultivariate.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/LowerCaseMultivariate.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -9,7 +9,7 @@
w <- new("HampelWeight")
if(is.null(z.start)) z.start <- numeric(ncol(trafo))
- if(is.null(A.start)) A.start <- trafo%*%solve(as.matrix(Finfo))
+ if(is.null(A.start)) A.start <- trafo%*%distr::solve(as.matrix(Finfo))
if(is.null(A.comp))
A.comp <- matrix(TRUE, nrow = nrow(trafo), ncol = ncol(trafo))
if(is.null(z.comp))
@@ -109,7 +109,7 @@
w <- new("BdStWeight")
k <- ncol(trafo)
- if(is.null(A.start)) A.start <- trafo%*%solve(Finfo)
+ if(is.null(A.start)) A.start <- trafo%*%distr::solve(Finfo)
pos.fct <- function(x, L2, stand){
X <- evalRandVar(L2, as.matrix(x))[,,1]
Modified: branches/robast-1.2/pkg/ROptEst/R/getComp.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getComp.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getComp.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -18,7 +18,7 @@
A.comp[i,j] <- FALSE
}
}
- A.comp[col(A.comp) < row(A.comp)] <- A.comp[col(A.comp) > row(A.comp)]
+ A.comp[col(A.comp) < row(A.comp)] <- FALSE
return(list(A.comp = A.comp, z.comp = z.comp))
}
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfLM.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfLM.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfLM.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -47,7 +47,7 @@
# print(c("z"=z))
if(is(neighbor,"TotalVarNeighborhood")){
a <- z
- z <- as.numeric(solve(A,a))
+ z <- as.numeric(distr::solve(A,a))
zc <- numeric(ncol(trafo))
}else if(is(neighbor,"ContNeighborhood")) {
zc <- z
@@ -159,7 +159,7 @@
# print(list(A0vecA,A0,a0))
- z0 <- as.numeric(solve(A0,a0))
+ z0 <- as.numeric(distr::solve(A0,a0))
std0 <- stdC
w0 <- w1
risk0 <- risk1
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asAnscombe.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asAnscombe.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -132,9 +132,9 @@
stop("Not yet implemented")
## non-standard norms
- FI1 <- trafo%*%solve(Finfo)
+ FI1 <- trafo%*%distr::solve(Finfo)
FI0 <- FI1%*%t(trafo)
- FI <- solve(FI0)
+ FI <- distr::solve(FI0)
if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") ){
QuadForm(normtype) <- PosSemDefSymmMatrix(FI)
normtype(risk) <- normtype
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asBias.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asBias.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asBias.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -82,7 +82,7 @@
}
# }
- FI <- solve(trafo%*%solve(Finfo)%*%t(trafo))
+ FI <- distr::solve(trafo%*%distr::solve(Finfo)%*%t(trafo))
if(is(normtype,"QFNorm"))
{QuadForm(normtype) <- PosSemDefSymmMatrix(FI);
normtype(risk) <- normtype}
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asCov.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asCov.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asCov.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -10,7 +10,7 @@
verbose <- getRobAStBaseOption("all.verbose")
info <- c("optimal IC in sense of Cramer-Rao bound")
- A <- trafo %*% solve(Finfo)
+ A <- trafo %*% distr::solve(Finfo)
b <- abs(as.vector(A))*max(abs(q.l(L2deriv)(1)),abs(q.l(L2deriv)(0)))
@@ -43,7 +43,7 @@
verbose <- getRobAStBaseOption("all.verbose")
info <- c("optimal IC in sense of Cramer-Rao bound")
- A <- trafo %*% solve(Finfo)
+ A <- trafo %*% distr::solve(Finfo)
b <- abs(as.vector(A))*(q.l(L2deriv)(1)-q.l(L2deriv)(0))
a <- -abs(as.vector(A))*q.l(L2deriv)(0)
asCov <- A %*% t(trafo)
@@ -80,7 +80,7 @@
if(! Cont && p>1)
stop("Not yet implemented")
info <- c("optimal IC in sense of Cramer-Rao bound")
- A <- trafo %*% solve(Finfo)
+ A <- trafo %*% distr::solve(Finfo)
IC <- A %*% L2deriv
if(is(Distr, "UnivariateDistribution")){
lower <- ifelse(is.finite(q.l(Distr)(0)), q.l(Distr)(1e-8), q.l(Distr)(0))
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asGRisk.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asGRisk.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asGRisk.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -21,7 +21,7 @@
k <- ncol(trafo)
## non-standard norms
- FI <- solve(trafo%*%matrix(1/Finfo,1,1)%*%t(trafo))
+ FI <- distr::solve(trafo%*%matrix(1/Finfo,1,1)%*%t(trafo))
if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") ){
QuadForm(normtype) <- PosSemDefSymmMatrix(FI)
normtype(risk) <- normtype
@@ -221,7 +221,7 @@
stop("Not yet implemented")
## non-standard norms
- FI <- solve(trafo%*%solve(Finfo)%*%t(trafo))
+ FI <- distr::solve(trafo%*%distr::solve(Finfo)%*%t(trafo))
if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") ){
QuadForm(normtype) <- PosSemDefSymmMatrix(FI)
normtype(risk) <- normtype
@@ -232,7 +232,7 @@
## starting values
if(is.null(z.start)) z.start <- numeric(k)
- if(is.null(A.start)) A.start <- trafo %*% solve(Finfo)
+ if(is.null(A.start)) A.start <- trafo %*% distr::solve(Finfo)
a.start <- as.numeric(A.start %*% z.start)
## sort out upper solution if radius = 0
@@ -288,7 +288,7 @@
verbose = verbose)
lower <- lowBerg$b}
#if(is.null(upper))
- upper <- 5*max(solve(Finfo))
+ upper <- 5*max(distr::solve(Finfo))
OptIterCall <- numeric(1)
Cov <- 0
@@ -510,7 +510,7 @@
b <- res$b
res <- c(res, list(biastype = biastype, normtype = normtype))
if(!is(risk, "asMSE")){
- FI <- trafo%*%solve(Finfo)%*%t(trafo)
+ FI <- trafo%*%distr::solve(Finfo)%*%t(trafo)
FI <- sum(diag(QuadForm %*% FI))
Risk <- getAsRisk(risk = risk, L2deriv = L2deriv, neighbor = neighbor,
biastype = biastype, normtype = normtype,
@@ -549,7 +549,7 @@
verbose = verbose)
normtype(risk) <- res$normtype
if(!is(risk, "asMSE")){
- FI <- trafo%*%solve(Finfo)%*%t(trafo)
+ FI <- trafo%*%distr::solve(Finfo)%*%t(trafo)
FI <- sum(diag(QuadForm %*% FI))
Risk <- getAsRisk(risk = risk, L2deriv = L2deriv, neighbor = neighbor,
biastype = biastype(risk), normtype = normtype(risk),
@@ -665,14 +665,14 @@
# if( is(neighbor,"TotalVarNeighborhood") && p>1)
# stop("Not yet implemented")
#
-# FI <- solve(trafo%*%solve(Finfo)%*%t(trafo))
+# FI <- distr::solve(trafo%*%distr::solve(Finfo)%*%t(trafo))
# if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") )
# {QuadForm(normtype) <- PosSemDefSymmMatrix(FI);
# normtype(risk) <- normtype}
# QF <- if(is(normtype,"QFNorm")) QuadForm(normtype) else diag(nrow(trafo))
#
# if(is.null(z.start)) z.start <- numeric(ncol(trafo))
-# if(is.null(A.start)) A.start <- trafo %*% solve(Finfo)
+# if(is.null(A.start)) A.start <- trafo %*% distr::solve(Finfo)
#
# radius <- neighbor at radius
# if(identical(all.equal(radius, 0), TRUE)){
@@ -791,7 +791,7 @@
#
# if(is(neighbor,"TotalVarNeighborhood")){
# a <- z
-# z <- solve(A,a)
+# z <- distr::solve(A,a)
# zc <- numeric(ncol(trafo))
# }else if(is(neighbor,"ContNeighborhood")) {
# zc <- z
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asHampel.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asHampel.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfRobIC_asHampel.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -170,7 +170,7 @@
stop("Not yet implemented")
## non-standard norms
- FI <- solve(trafo%*%solve(Finfo)%*%t(trafo))
+ FI <- distr::solve(trafo%*%distr::solve(Finfo)%*%t(trafo))
if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") ){
QuadForm(normtype) <- PosSemDefSymmMatrix(FI)
normtype(risk) <- normtype
@@ -180,7 +180,7 @@
## starting values
if(is.null(z.start)) z.start <- numeric(k)
- if(is.null(A.start)) A.start <- trafo%*%solve(Finfo)
+ if(is.null(A.start)) A.start <- trafo%*%distr::solve(Finfo)
a.start <- as.numeric(A.start %*% z.start)
## initialize
@@ -334,7 +334,7 @@
tol, QuadForm, verbose, nrvalpts, warn){
if(missing(warn)|| is.null(warn)) warn <- FALSE
- ClassIC <- trafo %*% solve(Finfo) %*% L2deriv
+ ClassIC <- trafo %*% distr::solve(Finfo) %*% L2deriv
lower.x <- getLow(Distr)
upper.x <- getUp(Distr)
@@ -405,7 +405,7 @@
# if(! is(neighbor,"ContNeighborhood") && p>1)
# stop("Not yet implemented")
#
-# FI <- solve(trafo%*%solve(Finfo)%*%t(trafo))
+# FI <- distr::solve(trafo%*%distr::solve(Finfo)%*%t(trafo))
# if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") )
# {QuadForm(normtype) <- PosSemDefSymmMatrix(FI); normtype(risk) <- normtype}
#
@@ -416,7 +416,7 @@
# b <- risk at bound
#
# if(checkBounds){
-# ClassIC <- trafo %*% solve(Finfo) %*% L2deriv
+# ClassIC <- trafo %*% distr::solve(Finfo) %*% L2deriv
# lower.x <- getLow(Distr)
# upper.x <- getUp(Distr)
# x <- seq(from = lower.x, to = upper.x, length = 5000)
@@ -502,7 +502,7 @@
#
# if(is(neighbor,"TotalVarNeighborhood")){
# a <- z
-# z <- solve(A,a)
+# z <- distr::solve(A,a)
# zc <- numeric(ncol(trafo))
# }else if(is(neighbor,"ContNeighborhood")) {
# zc <- z
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfStand.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfStand.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfStand.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -42,7 +42,7 @@
erg[col(erg) < row(erg)] <- erg[col(erg) > row(erg)]
- return(trafo %*% solve(erg))
+ return(trafo %*% distr::solve(erg))
})
###############################################################################
## standardizing constant for one-sided bias
Modified: branches/robast-1.2/pkg/ROptEst/R/getInfV.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getInfV.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getInfV.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -36,12 +36,12 @@
(weight(w)(evalRandVar(L2deriv, as.matrix(x)) [,,1]))^2
}
- .solve <- function(A0, b0) solve(A0,b0)
- if(is.matrix(stand)){
- if(nrow(stand)!=ncol(stand))
- .solve <- function(A0,b0) MASS::ginv(A0)%*%b0
- }
- cent0 <- .solve(stand, cent)
+# .solve <- function(A0, b0) distr::solve(A0,b0)
+# if(is.matrix(stand)){
+# if(nrow(stand)!=ncol(stand))
+# .solve <- function(A0,b0) MASS::ginv(A0)%*%b0
+# }
+ cent0 <- distr::solve(stand, cent, generalized = TRUE)
integrandV <- function(x, L2.i, L2.j, i, j){
Modified: branches/robast-1.2/pkg/ROptEst/R/getMaxIneff.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getMaxIneff.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getMaxIneff.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -21,7 +21,7 @@
Finfo <- L2Fam at FisherInfo
L2derivDim <- numberOfMaps(L2Fam at L2deriv)
- FI0 <- trafo%*%solve(Finfo)%*%t(trafo)
+ FI0 <- trafo%*%distr::solve(Finfo)%*%t(trafo)
std <- if(is(normtype,"QFNorm"))
QuadForm(normtype) else diag(nrow(trafo))
s0 <- sum(diag(std%*%FI0))
Modified: branches/robast-1.2/pkg/ROptEst/R/getRadius.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getRadius.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getRadius.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -23,7 +23,7 @@
}else{
L2deriv <- diag(dimension(L2Fam at L2deriv)) %*% L2Fam at L2deriv
}
- z <- solve(stand(IC),cent(IC))
+ z <- distr::solve(stand(IC),cent(IC))
r <- getInfRad(clip = clip(IC), L2deriv = L2deriv,
risk = risk, neighbor = neighbor, biastype = biastype(risk),
Distr = L2Fam at distribution, stand = stand(IC),
Modified: branches/robast-1.2/pkg/ROptEst/R/getStartIClcsc.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/getStartIClcsc.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/getStartIClcsc.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -101,7 +101,7 @@
LM0$Aw <- LM0$A <- (Aa+Ai+t(Ai)+t(Aa))/4
ai <- Ai %*% zi
LM0$a <- (ai+aa)/2
- LM0$aw <- solve(LM0$A, LM0$a)
+ LM0$aw <- distr::solve(LM0$A, LM0$a)
}
.xiMkLM <- function(LMset,xi){
Modified: branches/robast-1.2/pkg/ROptEst/R/leastFavorableRadius.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/leastFavorableRadius.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/leastFavorableRadius.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -20,8 +20,8 @@
normtype <- normtype(risk)
trafo <- trafo(L2Fam at param)
- FI0 <- trafo%*%solve(L2Fam at FisherInfo)%*%t(trafo)
- FI <- solve(FI0)
+ FI0 <- trafo%*%distr::solve(L2Fam at FisherInfo)%*%t(trafo)
+ FI <- distr::solve(FI0)
if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") )
{QuadForm(normtype) <- PosSemDefSymmMatrix(FI);
normtype(risk) <- normtype}
Modified: branches/robast-1.2/pkg/ROptEst/R/optRisk.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/optRisk.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/optRisk.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -3,7 +3,7 @@
###############################################################################
setMethod("optRisk", signature(model = "L2ParamFamily", risk = "asCov"),
function(model, risk){
- return(list(asCov = solve(model at FisherInfo)))
+ return(list(asCov = distr::solve(model at FisherInfo)))
})
###############################################################################
Modified: branches/robast-1.2/pkg/ROptEst/R/radiusMinimaxIC.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/radiusMinimaxIC.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/radiusMinimaxIC.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -123,10 +123,10 @@
Finfo <- L2Fam at FisherInfo
p <- nrow(trafo)
- FI0 <- trafo%*%solve(Finfo)%*%t(trafo)
+ FI0 <- trafo%*%distr::solve(Finfo)%*%t(trafo)
if(is(normtype,"InfoNorm") || is(normtype,"SelfNorm") )
- {QuadForm(normtype) <- PosSemDefSymmMatrix(solve(FI0));
+ {QuadForm(normtype) <- PosSemDefSymmMatrix(distr::solve(FI0));
normtype(risk) <- normtype}
std <- if(is(normtype,"QFNorm")) QuadForm(normtype) else diag(p)
loRisk <- sum(diag(std%*%FI0))
Modified: branches/robast-1.2/pkg/ROptEst/R/roptest.new.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/roptest.new.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/roptest.new.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -318,7 +318,7 @@
es.list0$fsCor <- eval(es.list0$fsCor)
es.list0$OptOrIter <- eval(es.list0$OptOrIter)
- if(debug) {cat("\n\n\n::::\n\n")
+ if(debug) {cat("\n\n\n:::: args for getStartIC\n\n")
argList <- c(list(model=L2Fam,risk=risk,neighbor=neighbor,
withEvalAsVar = withEvalAsVarSIC, withMakeIC = withMakeICSIC,
modifyICwarn = modifyICwarnSIC), es.list0)
@@ -337,7 +337,7 @@
if(debug){
- ICstart <- "BUL"
+ ICstart <- "ICstart-result-debug"
argList <- list(x, IC = ICstart, start = initial.est, steps = steps,
useLast = kStepCtrl$useLast,
withUpdateInKer = kStepCtrl$withUpdateInKer,
Modified: branches/robast-1.2/pkg/ROptEst/R/updateNorm.R
===================================================================
--- branches/robast-1.2/pkg/ROptEst/R/updateNorm.R 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/R/updateNorm.R 2018-08-11 14:34:06 UTC (rev 1122)
@@ -6,7 +6,7 @@
{Cv <- getInfV(L2deriv = L2, neighbor = neighbor,
biastype = biastype, Distr = Distr,
V.comp = V.comp, cent = cent, stand = stand, w = w)
- QuadForm(normtype) <- PosSemDefSymmMatrix(solve(Cv))
+ QuadForm(normtype) <- PosSemDefSymmMatrix(distr::solve(Cv))
normtype})
Modified: branches/robast-1.2/pkg/ROptEst/inst/NEWS
===================================================================
--- branches/robast-1.2/pkg/ROptEst/inst/NEWS 2018-08-11 13:54:52 UTC (rev 1121)
+++ branches/robast-1.2/pkg/ROptEst/inst/NEWS 2018-08-11 14:34:06 UTC (rev 1122)
@@ -29,6 +29,11 @@
+ inserted code for time checking (which is inactive usually; only if in kStepEstimator.R in
RobAStBase, the respective ##-t-## lines are de-commented the timings are visible as
attribute "kStepTimings" in the result of roptest ...)
++ now specified that we want to use distr::solve
++ internal function .getComp, determining by symmetry slots which entries in LMs a and A
+ have to be computed, now fills the lower triangle of A with FALSE (was not used so far,
+ but can be used in a faster computation method for checkIC makeIC to determine whether
+ it is cleverer to integrate in k or in p space)
#######################################
version 1.1
More information about the Robast-commits
mailing list