UA MATH571B 试验设计 QE练习题1

  • 第一题
  • 第二题
  • 第三题

这是2014年一月Methodology的1-3题。

第一题

希望今年的考试没有这种概念题!

第一个是对的,比较多正态总体的方差就是FFF检验做的事情。
第二个不对,选择Bonferroni方法还是Scheffe方法取决于对多少个统计量做联合推断,如果是对比较多的统计量做联合推断就用Scheffe方法;否则就用Bonferroni方法。
第三个不对,Type I SS是sequential SS,要考虑进入模型的顺序;Type III SS才是不考虑顺序的。
第四个不对,虽然满足了ar=bkar=bkar=bk的条件,但这只是必要条件。
第五个是对的,这是教材上的原话Pg304:
In many situations it is impossible to perform all of the runs in a 2 k factorial experiment under
homogeneous conditions. For example, a single batch of raw material might not be large
enough to make all of the required runs. In other cases, it might be desirable to deliberately
vary the experimental conditions to ensure that the treatments are equally effective (i.e.,
robust) across many situations that are likely to be encountered in practice. For example, a
chemical engineer may run a pilot plant experiment with several batches of raw material
because he knows that different raw material batches of different quality grades are likely to
be used in the actual full-scale process.
The design technique used in these situations is blocking.


这一张里面的第一个是错的,whole-plot factor是hard-to-change的,sub-plot factor是easy-to-change的;
关于Graeco-Latin设计的前三个都是错的,第四个是对的。Graeco-Latin设计应对的是有三个nuisance factor的情况,Latin Square设计才是两个;a=5a=5a=5对应的总试验单位数目与Latin Square相同,仍然是525^252;残差的自由度是(a−1)(a−1−1−1)=8(a-1)(a-1-1-1)=8(a−1)(a−1−1−1)=8,三个1分别是贡献给三个nuisance factor的;因此对model的F检验自由度是(52−1−8,8)=(16,8)(5^2-1-8,8) = (16,8)(52−1−8,8)=(16,8)。这里附一个Graeco-Latin的ANOVA table帮助理解:

(h)是正确的,确定confounding factor的时候用的factor越多越好,所以从I=ABCDI = ABCDI=ABCD开始;
(i)不对,根据选最长的原理,应该是选择4,4 ,6那个。

第二题


第一小问的提示在Computer Output里面,Job和Operator都是random factor,并且Operator标识是Operator(Job)说明Operator这个factor是嵌套在Job这个factor中的,因此这是一个Nested Design。

第二小问:Statistical Model of this nested design is
yijk=μ+τi+βj(i)+ϵ(ij)kϵ(ij)k∼iidN(0,σ2);i=1,⋯,6;j=1,⋯,3;k=1,⋯,2y_{ijk} = \mu + \tau_i + \beta_{j(i)} + \epsilon_{(ij)k}\\ \epsilon_{(ij)k} \sim_{iid}N(0,\sigma^2);i=1,\cdots,6;j=1,\cdots,3;k=1,\cdots,2 \\ yijk​=μ+τi​+βj(i)​+ϵ(ij)k​ϵ(ij)k​∼iid​N(0,σ2);i=1,⋯,6;j=1,⋯,3;k=1,⋯,2

where τi\tau_iτi​ indicates treatment effect of job iii and βj(i)\beta_{j(i)}βj(i)​ indicates treatment effect of operator jjj nested in job iii. (参考UA MATH571B 试验设计VI 随机效应与混合效应1写出假设)
τi∼iidN(0,στ2),βj(i)∼iidN(0,σβ2)\tau_i \sim_{iid} N(0,\sigma_{\tau}^2),\beta_{j(i)} \sim_{iid} N(0,\sigma^2_{\beta})τi​∼iid​N(0,στ2​),βj(i)​∼iid​N(0,σβ2​)
第三小问:将ANOVA table补充完整;
Job的自由度为6−1=56 - 1=56−1=5,Operator的自由度为6×(3−1)=126 \times (3-1)=126×(3−1)=12,残差的自由度是3×6×(2−1)=183 \times 6 \times (2-1)=183×6×(2−1)=18;剩下的问题就很简单了

第四小问:要计算Job的treatment factor的方差,可以使用教材上的公式(14.12)


第五小问:要计算operator的方差,可以使用教材上的公式(14.11)


第六小问:我们可以看一下第三小问算出来的p值,然后得出结论,
P value of factor job is near zero, so we can reject null hypothesis that variability of effect of job is zero and conclude that the effects of jobs differ significantly. While p-value of factor operator nested in job is greater than .05, so we conclude the effects of different operators working for the same job are not significantly different.

第三题




第一问:两因子的析因设计自由度分析可以参考这张图

因此压强的自由度是4−1=34-1=34−1=3,温度的是3−1=23-1=23−1=2,第一个glm procedure没有用交互项,因此残差的自由度是4×3−1−2−3=64\times 3 - 1 - 2 -3 = 64×3−1−2−3=6

第二问:有两个及以上因子的析因设计都需要检验因子效应的可加性,第二个glm procedure就是Tukey’s 检验的简易操作,用来检验可加性,SAS output的结果是q的p值超过了.05,因此拒绝原假设,得出因子效应具有可加性的结论,所以答案是
check the non-additivity of the model

