Displaying 1 result from an estimated 1 matches for "mult0".
Did you mean:
multi
2016 Apr 02
0
BCa Bootstrap confidence intervals
...s a confidence interval for a semipartial correlation
# Arguments:
# alpha: alpha value for 1-alpha confidence
# part: sample semipartial correlation
# mult: sample squared multiple correlation in full model
# n: sample size
# Returns:
# confidence interval
z <- qnorm(1 - alpha/2)
mult0 <- mult - part^2
zr <- log((1 + part)/(1 - part))/2
a <- (mult^2 - 2*mult + mult0 - mult0^2 + 1)/(1 - part^2)^2
se <- sqrt(a/(n - 3))
LL0 <- zr - z*se
UL0 <- zr + z*se
LL <- (exp(2*LL0) - 1)/(exp(2*LL0) + 1)
UL <- (exp(2*UL0) - 1)/(exp(2*UL0) + 1)
CI <- c(LL, UL)
return(C...