模擬實驗我想做一個模擬實驗，生成12個隨機變量，每個變量包含100個觀測值，這樣變量就受到異質性問題，並且服從正態分佈
該實驗重複 10,000 次。
# 我不知道如何使用 JMPpro 做到這一點
謝謝
高級統計建模視窗回复：模擬實驗“如何使變量的方差異質” - 這是否意味著每列需要不同的西格瑪？ 您需要 12 個隨機選擇的西格瑪值嗎？
“異方差問題” - 這是否意味著您有一個基於 mu 改變 sigma 的函數？ 還有一個選擇 mu 的函數，可能是隨機的？ （我沒有任何專業知識，只是讀了 30 秒後猜測。）
我想你會想使用嵌套的 for 循環來做到這一點。 我原來的建議太簡單了。 您還可以使用列公式，但 JSL 腳本更易於管理。 嵌套循環看起來像這樣
m = J( 100 * 10000, 12, 0 );
for(row=1, row<=nrows(m), row+=1,
for(col=1, col<=ncols(m), col+=1,
m[row,col] = HSD(row,col);
);
);<BR /><BR />dt = astable(m)<BR /><BR />...
您必須用您實際想要的內容編寫“HSD”用戶編寫的函數。
您可能還想將 10000 次重複更改為 2 或 3 次重複，直到您決定分析應該做什麼。 我使用了 Summary，因為它將所有 10K 分析打包到一個數據表中。 您可能想要運行一個平台 10K 次，獲取一些結果，然後關閉該平台。 您不希望一次打開 10K 平台，因為 JMP 會耗盡內存。
回复：模擬實驗尊敬的先生
首先，非常感謝您的幫助。
事實上，你上面的腳本很接近我想要的，
但我不知道如何使變量的方差異質。
我不知道如何修改你的腳本，例如：
隨機變量數 = 12
每個變量的觀測值數量 = 100
由此產生的變量正遭受異方差問題
可重複實驗10,000次
問候
回复：模擬實驗看看這樣的搜索https://www.google.com/search?q=site%3Ajmp.com+Heteroscedasticity查找以前在社區中可能有好的答案的問題。 在谷歌搜索 GUI 中你可以這樣做網站：jmp.com限制結果：
site:jmp.com Heteroscedasticitysite:jmp.com Heteroscedasticity網站：jmp.com 異方差
我不是統計學家。 上面的數據是隨機的。 您指出的來源表明數據應該在各個季節有某種重複......而不是完全隨機的。 如果 1..12 代表一年中的月份，那麼擁有一月列、二月列……可能沒有意義。 擁有日期列和值列可能更有意義。
| 日期 |
價值 |
| 2020/01 |
7 |
| 2020/02 |
6 |
| 2020/03 |
8 |
我抬頭異方差性並懷疑您需要一種不同的方法來生成隨機數，以便它們實際上表現出該特徵。 否則它只是隨機的，沒有什麼有趣的，但 mu 和 sigma 傳遞到隨機法線，如上圖所示。
只是猜測，並且出路......我認為你必須生成一個隨機 mu （正常？均勻？），以及基於該 mu 的 sigma，然後使用 RandomNormal(mu,sigma)。 或者我可能完全錯了。
回复：模擬實驗
來源：
隨機係數週期性自回歸 (repec.org)
回复：模擬實驗尊敬的先生克雷格·黑爾斯
非常感謝你對我的幫助，
我想做一個模擬實驗：
隨機變量數 = 12
每個變量的觀測值數量 = 100
由此產生的變量正遭受異方差問題
# 可以重複實驗10,000次
回复：模擬實驗總結可能是一個好的開始
</CODE></P>
<PRE><CODE class=" language-jsl">dt = As Table( J( 100 * 10000, 12, Random Normal( 0, 1 ) ) );
dt << New Column( "group", formula( Floor( (Row() - 1) / 100 ) + 1 ) );
dt << runformulas;
dt:group << deleteformula;
dt << Go To( group );
dt << Move Selected Columns( To first );
dtsummary = dt << Summary(
Group( :group ),
Mean( :Col1 ),
Min( :Col1 ),
Max( :Col1 ),
Std Dev( :Col1 ),
Mean( :Col2 ),
Min( :Col2 ),
Max( :Col2 ),
Std Dev( :Col2 ),
Freq( "None" ),
Weight( "None" )
);
dtsummary << Graph Builder(
Size( 653, 452 ),
Show Control Panel( 0 ),
Variables(
X( :group ),
Y( :"Mean(Col1)"n ),
Y( :"Min(Col1)"n, Position( 1 ) ),
Y( :"Max(Col1)"n, Position( 1 ) ),
Y( :"Std Dev(Col1)"n, Position( 1 ) ),
Y( :"Mean(Col2)"n, Position( 1 ) ),
Y( :"Min(Col2)"n, Position( 1 ) ),
Y( :"Max(Col2)"n, Position( 1 ) ),
Y( :"Std Dev(Col2)"n, Position( 1 ) )
),
Elements(
Points(
X,
Y( 1 ),
Y( 2 ),
Y( 3 ),
Y( 4 ),
Y( 5 ),
Y( 6 ),
Y( 7 ),
Y( 8 ),
Legend( 7 )
)
)
);
</CODE></P>
<P><SPAN class="lia-inline-image-display-wrapper lia-image-align-inline" image-alt="Color!" style="width: 812px;"><SPAN class="lia-inline-image-display-wrapper" image-alt="Color!" style="width: 812px;"><span class="lia-inline-image-display-wrapper" image-alt="Color!" style="width: 812px;"><img src="https://us.v-cdn.net/6038925/uploads/images/41991iFA52067AFA73F318/Untitled%204%20By%20%28group%29%20-%20Graph%20Builder.png" role="button" title="Untitled 4 By (group) - Graph Builder.png" alt="Color!" /><span class="lia-inline-image-caption" onclick="event.preventDefault();">Color!</span></span><SPAN class="lia-inline-image-caption" onclick="event.preventDefault();">Color!</SPAN></SPAN><SPAN class="lia-inline-image-caption" onclick="event.preventDefault();">顏色！</SPAN></SPAN></P></LINGO-BODY><LINGO-SUB id="lingo-sub-481697" slang="en-US" mode="NONE">回复：模擬實驗</LINGO-SUB><LINGO-BODY id="lingo-body-481697" slang="en-US" mode="NONE"><P>我不知道異質性問題是什麼，但其餘的可能看起來像這樣：</P>
<P> </P>
<PRE><CODE class=" language-jsl">dt = As Table( J( 100 * 10000, 12, Random Normal( 0, 1 ) ) );
dt << New Column( "group", formula( Floor( (Row() - 1) / 100 ) + 1 ) );
dt << runformulas;
dt:group << deleteformula;
dt << Go To( group );
dt << Move Selected Columns( To first );
使用組變量進行分析。
12 million random normal numbers12 million random normal numbers1200 萬個隨機正態數
