[Distr-commits] r179 - branches/distr-2.0/pkg/distr/R

noreply at r-forge.r-project.org noreply at r-forge.r-project.org
Mon Jun 16 11:21:37 CEST 2008


Author: stamats
Date: 2008-06-16 11:21:37 +0200 (Mon, 16 Jun 2008)
New Revision: 179

Modified:
   branches/distr-2.0/pkg/distr/R/internalUtils.R
Log:
there is a bug in convpow respectively .makePnew
Consider
D1 <- Norm()
D2 <- Exp()
D3 <- D1 + D2
convpow(D3, 5)

I already modified x[l2+1,l] to x[l2+1:l]. Is this modification correct? Peter?

Modified: branches/distr-2.0/pkg/distr/R/internalUtils.R
===================================================================
--- branches/distr-2.0/pkg/distr/R/internalUtils.R	2008-06-16 09:17:47 UTC (rev 178)
+++ branches/distr-2.0/pkg/distr/R/internalUtils.R	2008-06-16 09:21:37 UTC (rev 179)
@@ -673,9 +673,9 @@
          if ((l%%2==0)&& is.null(myPf)){
                l2 <- l/2
                if (is.null(pxl))
-                   x.l <- c(x[1:l2],(x[l2]+x[l2+1])/2,x[l2+1,l])
+                   x.l <- c(x[1:l2],(x[l2]+x[l2+1])/2,x[l2+1:l])
                if (is.null(pxu))
-                   x.u <- c(x[1:l2],(x[l2]+x[l2+1])/2,x[l2+1,l])
+                   x.u <- c(x[1:l2],(x[l2]+x[l2+1])/2,x[l2+1:l])
                l <- l+1
                }
          cfun <- .csimpsum



More information about the Distr-commits mailing list