Nothing
clusOpt2 <- function(C1, C2, delta, unit.rv, k=1, CV0=NULL, tot.cost=NULL, cal.sw){
if (!is.null(CV0) & !is.null(tot.cost))
stop("CV0 and tot.cost cannot both be non-null.\n")
if (cal.sw==1 & is.null(tot.cost))
stop("When cal.sw=1, tot.cost cannot be null.\n")
if (cal.sw==2 & is.null(CV0))
stop("When cal.sw=2, CV0 cannot be null.\n")
c.ratio <- C1/C2
n.opt <- sqrt(c.ratio * (1-delta)/delta)
if (cal.sw == 1){
m.opt <- tot.cost / (C1 + C2*n.opt)
if (any(m.opt < 0)) stop(paste("m.opt is negative. Check inputs. m.opt=",m.opt,"\n"))
CV <- sqrt(unit.rv/m.opt/n.opt*k*(1 + delta*(n.opt-1)))
output <-
structure(list(C1 = C1,
C2 = C2,
delta = delta,
"unit relvar" = unit.rv,
k = k,
cost = tot.cost,
m.opt = round(m.opt,1),
n.opt = round(n.opt,1),
CV = round(CV,4)),
class = "power.htest")
}
if (cal.sw == 2) {
m.opt <- unit.rv * k * (1 + delta*(n.opt-1)) / n.opt / CV0^2
if (any(m.opt < 0)) stop(paste("m.opt is negative. Check inputs. m.opt=",m.opt,"\n"))
cost <- C1*m.opt + C2*m.opt*n.opt
output <-
structure(list(C1 = C1,
C2 = C2,
delta = delta,
"unit relvar" = unit.rv,
k = k,
cost = cost,
m.opt = round(m.opt,1),
n.opt = round(n.opt,1),
CV = CV0),
class = "power.htest")
}
output
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.