复发分析累积公式（原为JSL编码）大家好
我想保存复发分析模型中每个原因级别的累积公式。以下是我的代码。
ra = Recurrence Analysis(
Y( :age ),
Label( :BOT ID ),
Cost( :Censor ),
Grouping( :Bot Series ),
Event Plot( 0 ),
Fit Model(
Scale Effects( :Bot Series ),
Shape Effects( :Bot Series ),
Run Model,
Model Type( "Proportional Intensity Poisson Process" )
)
);
// Get fitted recurrence model object
fit = Report(ra)[OutlineBox("Fitted Recurrence Model")] << Get Scriptable Object;
// Save intensity and cumulative formulas from fit
intensity_col = fit << Save Intensity Formula;
cumul_col = fit << Save Cumulative Formula (
If( :Part Grouping2 == "Caster Assy-Left"",
full_formula
)
);
I have Part Grouping2 as Cause and there are 13 levels in Part Grouping2, I want to generate fit value for each level in Part Grouping2. I tried some code but did not work. In the code Part Grouping2 == "Caster Assy-Left" . That is just trying to use the formular for one Part Grouping2 level Caster Assy-Left. There are 13 total levels and I want to automate this so I don't need to repeat 13 times.
Thanks.
</CODE></P>
<DIV>
</DIV>
</LINGO-BODY>
<LINGO-LABS id="lingo-labs-912116" slang="en-US" mode="UPDATE">
<LINGO-LABEL>自动化和脚本编写</LINGO-LABEL>
</LINGO-LABS>
<LINGO-SUB id="lingo-sub-912154" slang="en-US" mode="CREATE" translate="no">回复：JSL编码</LINGO-SUB>
<LINGO-BODY id="lingo-body-912154" slang="en-US" mode="CREATE">
<P>你最终想要提取或实现什么？你的脚本会为你生成公式，公式看起来大致如下：</P>
<P><SPAN class="lia-inline-image-display-wrapper lia-image-align-inline" image-alt="jthi_0-1762581573217.png" style="width: 999px;"><SPAN class="lia-inline-image-display-wrapper" image-alt="jthi_0-1762581573217.png" style="width: 999px;"><SPAN class="lia-inline-image-display-wrapper" image-alt="jthi_0-1762581573217.png" style="width: 999px;"><SPAN class="lia-inline-image-display-wrapper" image-alt="jthi_0-1762581573217.png" style="width: 999px;"><SPAN class="lia-inline-image-display-wrapper" image-alt="jthi_0-1762581573217.png" style="width: 999px;"><SPAN class="lia-inline-image-display-wrapper" image-alt="jthi_0-1762581573217.png" style="width: 999px;"><span class="lia-inline-image-display-wrapper" image-alt="jthi_0-1762581573217.png" style="width: 999px;"><img src="https://us.v-cdn.net/6038925/uploads/images/86635iE72CED4B36533B07/jthi_0-1762581573217.png" role="button" title="jthi_0-1762581573217.png" alt="jthi_0-1762581573217.png" /></span></SPAN></SPAN></SPAN></SPAN></SPAN></SPAN></P>
<PRE><CODE class=" language-jsl">Names Default To Here(1);
dt = Open("$SAMPLE_DATA/Reliability/Bladder Cancer.jmp");
ra = dt << Recurrence Analysis(
Y(:Age),
Cost(:Cost),
Grouping(:Treatment Group),
Label(:Patient Number),
Fit Model(
Scale Effects(:Treatment Group),
Shape Effects(:Treatment Group),
Run Model,
Model Type("Proportional Intensity Poisson Process")
)
);
fit = Report(ra)[Outline Box("Fitted Recurrence Model")] << Get Scriptable Object;
intensity_col = fit << Save Intensity Formula;
cumul_col = fit << Save Cumulative Formula;
回复：JSL编码你好，
谢谢你的回复。
我想预测每台机器90天后的累计值（本研究中的维修次数）。我已经算出来了。
我只需要使用“fit = Report(ra)[OutlineBox("Fitted Recurrence Model Cause= Caster Assay-Left)] << Get Scriptable Object;”。这样，“Cause ="”就会生成我想要的因果级别的公式。
为了实现对 13 个原因级别的自动化处理，我创建了一个原因列表，并使用 for 循环遍历每个级别。
谢谢。
回复：JSL编码您使用的是“按列”吗？循环可能并非必要，但这也不失为一个好办法，只需确保您不会意外覆盖任何列即可。
如果您使用的是“按列”功能，那么广播功能如果能满足您的需求，也可能有效。
Names Default To Here(1);
dt = Open("$SAMPLE_DATA/Reliability/Bladder Cancer.jmp");
ra = dt << Recurrence Analysis(
Y(:Age),
Cost(:Cost),
Grouping(:Treatment Group),
Label(:Patient Number),
Fit Model(
Scale Effects(:Treatment Group),
Shape Effects(:Treatment Group),
Run Model,
Model Type("Proportional Intensity Poisson Process")
),
By(:Cause of Death)
);
fit = Report(ra[1])[Outline Box("Fitted Recurrence Model")] << Get Scriptable Object;
fit << Broadcast(Save Intensity Formula);
fit << Broadcast(Save Cumulative Formula);
