[Gsdesign-commits] r374 - in pkg/gsDesign: . R

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

Author: keaven
Date: 2014-08-10 21:27:14 +0200 (Sun, 10 Aug 2014)
New Revision: 374

Modified:
   pkg/gsDesign/DESCRIPTION
   pkg/gsDesign/R/ssrCP.R
Log:
gsCP update and documentation update

Modified: pkg/gsDesign/DESCRIPTION
===================================================================
--- pkg/gsDesign/DESCRIPTION	2014-08-10 19:26:11 UTC (rev 373)
+++ pkg/gsDesign/DESCRIPTION	2014-08-10 19:27:14 UTC (rev 374)
@@ -1,5 +1,5 @@
 Package: gsDesign
-Version: 2.9-1
+Version: 2.9-2
 Title: Group Sequential Design
 Author: Keaven Anderson 
 Maintainer: Keaven Anderson <keaven_anderson at merck.com>

Modified: pkg/gsDesign/R/ssrCP.R
===================================================================
--- pkg/gsDesign/R/ssrCP.R	2014-08-10 19:26:11 UTC (rev 373)
+++ pkg/gsDesign/R/ssrCP.R	2014-08-10 19:27:14 UTC (rev 374)
@@ -244,12 +244,13 @@
                             target=1-x$beta, x=x$x, 
                             n2=x$maxinc*x$x$n.I[2]-x$x$n.I[1],
                             z2=x$z2fn, theta=x$theta)$root
+  
     if (changepoint3 >= changepoint2){
       en <- en + x$maxinc*x$x$n.I[2]*(pnorm(changepoint2-mu)-
                                         pnorm(changepoint-mu))
       a <- changepoint
       b <- changepoint2
-      n2 <- x$x$n.I[2]-x$x$n.I[1]
+      n2 <- x$maxinc*x$x$n.I[2]-x$x$n.I[1]
       grid <- normalGrid(mu=(a+b)/2,bounds=c(a,b),r=r)
       for(i in 1:length(theta)) Power[i] <-  Power[i] +
         sum(dnorm(grid$z-sqrt(x$x$n.I[1])*theta[i]) *