Necesito ayuda para editar un script desarrollado originalmente para extraer parámetros de crecimiento microbiano de un modelo Logistic 4P (Fit Curve Platform). Me gustaría editar el script para extraer los parámetros de cada tratamiento "cepa, medios, replicar", ver el resultado deseado para lo que me gustaría. Obviamente me falta algún lenguaje importante en el guión, pero no tengo claro cómo hacerlo. ¡Una respuesta de secuencias de comandos sería genial! ¡Por favor ayuda!
Este es el error que recibo:
"Problema de subíndice en el acceso o la evaluación de 'r [Cuadro de esquema (Estimaciones de parámetros"), Cuadro de número de columnas ("Estimación", 1)]', r] / * ### * / Cuadro de esquema ("Estimaciones de parámetros"), Número Col Box ("Estimación", 1)]
Guión que me gustaría editar: "Extraer valores críticos"
Script para la tabla de salida deseada - "Salida deseada"
Extraer valores críticos
Names Default To Here(1);
// Determine coefficients of Logistic 4P model, for estimation of growth parameters.
// This assumes that the current data table contains the data of interest, and that:
// the column containing the time values is labeled "Time (h)"
// the column containing the data values is labeled "Ln OD600"
// the column containing “strain” and “media” and “replicate” are grouping columns
//
// The names can be changed in the lines below,
// or the code can be changed to add a file dialog and a column selection,
// or the code can be changed to loop through multiple Y columns.
// First step: Fit Logistic 4P curve to data, store results in object Test.
Test = Fit Curve(
Y( :Ln OD600 ),
X( :Name( "Time (h)" ) ),
By( :Strain, :Media, :Replicate ),
Fit Logistic 4P
);
// 2nd step: Send report to variable r.
r = Test<<REPORT></REPORT>
Creo que hay un nivel innecesario de suscripción en el árbol de visualización que está causando el error. Aquí hay un ejemplo de lo que quiero decir:
a = GrowthRate = r[outlinebox("Parameter Estimates"),NumberColBox("Estimate",1)] << get(1);
Creo que eliminar el último nivel evitará el error:
a = GrowthRate = r[outlinebox("Parameter Estimates"),NumberColBox("Estimate")] << get(1);Hola @markbailey
Gracias por la ayuda, pero sigo recibiendo el mismo error. Eliminé el último nivel o el "1" después de "Estimación". ¡Cualquier otra sugerencia es bienvenida!
Aquí hay un ejemplo de cómo lo hago.
Names Default to Here( 1 );
dt = Open( "$SAMPLE_DATA/Big Class.jmp" );
fc = dt << Fit Curve(
Y( :weight ),
X( :height ),
Fit Exponential 2P
);
fcr = fc << Report;
{ a, b } = fcr["Parameter Estimates"][NumberColBox("Estimate")] << Get;