求可自动计算药物筛选质控指标Z'prime值的R包/函数

一、现成第三方包函数 #

• screening包：内置zprime()函数专门用于计算该指标，传入阳性对照和阴性对照的数值向量即可：

library(screening)
zprime(pos = positive_control_data, neg = negative_control_data)

• hitandrun包：提供的z.prime()函数同样支持计算，调用逻辑一致：

library(hitandrun)
z.prime(positive_control_data, negative_control_data)

二、自定义函数（轻量方案） #

如果不想安装额外包，可使用自定义函数，这里在你提供的基础上添加了缺失值处理，适配实验数据常见场景：

Z_prime <- function(pos, neg) {
  pos_mean <- mean(pos, na.rm = TRUE)
  neg_mean <- mean(neg, na.rm = TRUE)
  pos_sd <- sd(pos, na.rm = TRUE)
  neg_sd <- sd(neg, na.rm = TRUE)
  res <- 1 - 3 * (pos_sd + neg_sd) / abs(pos_mean - neg_mean)
  return(res)
}

# 调用示例
# Z_prime(pos = c(10.2, 9.8, 10.5), neg = c(1.1, 0.9, 1.2))

内容的提问来源于stack exchange，提问作者Sophia M.Μ.