Fórmula acumulativa del análisis de recurrencia (codificación JSL)Hola a todos
Quiero guardar la fórmula acumulativa de un modelo para cada nivel de causa en un modelo de análisis de recurrencia. A continuación se muestra el código que tengo.
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>Automatización y scripting</LINGO-LABEL>
</LINGO-LABS>
<LINGO-SUB id="lingo-sub-912154" slang="en-US" mode="CREATE" translate="no">Re: Codificación JSL</LINGO-SUB>
<LINGO-BODY id="lingo-body-912154" slang="en-US" mode="CREATE">
<P>¿Qué intentas extraer o hacer al final? Tu script creará las fórmulas por ti, que parecen tener un aspecto similar a este.</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;
Re: Codificación JSLHola,
Gracias por su respuesta.
Quiero predecir los valores acumulados (número de reparaciones en este estudio) 90 días después para cada máquina. Ya lo he resuelto.
Solo necesito usar "fit = Report(ra)[OutlineBox("Fitted Recurrence Model Cause= Caster Assay-Left)] << Get Scriptable Object;". Así, "Cause =" generará la fórmula para el nivel de causa que deseo.
Para automatizar esto para 13 niveles de causa, simplemente creé una lista de causas y utilicé un bucle for para recorrer cada nivel.
Gracias.
Re: Codificación JSL¿Estás usando la opción "Por columna"? Probablemente no sea necesario usar un bucle, pero es una buena opción. Solo asegúrate de no sobrescribir ninguna columna por accidente.
Si utiliza la opción "Por columna", la difusión también podría funcionar si cumple con sus expectativas.
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);