第三问:这个我就不写了,参考课件Topic 11第29页,第四问问检验的名字可以说是Tukey’s test of nonadditivity也可以说Tukey’s one degree of freedom test,但不能只说Tukey’s test,因为treatment effect的配对比较还有Tukey检验和Tukey-Kramer检验。


第五问&第六问:看那个SAS output,有q的时候两个效应都是不显著的,并且q的p值说明因子是可加的,因此我们可以只关注第一个glm procedure的结果即可。

UA MATH571B 试验设计 QE练习题1相关推荐

  1. UA MATH571B 试验设计 QE练习题 不使用代码分析试验结果I

    UA MATH571B 试验设计 QE练习题 不使用代码分析试验结果I 2014年5月第一题 2015年5月第一题 2016年5月第二题 2017年1月第一题 不使用代码分析试验结果考察的是对试验设计 ...

  2. UA MATH571B 试验设计 QE练习题 平衡非完全区组设计BIBD

    UA MATH571B 试验设计 QE练习题 平衡非完全区组设计BIBD 2015年1月第一题 2018年1月第一题 2015年1月第一题 Part a 显然这是一个BIBD(Balanced Inc ...

  3. UA MATH571B 试验设计 2k析因设计理论上

    UA MATH571B 试验设计 2k析因设计理论上 2k2^k2k析因设计 因子效应的计算 ANOVA table Single-replicate 2k2^k2k析因设计 2k2^k2k析因设计( ...

  4. UA MATH571B 试验设计 Quarter 2-level析因设计

    UA MATH571B 试验设计 Quarter 2-level析因设计 2k−22^{k-2}2k−2设计的基本概念 2k−22^{k-2}2k−2试验结果的SAS分析 Quarter 2-leve ...

  5. UA MATH564 概率论 QE练习题 信封问题

    UA MATH564 概率论 QE练习题 信封问题 2015年1月的第二题 2015年5月的第一题 这一篇介绍QE理论中出现了两个信封问题相关的题目:2015年1月的第二题和2015年5月的第一题.信 ...

  6. UA MATH571B 试验设计 总结 判断试验类型的例题

    UA MATH571B 试验设计 总结 判断试验类型的例题 2015年5月第二题:Blocking factorial design 2015年5月第三题:Split-plot Design 2017 ...

  7. UA MATH571B 试验设计 总结 试验的类型与选择

    UA MATH571B 试验设计 总结 试验的类型与选择 试验的意义,试验设计的意义 单个treatment factor的试验 多个treatment factor的试验 多个treatment f ...

  8. UA MATH571B 试验设计 2k析因设计理论下

    UA MATH571B 试验设计 2k析因设计理论下 Confounding Blocking 222^222析因设计 Blocking 232^323析因设计 principal block Fra ...

  9. UA MATH564 概率论 QE练习题6

    UA MATH564 概率论 QE练习题6 第二题 第三题 这是2016年5月的2-3题. 第二题 Part a Denote O as center of the circle. Randomly ...

最新文章

  1. 简单介绍实体类或对象序列化时,忽略为空属性的操作
  2. 在Docker中体验数据库之MySql
  3. 2.4 Add a dynamic map
  4. Vue使用axios,设置axios请求格式为form-data
  5. MATLAB Robotic Toolbox 机器人工具箱示例
  6. ArcGIS License Server Administrator 10.2 无法启动许可的解决办法
  7. mysql查询所有姓王的信息_MySQL的查询练习
  8. 太原警方通过微博提醒您手机丢失如何保微信安全
  9. 篮球计时计分器c语言程序,C51的篮球计时计分器程序
  10. 中国石油大学c语言程序设计答案,中国石油大学C语言程序设计第一册在线作业答案...
  11. 【MTK AF】Acce/Gyro/PD/Laser Driver Check
  12. 金仓数据库KingbaseES blob类型数据导入导出
  13. javaScript中什么时候用分号;
  14. lua的演进 lua的历史
  15. CROSSFORMER: A VERSATILE VISION TRANSFORMER BASED ON CROSS-SCALE ATTENTION
  16. java.lang.exceptionininitializererror_解决出现 java.lang.ExceptionInInitializerError错误问题
  17. 三星android文件传输,最好的三星Galaxy S8管理器:如何将文件传输到三星Galaxy S8...
  18. 四十七、Fluent近壁面处理
  19. Docker复习07——Docker 网络
  20. mybatis报错:java.sql.SQLException: 无效的列类型: 1111

热门文章

  1. Leetcode 136. 只出现一次的数字 解题思路及C++实现
  2. Ethernet/IP 学习笔记五
  3. ubuntu11.10 安装tpc_w
  4. JVM锁和分布式锁是什么关系
  5. JdbcTemplate的CRUD操作
  6. Python 技术篇-基于PyHook3+threading多线程实现鼠标单击事件和双击事件的识别实例演示
  7. C语言用DFS实现找到图的所有路径(邻接矩阵实现)
  8. 如何:为 Windows 窗体 RichTextBox 控件设置字体属性
  9. PCB布线技术 很好很强大
  10. CTFshow php特性 web93