计算qnorm值的倒数

Question

我想计算一下功率。在一致性研究中，我采用了Walter的一个公式来计算样本。然而，当我有科目数量时，我需要使用它来计算幂。

calcul_f <- function(alpha,puissance,ccc0,ccc1,nb){
cat("ccc0 correpond to the concordance to compare et ccc1 the new concordance et nb the number of replicate")
U_alpha <- abs(qnorm(alpha))
U_beta <- qnorm(puissance)

A <- ccc0/(1-ccc0)
B <- ccc1/(1-ccc1)
C <- (1+(nb*A)) / (1+(nb*B))


NSN <- 1+((2*((U_alpha+U_beta)^2)*nb)/(((log(C))^2)*(nb-1)))
NSN
}
calcul_f (alpha = 0.05,
puissance = 0.8,
ccc0 = 0.6,
ccc1 = 0.75,
nb = 2)

Answer 1

qnorm的逆函数是pnorm：

> pnorm(qnorm(0.2))
[1] 0.2
> qnorm(pnorm(4.2))
[1] 4.2

这适用于每个概率分布，因为分位数(前缀为q的R函数(被定义为累积分布函数的逆函数(前缀为p(。