时间戳

首页 > 解决方案 > 计算 qnorm 值的倒数

r - 计算 qnorm 值的倒数

问题描述

我想计算一下功率。我采用了 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)

标签： r

解决方案

qnorm的反函数是pnorm：

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

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

推荐阅读