search for: d1d4pos

...quot;numeric")     for (i in 1:1000)   { cdf  <- integrate(fprime, lower = 0, upper = z)$value     f    <- cdf - r     # Newton method     z    <- z - f/fprime         if (any(f < tol)) break    }   cbind(z,cdf) }   alp  <- 2.0165 bt1  <- 29.107 ; bt2 <- 41.517 x1   <- d1d4pos[,1];x1[1:10] x2   <- d1d4pos[,2];x2[1:10]   > x1   <- d1d4pos[,1];x1[1:10]  [1] 28.4 53.6  1.3 29.5 52.1 65.9 72.6 67.6 58.7 34.5 > x2   <- d1d4pos[,2];x2[1:10]  [1]  43.5  56.2   0.3  16.6  71.1  86.3 172.8 111.8  89.9  70.2   z <- (x1/bt1)+(x2/bt2); z   newton_gam(min(z)) newton...