[Distr-commits] r934 - in branches/distr-2.6/pkg/distrMod: R man

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Tue Apr 15 23:28:56 CEST 2014


Author: ruckdeschel
Date: 2014-04-15 23:28:56 +0200 (Tue, 15 Apr 2014)
New Revision: 934

Modified:
   branches/distr-2.6/pkg/distrMod/R/existsPIC.R
   branches/distr-2.6/pkg/distrMod/man/existsPIC.Rd
Log:
distrMod: fixed some minor issue in existsPIC (in case we get 0 matrix, and less strict tolerance)

Modified: branches/distr-2.6/pkg/distrMod/R/existsPIC.R
===================================================================
--- branches/distr-2.6/pkg/distrMod/R/existsPIC.R	2014-04-12 13:07:02 UTC (rev 933)
+++ branches/distr-2.6/pkg/distrMod/R/existsPIC.R	2014-04-15 21:28:56 UTC (rev 934)
@@ -5,7 +5,7 @@
     {kerA.perp <- A.svd$v[,A.svd$d > tol * max(A.svd$d)]
      ## projector to ker A perp
      Pi.kerA.perp <- kerA.perp%*%solve(t(kerA.perp)%*%kerA.perp, tol = tol)%*%t(kerA.perp)
-}else{Pi.kerA.perp <- 0}
+}else{Pi.kerA.perp <- 0*A}
 
 B <- as.matrix(B)
 B.svd <- svd(B, nu = 0)
@@ -13,11 +13,11 @@
     {kerB.perp <- B.svd$v[,B.svd$d > tol * max(B.svd$d)]
      ## projector to ker B perp
      Pi.kerB.perp <- kerB.perp%*%solve(t(kerB.perp)%*%kerB.perp, tol = tol)%*%t(kerB.perp)
-}else{Pi.kerB.perp <- 0}
+}else{Pi.kerB.perp <- 0*B}
 isTRUE(all.equal(Pi.kerB.perp%*%Pi.kerA.perp, Pi.kerB.perp, tolerance = tol ))
 }
 
-setMethod("existsPIC", "L2ParamFamily", function(object, warning = TRUE, tol = .Machine$double.eps){
+setMethod("existsPIC", "L2ParamFamily", function(object, warning = TRUE, tol = .Machine$double.eps^.5){
 if(!isKerAinKerB(object at FisherInfo, trafo(object), tol = tol))
   {if(warning)
       warning("trafo of parameter is not (locally) identifyable for parameter theta.")

Modified: branches/distr-2.6/pkg/distrMod/man/existsPIC.Rd
===================================================================
--- branches/distr-2.6/pkg/distrMod/man/existsPIC.Rd	2014-04-12 13:07:02 UTC (rev 933)
+++ branches/distr-2.6/pkg/distrMod/man/existsPIC.Rd	2014-04-15 21:28:56 UTC (rev 934)
@@ -8,13 +8,13 @@
 (partial) influence curves to Trafo \eqn{D_\theta}{D_theta}.}
 
 \usage{existsPIC(object, ...)
-\S4method{existsPIC}{L2ParamFamily}(object, warning = TRUE, tol = .Machine$double.eps)
+\S4method{existsPIC}{L2ParamFamily}(object, warning = TRUE, tol = .Machine$double.eps^.5)
 }
 \arguments{
   \item{object}{L2ParamFamily}
   \item{...}{further arguments used by specific methods.}
   \item{warning}{logical: should a warning be issued if there exist no (partial) influence curves?}
-  \item{tol}{ the tolerance the linear algebraic operations.  Default is \code{.Machine$double.eps}.}
+  \item{tol}{ the tolerance the linear algebraic operations.  Default is \code{.Machine$double.eps^.5}.}
 }
 \details{To check the existence of (partial) influence curves and, simultaneously, for bounded
 (partial) influence curves, by Lemma 1.1.3 in Kohl(2005) [resp. the fact that \eqn{{\rm ker} I={\rm ker} J}{ker I = ker J}



More information about the Distr-commits mailing list