Anova-n输出满秩。返回NaN假定值。

23日视图(30天)

显示旧的评论

运气哈维兰德 2022年12月12日

0
链接

这个问题直接联系

//www.tianjin-qmedu.com/matlabcentral/answers/1876737-anova-n-outputs-as-not-full-rank-returns-nan-p-value

评论道: 运气哈维兰德2022年12月13日

答:接受杰夫•米勒

我想分析一组感兴趣的标本用三种不同的工艺参数。停留时间(按多长时间和加热),最大温度实现,在制造和最大压力。为例。样品可以用168摄氏度,300秒的停留时间,和483 kPa的压力,或204摄氏度,800秒的停留时间,和689 kPa的压力,或204摄氏度,300秒的停留时间,和483 kPa的压力,等等。

这些参数将对试样的强度有影响。每个试样的强度被发现通过实验测试和“y0”中列出。我想确定哪些工艺参数影响强度(强度越高越好)。

ANOVA-N当我运行代码,它返回NaN的假定值说没有我的测量是满秩。

gX1 = [“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”]

gX1 =1×30字符串数组

“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“168”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”“204”

gX2 = [“483”“483”“483”“483”“483”“483”“483”“689”“689”“689”“689”“689”“689”“689”“689”“483”“483”“483”“483”“483”“483”“483”“689”“689”“689”“689”“689”“689”“689”“689”]

gX2 =1×30字符串数组

“483”“483”“483”“483”“483”“483”“483”“689”“689”“689”“689”“689”“689”“689”“689”“483”“483”“483”“483”“483”“483”“483”“689”“689”“689”“689”“689”“689”“689”“689”

gX3 = [“300”“300”“300”“300”“300”“300”“300”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“300”“300”“300”“300”“300”“300”“300”“300”]

gX3 =1×30字符串数组

“300”“300”“300”“300”“300”“300”“300”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“600”“300”“300”“300”“300”“300”“300”“300”“300”

y0 = (18.2919831666125;

12.5641727939410;

15.3621391423407;

14.2489632470725;

13.3879424625516;

16.5153213712401;

16.4715840677642;

23.1366242658033;

22.3254273493200;

22.5234847436383;

23.4593859332890;

23.1149682822853;

24.1407966949033;

25.0761899528730;

23.8484621339037;

19.8895058992870;

18.5343968548295;

19.8404597873812;

19.0451878113906;

17.3858461481879;

17.0566716618079;

17.2640841888405;

22.7626298124232;

22.8208067604289;

21.9108601449323;

20.8905709384620;

21.5143053296674;

19.6368344554986;

18.3404852458128;

18.2286716918400;】

y0 = 30×1

18.2920 12.5642 15.3621 14.2490 13.3879 16.5153 16.4716 23.1366 22.3254 22.5235

(yy, zz, Xstats) = anovan (y0, {gX1、gX2 gX3},“模型”,“互动”,“varnames”,{“临时”,“压力”,“时间”});

我怎么解决这个问题?

谢谢,

运气

0评论
显示隐藏1以上的评论

登录置评。

在回答这个问题。

接受的答案

杰夫•米勒 2022年12月13日

0
链接

直接链接到这个答案

//www.tianjin-qmedu.com/matlabcentral/answers/1876737-anova-n-outputs-as-not-full-rank-returns-nan-p-value answer_1126657

简短的回答是,你必须收集更多的数据来解决这个问题。看这里:

                              > >交叉表(gX1 gX2 gX3)
                             
                              ans (:: 1) =
                             
                              7 0
                             
                              0 8
                             
                              ans (:: 2) =
                             
                              0 8
                             
                              7 0

你只有4中数据的8细胞2 x2x2设计。这样的安排,每个因素的主要作用是困惑的双向互动的剩下的两个因素,所以没有办法估计所有的主要影响和双向互动的同时。以矩阵计算,预测矩阵包括所有这些条款不是满秩。

如果你不想收集更多的数据,那么你只能估计一个子集的主要效应和双向互动的条件。这样做的方法之一是遗漏因素——只是分析他们两个。另一种方法就是使用条款矩阵来表示因素和交互包括与忽视。但与这些选项中的任何一个,你应该记住,你们遗漏的影响方式仍在那里,导致包括条款的明显影响。

3评论
显示隐藏特殊照顾年长的评论

运气哈维兰德 2022年12月13日

谢谢杰夫。你一个真正的一个。

给你的建议: