原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第七章
1．
a．
lamp-data.frame(
+ X=c(115,116,98,83,103,107,118,116,73,89,85,97),
+ A=factor(rep(1:3,c(4,4,4)))
+ )
lamp.aov-aov(X ~ A, data=lamp)
source(anova.tab.R);anova.tab(lamp.aov)
Df Sum Sq Mean Sq F value Pr(F)
A 2 1304.00 652.00 4.9228 0.03595 *
Residuals 9 1192.00 132.44
Total 11 2496.00
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘.’ 1
P值大于0.01，接受原假设。没有显著差异
b．
source(interval_estimatel.R)
x-c(115,116,98,83)
interval_estimatel(x)
mean df a b
1 103 3 78.04264 127.9574
均值为103，区间估计为（78.04264，127.9574）
source(interval_estimatel.R)
x-c(103,107,118,116)
interval_estimatel(x)
mean df a b
1 111 3 99.59932 122.4007
均值为111，区间估计为（99.59932，122.4007）
source(interval_estimatel.R)
x-c(73,89,85,97)
interval_estimatel(x)
mean df a b
86 3 70.08777 101.9122
均值为86，区间估计为（70.08777，101.9122）
c．
pairwise.t.test(X, A, p.adjust.method = none)
Pairwise comparisons using t tests with pooled SD
data: X and A
1 2
2 0.351 -
3 0.066 0.013
P value adjustment method: none
水平
均值
P
1
103
1.000
0.351
0.066
2
111
0.351
1.000
0.013
3
86
0.066
0.013
1.000
X=c(115,116,98,83,103,107,118,116,73,89,85,97)
A=factor(rep(1:3,c(4,4,4)))
pairwise.t.test(X, A, p.adjust.method = holm)
Pairwise comparisons using t tests with pooled SD
data: X and A
1 2
2 0.35 -
3 0.13 0.04
P value adjustment method: holm
X=c(115,116,98,83,103,107,118,116,73,89,85,97)
A=factor(rep(1:3,c(4,4,4)))
pairwise.t.test(X, A, p.adjust.method = bonferroni)
Pairwise comparisons using t tests with pooled SD
data: X and A
1 2
2 1.00 -
3 0.20 0.04
P value adjustment method: bonferroni
2．
lamp-data.frame(
+ X=c(20,18,18,17,15,16,13,18,22,17,26,19,26,28,23,25,24,25,18,22,27,24,12,14),
+ A=factor(rep(1:4,c(10,6,6,2)))
+ )
lamp.aov-aov(X ~ A, data=lamp)
source(anova.tab.R);anova.tab(lamp.aov)
文档评论（0）