用R语言计算统计学习题(下)

上集:用R语言计算统计学习题(上)

第8章 置信区间估计

8.25,电话公司调查安装费用,随机选取了500个家庭用户显示,其中135个愿意增购。

问题:求愿意增购的99%置信区间估计

这道题我搜索了一阵子,最后终于发现prop.test函数,方法就很简单了:

代码:prop.test(135,500,0.99)  135代表样本数X,500代表总数N,0.99代表置信区间。

答案是在0.23-0.31之间。


第9章 假设检验基础:单侧检验

9.53 据估计,约50%的消费者选择商品后又取消了订单。一家电商在新系统下选取了500个样本,其中210个在选择商品后取消了订单。

问题:在0.01的显著性水平下,是否有证据表明这类消费者比例低于50%?

假设样本容量为100,且有42个取消了订单?

因为情况只有取消或者不取消两种情况,参考《统计建模与R软件》P221种子发芽例题,可以用二项分布检验函数来做

代码为binom.test(210,500,p=0.5)

从结果可以看出,在500个样本测试中,p-value<0.01,所以拒绝原假设,比例低于原来的50%;而在100个样本测试中,p-value>0.01,所以同意原假设。

9.68 某银行在12-1点选取了15名顾客作为样本,记录了他们的排队等待时间(数据略)

问题:在0.05显著性水平上,有没有证据表明平均等待时间少于5分钟?

这道题用 单样本的t检验就可以了。代码是:

t.test(time,alternative = "greater", mu = 5)

从运算结果可以看出,P值拒绝原假设。


第10章 两个样本数值数据假设检验和单向方差分析

10.23 两组数据分别采集自接受某种治疗的病人前后情况,是否有证据表明治疗前后有所提高。

这道题当然是还用万能的 t检验。

概率p值大于0.05,无法拒绝原假设H0=H1,即移植前后密度没有显著变高。

这道题也可以用wilcox检验(R语言实战P150, 统计建模与R语言 P245),它是组间差异的非参数检验。

结论也是一样的。

不过诡异的是,10.82的R语言运行结果就和书本答案不同...各位可以自己run一下。


第14章 质量管理中的统计应用

某铁工厂连续32天,各检查500个样品,次品数如列表(数据略)。

问题:建立相应的P图,以及该过程是否处于统计控制状态。

这道题其实也是先用ts函数,画出时间序列图,

R语言里好像没有现成的P值函数,所以我按照《商务统计学》P446的公式自己写了一个。算出上下控制限值后,以此用abline画出参考线。

tproduct<-ts(steel,start=c(2016,1,1),frequency=365)

p=sum(steel)/500/32

UCL=p+3*sqrt((p*(1-p))/500)

LCL=p-3*sqrt((p*(1-p))/500)

plot(tproduct/500,ylim=c(0.018,0.08))

abline(h=mean(steel/500),col="red",lty=6)

abline(h=UCL,col="blue",lty=5)

可以看出,波动本身是在可允许范围内。

R语言本身的时间序列包也有波动值函数,主要有pp检验和adf检验。

install.packages("tseries");library(tseries)

pp.test(tproduct)

adf.test(tproduct)

但素,pp检验和adf检验的结果是相反的,两者的p值相差很远。pp检验和书的结论一致。搜索网页其他人好像也有类似情况,原因如何,还需探究。

《商务统计学》还有很多习题,但个人能力有限,很多未能和R语言函数有效联系。就写到这儿吧。希望各位不吝指教。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 219,928评论 6 509
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 93,748评论 3 396
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 166,282评论 0 357
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 59,065评论 1 295
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 68,101评论 6 395
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 51,855评论 1 308
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 40,521评论 3 420
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 39,414评论 0 276
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 45,931评论 1 319
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 38,053评论 3 340
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 40,191评论 1 352
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 35,873评论 5 347
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 41,529评论 3 331
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 32,074评论 0 23
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 33,188评论 1 272
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 48,491评论 3 375
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 45,173评论 2 357

推荐阅读更多精彩内容