JMP 中的 SAS 代碼?我想知道是否有辦法在 JMP 中運行 SAS 代碼(如下)。 或者我如何在 JMP 中模擬下面的 SAS 代碼。
proc mixed
Class year week treatment rep;
model stdmisummax=year|week|treatment/ddfm=kr residual;
random rep(year) treatment*rep(year);
repeated week/type=ar(1) subject=treatment*rep(year);
run;
我的因素是:治療、週(重複測量)、年份、重複(重複)。
先感謝您。
PS:我嘗試使用添加的重複測量在 JMP 中運行類似的操作,但它不會成功,因為我無法交叉三個因素來生成主題 ID。
我認為我的 JMP 版本不是專業版,或者是大學的某種專業版。
我沒有 SAS 軟件。
回复:JMP 中的 SAS 代碼?是的,隨機效應、年份(受試者)或第一個效應是受試者變異之間的估計。
回复:JMP 中的 SAS 代碼?@jiancao謝謝你的解釋。
每次我包含多個隨機效應時,我都會得到“0”(估計值)、標準誤差(參見屏幕截圖)(其中一種隨機效應)。 為什麼會這樣呢? 是因為築巢的緣故嗎?
感謝您的幫助,
伊莎多拉
回复:JMP 中的 SAS 代碼?您可以在重複測量的分析中添加多個隨機效應。 但是,您會收到一張紙條,上面寫著“重複主題列已嵌套在所有隨機效應中”。請參閱下面的屏幕截圖。 我認為這就是亞當建議放棄的原因代表(年)因為它會築巢在治療。

