Displaying 1 result from an estimated 1 matches for "deltau".
Did you mean:
delta
2006 Mar 30
1
warning message in hand-made function
...(Crawford &
Garthwaite, 2002, Neuropsychologia)
c <- (patient-mean.c)/sd.c
#finding the non central parameter of t distributions
f <- function(delta, pr, x, df) pt(x, df = df, ncp = delta) - pr
deltaL <- uniroot(f, lower=-37.62, upper=37.62, pr = 0.025, x = c*
(n^0.5), df = dl)
deltaU <- uniroot(f, lower=-37.62, upper=37.62, pr = 0.975, x = c*
(n^0.5), df = dl)
CI.U <- pnorm(deltaL$root/(n^0.5)) * 100 #upper bound of the
confidence interval
CI.L <- pnorm(deltaU$root/(n^0.5)) * 100 #lower bound of the
confidence interval
#output
output <- list(statistic=t.o...