Ho bisogno di aiuto per modificare uno script originariamente sviluppato per estrarre i parametri di crescita microbica da un modello Logistic 4P (Fit Curve Platform). Vorrei modificare lo script per estrarre i parametri da ciascun trattamento "strain, media, replicate", vedere l'output desiderato per quello che vorrei. Mi manca ovviamente un linguaggio importante nella sceneggiatura, ma non sono chiaro come farlo. Una risposta di scripting sarebbe fantastica! Per favore aiuto!
Questo è l'errore che sto riscontrando:
"Problema di sottoscrizione nell'accesso o nella valutazione di 'r [Casella di struttura (stime dei parametri"), Casella di numeri ("Stima", 1)] ", r] / * ### * / Casella di struttura (" Stime dei parametri "), Casella numerica ("Stima", 1)]
Script che vorrei modificare - "Estrai valori critici"
Script per la tabella di output desiderata - "Output desiderato"
Estrai valori critici
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>
Penso che ci sia un livello non necessario di sottoscrizione nell'albero di visualizzazione che causa l'errore. Ecco un esempio di cosa intendo:
a = GrowthRate = r[outlinebox("Parameter Estimates"),NumberColBox("Estimate",1)] << get(1);
Penso che la rimozione dell'ultimo livello eviterà l'errore:
a = GrowthRate = r[outlinebox("Parameter Estimates"),NumberColBox("Estimate")] << get(1);Ciao @markbailey
Grazie per l'aiuto, ma continuo a ricevere lo stesso errore. Ho rimosso l'ultimo livello o "1" dopo "Stima". Qualsiasi altro suggerimento è il benvenuto!
Ecco un esempio di come lo faccio.
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;