您是否嘗試過包含兩個隨機效應? 你看到了什麼?
回复:JMP 中的 SAS 代碼?你好@亞當莫里斯和@唐麥考馬克
我想知道您是否可以幫助我更好地理解為什麼您只使用一個術語作為隨機效應。 是不是因為我當我使用 AR1 時,是否允許在“隨機效應”下僅使用一個隨機項?
謝謝你,
伊莎多拉
回复:JMP 中的 SAS 代碼?@唐麥考馬克這非常有幫助!
感謝您與我分享,並就哪種解決方案是最佳方法向我提供意見。 我真的很感激。
伊莎多拉
回复:JMP 中的 SAS 代碼?我的解決方案如下:
Fit Model(
Y( :Response ),
Effects(
:Treatment,
:Week,
:Treatment * :Week,
:Year,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Rep[:Treatment, :Year] ),
NoBounds( 1 ),
Personality( "Mixed Model" ),
Subject( :Treatment, :Year, :Rep ),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Run( Variogram( AR1( 1 ) ) )
);
只有兩點不同。 為了獲得額外的隨機效果,我將 Rep 嵌套在 Year 中。 我使用治療、年份和重複來隱式定義主題。 不過,我發現 Adam 明確定義主體的方法在直觀上更有吸引力。
回复:JMP 中的 SAS 代碼?@亞當莫里斯,你真是太有幫助了!
非常感謝您分享您的知識和所有這些好信息! 我很感激。
回复:JMP 中的 SAS 代碼?@唐麥考馬克非常感謝您的回复。
您介意也分享您的解決方案嗎? 我想對它們進行比較,因為我有興趣了解有關 JMP 編碼和其他數據處理方法的更多信息。
謝謝你!
伊莎多拉
回复:JMP 中的 SAS 代碼?嗨伊莎多拉,
為了與閱讀此主題的其他人分享,建議的解決方案(來自 JMP 技術支持)是創建一個新的主題列,該列是標識主題的 3 個變量(治療、年份和代表)的串聯。
還需要創建 Week 列的連續副本以用作 Repeated 變量。
JMP Pro 無法擬合 SAS 代碼指定的完全相同的模型(對於 AR(1) 模型,JMP 將隨機效應嵌套在主題變量中)。但是,JMP 可以擬合省略重複(年份)項的模型。下面的腳本顯示瞭如何設置它。
謝謝,
亞當
New Column("Subject", Character, Formula( :Treatment || ", " || Char( :Year ) || ", " || Char( :Rep ) ));
New Column("Week2", numeric, continuous, formula(:Week));
Fit Model(
Y( :response ),
Effects(
:Treatment,
:Week,
:Year,
:Treatment * :Week,
:Treatment * :Year,
:Week * :Year,
:Treatment * :Week * :Year
),
Random Effects( :Treatment*:Year*:Rep ),
Subject(:Subject),
Repeated Effects( :Week2 ),
Repeated Structure( "AR(1)" ),
Personality( "Mixed Model" )
);回复:JMP 中的 SAS 代碼?伊莎多拉,
我和亞當·莫里斯討論過你的問題。 獨立地,我得出了與他相同的結果,儘管我以稍微不同的方式設置模型。 我更喜歡他的設置,所以我會使用那個。 如果您還有其他問題,請告訴我。
大學教師
回复:JMP 中的 SAS 代碼?
你好@唐麥考馬克,我附上了每個主題的數據表的屏幕截圖。
我進行了 4 次重複、6 次治療和 2 年。 在處理後的不同天收集樣品。 我在數據表中將這些採樣編號稱為採樣日,在 SAS 腳本中將其稱為“週”。
我有 JMP Pro 13。
先感謝您,
伊莎多拉
回复:JMP 中的 SAS 代碼?了解更多關於數據表的結構會很有幫助。 對於一個主題來說,它會是什麼樣子? 一周內是否會發生重複? 您需要 JMP Pro 才能使用 AR(1) 協方差結構。
回复:JMP 中的 SAS 代碼?我並不是混合效應模型方面的專家,尤其是 JMP Pro 平台方面的專家。 也許這樣的專家會加入討論。 如果沒有,您可以將其發送給 JMP 技術支持以獲得更多建議。
回复:JMP 中的 SAS 代碼?@馬克貝利感謝您分享此註釋。 我以前讀過。
我的 JMP 版本中有混合平台。 問題是我無法在這個平台中正確規定主題ID,因為我的主題是治療*代表(年)的交互。 該平台不允許我交叉和嵌套在主題 ID 框中添加主題“treatment*rep(year)”。 它允許您僅添加一個因素。
謝謝
回复:JMP 中的 SAS 代碼?對困惑感到抱歉。 不可以,無法在 JMP 中運行 SAS 程序。
JMP 能夠執行一些簡單的重複測量分析。 看到這個筆記。更複雜的情況(即不是簡單的協方差結構)需要 JMP Pro 中的混合模型平台。
回复:JMP 中的 SAS 代碼?@馬克貝利感謝您的回复。 問題是我不想在 SAS 中運行分析。 我想知道是否可以以某種方式在 JMP 中運行此代碼,或者如何在 JMP 腳本或擬合模型中模擬該代碼。
謝謝。
回复:JMP 中的 SAS 代碼?如果這適用於您的情況,您應該知道 JMP 是 SAS 的完整客戶。 正如 SAS 作為客戶端運行一樣,JMP 也可以。 它需要 SAS 服務器連接。 該服務器可能在同一台計算機上運行,就像我的工作筆記本電腦一樣,也可能通過網絡連接遠程運行。 您可以訪問遠程 SAS 服務器(公司、大學)嗎? 如果是這樣,請將 SAS 代碼粘貼到 JMP 腳本編輯器中。 在編輯器中右鍵單擊並選擇提交 SAS。
選擇“文件”>“首選項”並首先設置與 SAS 服務器的連接。
回复:JMP 中的 SAS 代碼?@uday_guntupalli感謝您的回复。 是的,我知道,但我不想進入 SAS。 我從一位同事那裡得到了這段代碼,因為我無法或無法找出如何使用 JMP 中規定的因素運行相同的分析。
謝謝 :)
回复:JMP 中的 SAS 代碼?@博迪尼,
我知道您說您沒有 SAS,但您知道可以免費下載 SAS 大學版嗎? 如果您確實需要測試 SAS 代碼,您可以下載 SAS 大學版並嘗試一下。
https://www.sas.com/en_us/software/university-edition/download-software.html
我不知道如何在 JMP 中執行等效步驟,所以我也會關注其他答案:)