[Robast-commits] r198 - pkg/ROptEst/inst/scripts

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

Author: stamats
Date: 2008-10-31 12:39:14 +0100 (Fri, 31 Oct 2008)
New Revision: 198

Modified:
   pkg/ROptEst/inst/scripts/ExponentialScaleModel.R
Log:
added computation of confidence intervals, seems to work correctly ...

Modified: pkg/ROptEst/inst/scripts/ExponentialScaleModel.R
===================================================================
--- pkg/ROptEst/inst/scripts/ExponentialScaleModel.R	2008-10-31 10:53:08 UTC (rev 197)
+++ pkg/ROptEst/inst/scripts/ExponentialScaleModel.R	2008-10-31 11:39:14 UTC (rev 198)
@@ -109,7 +109,15 @@
 estimate(E1.est22)
 estimate(E1.est23)
 
+## confidence intervals
+confint(E1.est11, symmetricBias())
+confint(E1.est12, symmetricBias())
+confint(E1.est13, symmetricBias())
+confint(E1.est21, symmetricBias())
+confint(E1.est22, symmetricBias())
+confint(E1.est23, symmetricBias())
 
+
 ## 4. one-step estimation: radius interval
 E1.IC11 <- radiusMinimaxIC(L2Fam=ExpScaleFamily(scale=estimate(E1.est01)),
                 neighbor=ContNeighborhood(), risk=asMSE(), loRad=0, upRad=Inf)
@@ -135,3 +143,11 @@
 estimate(E1.est41)
 estimate(E1.est42)
 estimate(E1.est43)
+
+## confidence intervals
+confint(E1.est31, symmetricBias())
+confint(E1.est32, symmetricBias())
+confint(E1.est33, symmetricBias())
+confint(E1.est41, symmetricBias())
+confint(E1.est42, symmetricBias())
+confint(E1.est43, symmetricBias())