[R语言]总体概率的区间估计学习

总体概率的区间估计没有总体均值的区间估计常用，但是用途也不小，在科研论文中，经常会使用该区间报告p值的意义，所以将做法记录如下：

First,remember that an interval for a proportion is given by:

p_hat +/- z * sqrt(p_hat * (1-p_hat)/n)

With that being said, we can use R to solvethe formula like so:

# Load Data
vehicleType = c("suv", "suv", "minivan", "car", "suv", "suv", "car", "car", "car", "car", "minivan", "car", "truck", "car", "car", "car", "car", "car", "car", "car", "minivan", "car", "suv", "minivan", "car", "minivan", "suv", "suv", "suv", "car", "suv", "car", "car", "suv", "truck", "truck", "minivan", "suv", "car", "truck", "suv", "suv", "car", "car", "car", "car", "suv", "car", "car", "car", "suv", "car", "car", "car", "truck", "car", "car", "suv", "suv", "minivan", "suv", "car", "car", "car", "car", "car", "minivan", "suv", "car", "car", "suv", "minivan", "car", "car", "car", "minivan", "minivan", "minivan", "car", "truck", "car", "car", "car", "suv", "suv", "suv", "car", "suv", "suv", "car", "suv", "car", "minivan", "car", "car", "car", "car", "car", "car", "car")

# Convert from string to factor
vehicleType = factor(vehicleType)

# Find the number of obs
n = length(vehicleType)

# Find number of obs per type
vtbreakdown = table(vehicleType)

# Get the proportion
p_hat = vtbreakdown['suv']/n#到这里就是仅仅是为了求出p值。

 

#下面的α指定完成以后，就可以用p值进行求置信区间了。加减变换一下就行。

# Set CI alpha level (1-alpha/2)*100%
alpha = 0.05

# Calculate the critical z-score
z = qnorm(1-alpha/2)

# Compute the CI
p_hat + c(-1,1)*z*sqrt(p_hat*(1-p_hat)/n)#CI上限

p_hat-c(-1,1)*z*sqrt(p_hat*(1-p_hat)/n)#CI下限

如果我的博客帮助到你或者你喜欢我的博客，请给我点个赞，你的支持是博主继续写下去动力。而如果你对我的博客有什么问题或建议，那可以给博主发邮件，博主会选择有代表性的给大家解答。

。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。

你好！我是幸福司马缸，一个幸福的思想家。

幸福司马缸——专注于探索幸福和高品质的生活，统计，健康，教育，投资。

目前在学习公共卫生，英语，统计和大数据，

业余时间研究如何幸福的生活，及如何让思考和知识变得有趣：包括人生的方方面面，家庭，生活，教育，投资，健康。

如果你想要联系到我，可以给我发邮件：[e-mail]

如果你想跟我一起学习，你也可以加入我的qq群：155454404，期待你和我一起讨论，一起思考。

践行有态度的思考与创作。

本文欢迎转发，转发请注明作者和出处。谢谢！

 

很高兴你能够看到我的博客，欢迎到我的其他博客去逛逛，希望能够给你更多的收获

http://blog.csdn.net/qq_36852770：Ash_Zhang的技术小屋；记录IT及代码的学习，以及极客生活

http://blog.sina.com.cn/happinessthinker：幸福司马缸；记录个人管理生活感悟