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 中执行等效步骤,所以我也会关注其他答案:)