95%置信区间

amlfly

求助各位大佬，阳性符合率、阴性符合率、总符合率的95%置信区间spss计算方法。比如后面的置信区间如何获取的，求大佬解答，不胜感激。

	对比试剂				符合率+置信区间
考核试剂		阳性	阴性	合计
	阳性	50	5	55	90.00%;（83.31%，98.51%）
	阴性	5	276	281	98.22%;（96.67%，99.77%）
	合计	55	281	336	97.02%;（95.21%，98.84%）

光脚大象

区间估计，可以翻一下书《参数估计》章节，很详细

瓶瓶

能用Excel实现的不用spss

负红颜

library(readxl)
df <- read_excel('考核试剂.xlsx',sheet= 2)

library(tidyverse)
table <- df %>% select(对比阳,对比阴) %>% as.matrix()

# 计算阳性符合率
positive_rate <- table[1, 1] / sum(table[1, ])

# 计算阴性符合率
negative_rate <- table[2, 2] / sum(table[2, ])

# 计算总符合率
total_rate <- sum(diag(table)) / sum(table)

# 计算阳性符合率的置信区间
positive_ci <- binom.test(table[1, 1], sum(table[1, ]))$conf.int

# 计算阴性符合率的置信区间
negative_ci <- binom.test(table[2, 2], sum(table[2, ]))$conf.int

# 计算总符合率的置信区间
total_ci <- binom.test(sum(diag(table)), sum(table))$conf.int

# 打印结果
print(paste("阳性符合率:", positive_rate))
print(paste("阳性符合率的95%置信区间:", positive_ci))
print(paste("阴性符合率:", negative_rate))
print(paste("阴性符合率的95%置信区间:", negative_ci))
print(paste("总符合率:", total_rate))
print(paste("总符合率的95%置信区间:", total_ci))

-----------------------------------------------以上是程序---------------下面是结果

> # 打印结果
> print(paste("阳性符合率:", positive_rate))[1] "阳性符合率: 0.909090909090909"
> print(paste("阳性符合率的95%置信区间:", positive_ci))[1] "阳性符合率的95%置信区间: 0.800460384229499" "阳性符合率的95%置信区间: 0.969819132154187"
> print(paste("阴性符合率:", negative_rate))[1] "阴性符合率: 0.98220640569395"
> print(paste("阴性符合率的95%置信区间:", negative_ci))[1] "阴性符合率的95%置信区间: 0.958966473902338" "阴性符合率的95%置信区间: 0.994197855216952"
> print(paste("总符合率:", total_rate))[1] "总符合率: 0.970238095238095"
> print(paste("总符合率的95%置信区间:", total_ci))[1] "总符合率的95%置信区间: 0.94594920046414"  "总符合率的95%置信区间: 0.985637629761318"

负红颜

可见，你数都算错了。

amlfly

负红颜 发表于 2023-4-24 10:05
library(readxl)
df % as.matrix()

这是用公式算的可能有点误差。好厉害的样子，我看不懂。参考区间那部分的操作步骤大佬能帮忙截个图或者简单用语言描述一下吗

负红颜

amlfly 发表于 2023-4-24 10:19
这是用公式算的可能有点误差。好厉害的样子，我看不懂。参考区间那部分的操作步骤大佬能帮忙截个图或者简 ...

查二项分布的置信区间

amlfly

负红颜 发表于 2023-4-24 10:26
查二项分布的置信区间

好的，万分感谢

wzzz2008

负红颜 发表于 2023-4-24 10:05
library(readxl)
df % as.matrix()

对于我看不懂的东西，只能说一句：好厉害

xinhaons

你这个结果应该是用二项分布近似正态分布估算的，用excel就可以计算。当然用SPSS或Minitab会更快而且可以使用精确计算。
近似的算法以第一行举例:
95%的正态分布k=1.96
标准差的估计为s=(p(1-p)/n)^0.5  即s=0.909091*(1-0.909091)/55再开根号
则置信区间为0.909091±1.96*s

判正确	判错误	p	s	ks	LL	UL
50	5	0.909091	0.038763767	0.075977	0.833114	0.985068
276	5	0.982206	0.007886423	0.015457	0.966749	0.997664
326	10	0.970238	0.00927043	0.01817	0.952068	0.988408

SPSS我没用过，用minitab精确计算可以用单比率的假设检验工具立即可以算出（当然也可以选择正态近似的方法）

glm1024

负红颜 发表于 2023-4-24 10:05
library(readxl)
df % as.matrix()

大佬用R计算就是牛逼。

负红颜

glm1024 发表于 2023-4-24 21:17
大佬用R计算就是牛逼。

别，您是大佬，我是跟着您进入R的。

胖橘12

有个计算的软件，发我邮箱，我发你