[Robast-commits] r408 - branches/robast-0.8/pkg/RobAStBase/R

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

Author: ruckdeschel
Date: 2010-06-11 00:04:52 +0200 (Fri, 11 Jun 2010)
New Revision: 408

Modified:
   branches/robast-0.8/pkg/RobAStBase/R/comparePlot.R
   branches/robast-0.8/pkg/RobAStBase/R/infoPlot.R
Log:
fixed a little bug in comparePlot if we want to plot pIC's

Modified: branches/robast-0.8/pkg/RobAStBase/R/comparePlot.R
===================================================================
--- branches/robast-0.8/pkg/RobAStBase/R/comparePlot.R	2010-05-28 20:35:02 UTC (rev 407)
+++ branches/robast-0.8/pkg/RobAStBase/R/comparePlot.R	2010-06-10 22:04:52 UTC (rev 408)
@@ -49,7 +49,7 @@
 
         trafO <- trafo(L2Fam at param)
         dims  <- nrow(trafO)
-        dimm <- length(L2Fam at param)
+        dimm <- ncol(trafO)
         
         to.draw <- 1:dims
         dimnms  <- c(rownames(trafO))
@@ -110,8 +110,8 @@
         dots$xlim <- NULL
 
         dims <- nrow(trafo(L2Fam at param))
-        IC1 <- as(diag(dimm) %*% obj1 at Curve, "EuclRandVariable")
-        IC2 <- as(diag(dimm) %*% obj2 at Curve, "EuclRandVariable")
+        IC1 <- as(diag(dims) %*% obj1 at Curve, "EuclRandVariable")
+        IC2 <- as(diag(dims) %*% obj2 at Curve, "EuclRandVariable")
 
 
         obj <- obj3
@@ -119,7 +119,7 @@
            {
            if(!identical(L2Fam1c,obj at CallL2Fam))
                stop("ICs need to be defined for the same model")
-           IC3 <- as(diag(dimm) %*% obj3 at Curve, "EuclRandVariable")
+           IC3 <- as(diag(dims) %*% obj3 at Curve, "EuclRandVariable")
            }
 
         obj <- obj4
@@ -127,7 +127,7 @@
            {
            if(!identical(L2Fam1c,obj at CallL2Fam))
                stop("ICs need to be defined for the same model")
-           IC4 <- as(diag(dimm) %*% obj4 at Curve, "EuclRandVariable")
+           IC4 <- as(diag(dims) %*% obj4 at Curve, "EuclRandVariable")
            }
 
       lineT <- NA

Modified: branches/robast-0.8/pkg/RobAStBase/R/infoPlot.R
===================================================================
--- branches/robast-0.8/pkg/RobAStBase/R/infoPlot.R	2010-05-28 20:35:02 UTC (rev 407)
+++ branches/robast-0.8/pkg/RobAStBase/R/infoPlot.R	2010-06-10 22:04:52 UTC (rev 408)
@@ -30,7 +30,7 @@
         
         trafO <- trafo(L2Fam at param)
         dims <- nrow(trafO)
-        dimm <- length(L2Fam at param)
+        dimm <- ncol(trafO)
         
         to.draw <- 1:(dims+1)
         dimnms <- rownames(trafO)