問題：我們觀察到使用不同數據子集運行時，分裂圖模型的行為不穩定。 對於六個響應變量列中的兩個，腳本為模型中的隨機因子生成錯誤的輸出。 具體來說，隨機因素佔REML方差分量估計表中總方差的0％，而最小二乘均值表顯示的是一列均值，而不是標準誤。
我們已經檢查了列的屬性，並替換了有問題的列中的數據，但無濟於事。 有兩行證據表明此行為不是由於腳本錯誤引起的：
- 該腳本可在其他四個響應變量列上正常運行
- 簡化的模型（僅包含隨機因子和有問題的數據列）為一個響應變量生成了錯誤的輸出，而對於另一個變量，輸出似乎是正確的
我們提供了實驗設計和分析方法的簡短摘要，以及有問題的數據，使用的腳本和模型輸出的示例。 我們還附帶了一個pdf文件，其中顯示了上述簡化模型的腳本和輸出。
感謝您為解決這個難題提供的幫助！
實驗設計：一個生態野外實驗，由四個（或三個）空間定義的塊（重複項）組成，在分塊圖設計中隨機分配了三個因素。 在整個地塊上，
將兩個相對較難更改的二級因子-帽毛去除（+/-）和水分增加（+/-）-隨機分配給每個塊中的四個圖。 在子圖規模上，將刀片移除（Blades +/-）隨機分配給每個整個圖中的兩個子圖（請參見下面的示例佈局）。
樣本數據集
大塊＃ | 跛行 | 濕 | 葉片 | 傑88 | 先生89 |
1個 | + L | + M | + B | 0 | -15 |
1個 | -L | + M | + B | -5.88235 | -2.94118 |
1個 | + L | -M | + B | 5.88235 | -11.7647 |
1個 | -L | -M | + B | -14.4118 | -38.2353 |
2 | + L | + M | + B | -2.94118 | 2.94118 |
2 | -L | + M | + B | 0 | -20.5882 |
2 | + L | -M | + B | 0 | -15 |
2 | -L | -M | + B | 0.29412 | -2.94118 |
3 | + L | + M | + B | 0 | 。 |
3 | -L | + M | + B | -32.3529 | 。 |
3 | + L | -M | + B | -17.6471 | 。 |
3 | -L | -M | + B | -15 | 。 |
4 | + L | + M | + B | 0 | -2.64706 |
4 | -L | + M | + B | -5.88235 | -14.4118 |
4 | + L | -M | + B | -5.58824 | -43.8235 |
4 | -L | -M | + B | -5.88235 | -35.2941 |
1個 | + L | + M | -B | -5.88235 | 11.47059 |
1個 | -L | + M | -B | 3.23529 | 5.29412 |
1個 | + L | -M | -B | -2.94118 | -8.82353 |
1個 | -L | -M | -B | -2.94118 | -11.7647 |
2 | + L | + M | -B | -5.88235 | 6.17647 |
2 | -L | + M | -B | 0 | -6.17647 |
2 | + L | -M | -B | 2.94118 | -26.4706 |
2 | -L | -M | -B | 3.23529 | -20.2941 |
3 | + L | + M | -B | -14.7059 | 。 |
3 | -L | + M | -B | 14.70588 | 。 |
3 | + L | -M | -B | -0.29412 | 。 |
3 | -L | -M | -B | -5.88235 | 。 |
4 | + L | + M | -B | -8.82353 | -26.4706 |
4 | -L | + M | -B | -26.4706 | -20.5882 |
4 | + L | -M | -B | -14.4118 | -44.1177 |
4 | -L | -M | -B | 2.94118 | -47.3529 |
分析。 為了簡化分析，我們選擇將四種全圖處理組合建模為單個因子的四個級別（以下腳本中的Concatenate [Limp，Moist]）。 我們計劃使用事後檢驗來區分治療手段。 隨機效應（整個圖的誤差項）是塊，帽貝和濕度代碼（例如，下面的：Name（“ Concatenate [Blk＃，Limp，Moist]”）[：Blk＃]）的串聯。
完整模型的示例腳本。
</CODE>Fit Model(<BR />Y( :Ja88Ap88 ),<BR />Effects(<BR />:Blk#,<BR />:Name( "Concatenate[Limp,Moist]" ),<BR />:Blades,<BR />:Blades * :Name( "Concatenate[Limp,Moist]" )<BR />),<BR />Random Effects( :Name( "Concatenate[Blk#,Limp,Moist]" )[:Blk#] ),<BR />NoBounds( 1 ),<BR />Personality( "Standard Least Squares" ),<BR />Method( "REML" ),<BR />Emphasis( "Minimal Report" ),<BR />Run(<BR />)<BR />),</PRE><P> </P><P><U>完整模型的部分輸出</U>。</P><P><STRONG>REML方差分量估計</STRONG>對於Ja88Apr88（n = 4個塊）</P><TABLE><TBODY><TR><TD><P><STRONG>隨機效應</STRONG></P></TD><TD><P><STRONG>可變比率</STRONG></P></TD><TD><P><STRONG>無功分量</STRONG></P></TD><TD><P><STRONG>標準錯誤</STRONG></P></TD><TD><P><STRONG>降低95％</STRONG></P></TD><TD><P><STRONG>上層95％</STRONG></P></TD><TD><P><STRONG>Wald p值</STRONG></P></TD><TD><P><STRONG>總計百分比</STRONG></P></TD><TD> </TD></TR><TR><TD><P>串聯[Blk＃，Limp，潮濕] [Blk＃]</P></TD><TD><P>-0.324822</P></TD><TD><P>-40.38441</P></TD><TD><P>27.376421</P></TD><TD><P>-94.04121</P></TD><TD><P>13.272392</P></TD><TD><P>0.1402</P></TD><TD><P>0.000</P></TD><TD> </TD></TR><TR><TD><P>剩餘的</P></TD><TD><P> </P></TD><TD><P>124.32778</P></TD><TD><P>50.756605</P></TD><TD><P>63.930876</P></TD><TD><P>338.78407</P></TD><TD><P> </P></TD><TD><P>100.000</P></TD><TD> </TD></TR><TR><TD><P>總</P></TD><TD><P> </P></TD><TD><P>124.32778</P></TD><TD><P>50.756605</P></TD><TD><P>63.930876</P></TD><TD><P>338.78407</P></TD><TD><P> </P></TD><TD><P>100.000</P></TD><TD> </TD></TR></TBODY></TABLE><P> </P><P><STRONG>固定效果測試</STRONG>對於Ja88Apr88<U>注意</U>分母的自由度似乎是正確的。</P><TABLE><TBODY><TR><TD><P><STRONG>資源</STRONG></P></TD><TD><P><STRONG>納帕姆</STRONG></P></TD><TD><P><STRONG>東風</STRONG></P></TD><TD><P><STRONG>DFDen</STRONG></P></TD><TD><P><STRONG>F比</STRONG></P></TD><TD><P><STRONG>概率> F</STRONG></P></TD><TD> </TD></TR><TR><TD><P>大塊＃</P></TD><TD><P>3</P></TD><TD><P>3</P></TD><TD><P>9</P></TD><TD><P>3.1201</P></TD><TD><P>0.0808</P></TD><TD> </TD></TR><TR><TD><P>級聯[Limp，潮濕]</P></TD><TD><P>3</P></TD><TD><P>3</P></TD><TD><P>9</P></TD><TD><P>0.2215</P></TD><TD><P>0.8791</P></TD><TD> </TD></TR><TR><TD><P>葉片</P></TD><TD><P>1個</P></TD><TD><P>1個</P></TD><TD><P>12</P></TD><TD><P>0.3675</P></TD><TD><P>0.5557</P></TD><TD> </TD></TR><TR><TD><P>刀片*級聯[Limp，潮濕]</P></TD><TD><P>3</P></TD><TD><P>3</P></TD><TD><P>12</P></TD><TD><P>1.0059</P></TD><TD><P>0.4238</P></TD><TD> </TD></TR></TBODY></TABLE><P> </P><P><STRONG>最小二乘均值表</STRONG>對於Ja88Apr88（n = 4個塊）</P><TABLE><TBODY><TR><TD><P><STRONG>水平</STRONG></P></TD><TD><P><STRONG>最小均方</STRONG></P></TD><TD><P><STRONG>標準錯誤</STRONG></P></TD><TD> </TD></TR><TR><TD><P>[1] 1，+ L，+ M</P></TD><TD><P>-2.28820</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[1] 1，+ L，-M</P></TD><TD><P>-9.03709</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[1] 1，-L，+ M</P></TD><TD><P>-8.62806</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[1] 1，-L，-M</P></TD><TD><P>8.48276</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[2] 2，+ L，+ M</P></TD><TD><P>7.78410</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[2] 2，+ L，-M</P></TD><TD><P>-1.69162</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[2] 2，-L，+ M</P></TD><TD><P>-3.73673</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[2] 2，-L，-M</P></TD><TD><P>-3.53222</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[3] 3，+ L，+ M</P></TD><TD><P>-11.31712</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[3] 3，+ L，-M</P></TD><TD><P>-6.88603</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[3] 3，-L，+ M</P></TD><TD><P>-11.93065</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[3] 3，-L，-M</P></TD><TD><P>-5.45444</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[4] 4，+ L，+ M</P></TD><TD><P>-14.25232</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[4] 4，+ L，-M</P></TD><TD><P>-2.45880</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[4] 4，-L，+ M</P></TD><TD><P>4.22192</P></TD><TD><P>。</P></TD><TD> </TD></TR><TR><TD><P>[4] 4，-L，-M</P></TD><TD><P>-19.56963</P></TD><TD><P>。</P></TD><TD> </TD></TR></TBODY></TABLE><P> </P></LINGO-BODY><LINGO-LABS id="lingo-labs-328571" slang="en-US"><LINGO-LABEL>數據表</LINGO-LABEL><LINGO-LABEL>腳本編寫</LINGO-LABEL></LINGO-LABS><LINGO-SUB id="lingo-sub-328695" slang="en-US">回复：分割圖模型的錯誤行為：受污染的數據，腳本錯誤或錯誤？</LINGO-SUB><LINGO-BODY id="lingo-body-328695" slang="en-US"><P>讓我作為我的評論的開頭，因為我不從科學或工程的角度理解您所描述的情況，因此我的評論可能完全不合時宜。 確實很難改變帽貝和濕度，這是隨機的。 我不知道您的假設是什麼，或者我不了解響應變量或與此相關的因素？ 關於您的實驗，我有一些疑問。 我遇到的困難是理解您的設計計劃（基於您提供的樣本數據）：</P><P>1.這句話對我來說毫無意義：<SPAN>在分割圖設計中隨機分配了3個因子”。 分割圖設計故意限制了隨機性，所以我不明白您在做什麼。 在樣本數據集中，看起來“您很難更改</SPAN>演員”正在發生很大變化！ 子圖因子似乎沒有像您的整個圖因子那樣變化？ 您有運行訂單嗎？</P><P>2.我也不理解您對積木的使用？ 我假設您正在使用這些塊來獲取整個繪圖錯誤的估計值？ 讓我看看我是否了解您在做什麼...您正在4個不同的“位置/空間”（塊）中運行2 ^ 3階乘實驗。 在每個塊內，您打算將隨機數限制在階乘矩陣中的3個因子（L和M）中的2個。 這聽起來更像是分割圖。 我只是在jmp文件中添加了設計，以查看這是否是您正在運行的。</P><P>3.為什麼要合併L＆M？ 如果執行此操作，您將無法量化L * M相互作用以及L和M的其他任何二階效應的能力？ 相反，您可以估計二次和三次效應？？？</P><P>4.您的測量系統看起來有一個辨別問題（它似乎不是連續的）。 查看樣本數據集時只有一些數字可用？ 這是真實數據嗎？</P><P>值得深思（請原諒我的漫步）...之所以需要注意並考慮隨機限制，是因為在整個實驗過程中，噪聲可能會發生意料之外的變化，進而影響實驗誤差和隨後進行統計檢驗。 思想是整個圖錯誤將不同於子圖錯誤，因此應針對每個圖執行不同的顯著性檢驗。 如果在計劃實驗時，如果您確定了噪聲變量並假設了這些噪聲變量的影響，則可以在設計實驗時考慮到這一點（例如，塊內部發生了什麼變化以及塊之間發生了什麼變化？） 。 如果在整個圖塊級別（在塊之間）變化的噪聲與在子電位差級別變化的噪聲之間幾乎沒有差異，則考慮該限制可能沒有影響。 但是，如果整個圖級的噪聲（以及估計誤差）與子圖級的噪聲存在顯著差異，則顯著性檢驗應反映出這一點。</P><P> </P></LINGO-BODY><LINGO-SUB id="lingo-sub-328867" slang="en-US">回复：分割圖模型的錯誤行為：受污染的數據，腳本錯誤或錯誤？</LINGO-SUB><LINGO-BODY id="lingo-body-328867" slang="en-US"><P>請閱讀<A href="https://community.jmp.com/t5/user/viewprofilepage/user-id/4358" target="_blank">@statman</A>的答复仔細。</P>
<P> </P>
<P>您的帖子也讓我感到困惑：</P>
<UL>
<LI>您已將Limpet和Moisture水平串聯起來，然後將它們的關節變化視為對響應的隨機影響？ 您不希望這些因素產生固定的影響嗎？</LI>
<LI>您所說的Blk＃實際上就是整個情節？ 也就是說，這些運行一次重置難以改變的因素？ 那麼它可能應該是模型中唯一的隨機效應。</LI>
<LI>似乎Blades是整個圖（很難更改）的因素，但是在同一整個圖中，如何將Blades設置為-B和+ B呢？</LI>
<LI>您的腳本表明固定和隨機效果都使用相同的因子。 那是故意的嗎？</LI>
<LI>您丟失了全部情節= 3？</LI>
</UL>
<P>也許此模型將以更穩定的估計提供您想要的信息：</P>
<P> </P>
<PRE><CODE class=" language-jsl">Fit Model(
Y( :Ja88Ap88 ),
Effects( :Limp, :Moist, :Blades, :Blades * :Limp, :Blades * :Moist ),
Random Effects( :Blk# ),
No Bounds( 1 ),
Personality( "Standard Least Squares" ),
Method( "REML" ),
Emphasis( "Minimal Report" ),
Run()
);
