Displaying 1 result from an estimated 1 matches for "qtwo".
Did you mean:
_two
2012 Sep 06
2
choose() function returning anomalous results (zero instead of one)
...ts of debug code).
Shayne Hodge
schodge@ieee.org
# VeridianDynamics Inventory Optimization Script
# Copyright 2012 Shayne Hodge
# Two-vendor Sensitivity Analysis
rm(list=ls()) # Clear previous session data
library("R.basic")
ptm <- proc.time()
ugly_function <- function(nn,k,Qone,Qtwo,f11,f12,f21,f22)
{
partial1 <- (((1-f11)*(1-f12))^(Qone-nn))*((f11+(1-f11)*f12)^nn)
partial2 <- (((1-f21)*(1-f22))^(Qtwo-k))* ((f21+(1-f21)*f22)^(k))
#ifelse ( ((Qone < 125)|((Qone-nn)<75)|(nn<75)), test1 <-
nChooseK(Qone,nn), test1 <- choose(Qone,nn))
#ifelse ( ((Qtwo < 125...