come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Ciao, sto cercando di ottenere jmp per calcolare una deviazione standard attorno alla mia linea di regressione lineare che sto usando come modello predittivo.In particolare, per ogni nuovo (futuro) punto di campionamento, vorrei aggiungere + e meno una deviazione standard al risultato della regressione per quel punto.So come calcolare gli intervalli di confidenza ma non una deviazione standard.Una volta che vedo come eseguire questo calcolo, vorrei inserirlo in uno script per automatizzare i risultati e la trama rispetto al tempo.Grazie in anticipo!
AnaliticaSceneggiaturaRe: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Ecco un esempio di un modo per farlo
names default to here(1);
// Open Data Table: big class.jmp
// → Data Table( "big class" )
dt = Open( "/C:/Program Files/SAS/JMPPRO/14/Samples/Data/big class.jmp" );
biv = dt << Bivariate( Y( :Weight ), X( :Height ), invisible );
biv << fit line( {Save Predicteds} );
biv << close window;
// New column: Column 7
Data Table( "big class" ) << New Column( "Upper",
Numeric,
"Continuous",
Format( "Best", 12 ),
Formula(
:Predicted weight + Col Std Dev( :weight )
));
// New column: Column 8
Data Table( "big class" ) << New Column( "lower",
Numeric,
"Continuous",
Format( "Best", 12 ),Formula(
:Predicted weight - Col Std Dev( :weight )
)
);
Graph Builder(
Size( 528, 454 ),
Show Control Panel( 0 ),
Variables(
X( :height ),
Y( :Predicted weight ),
Y( :Upper, Position( 1 ) ),
Y( :lower, Position( 1 ) )
),
Elements(
Points( X, Y( 1 ), Y( 2 ), Y( 3 ), Legend( 12 ) ),
Smoother( X, Y( 1 ), Y( 2 ), Y( 3 ), Legend( 13 ) )
)
);
Re: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Grazie txnelson.Ho notato che il tuo esempio mostra grafici rettilinei.Lo script funzionerebbe per i valori futuri previsti dalla mia formula di regressione?La mia formula di regressione prevede il valore di ecoli (dipendente) in base alla torbidità (indipendente) e sto cercando di creare un modello in cui utilizzo i futuri valori di torbidità in tempo reale per prevedere l'ecoli..... e mostrare anche +/ uno standard deviazione.La trama qui sotto è il mio obiettivo.
Una cosa che non mi è chiara è la parte "col std dev( :weight)" del tuo esempio.Ho bisogno di pensare a cosa viene calcolato std dev.Nel mio esempio, non sarebbe torbidità, né ecoli.Sarebbe invece una deviazione standard dell'adattamento di regressione, corretto?(Suppongo che sarebbe lo std dev della pendenza?) Apprezzo qualsiasi consiglio!
Re: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Le linee disegnate sono derivate dai valori previsti, quindi seguiranno qualsiasi lime di regressione che rientri. La scelta di quale colonna utilizzare nella probabilità col std dev sia la colonna prevista
Re: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Non dovresti usare la funzione Col Std Dev() per calcolare la deviazione standard sulla linea. Questa funzione restituisce la SD stimata incondizionata della variazione di risposta. Penso che tu voglia la SD condizionale. Questa quantità viene riportata come RMSE nella tabella Riepilogo dell'adattamento nel report di regressione. È condizionato dal predittore.
Re: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Grazie Marco.Nel mio diagramma di esempio (screenshot) sopra, l'intervallo di deviazione standard +/- uno fluttua di dimensioni intorno al valore di ecoli previsto (cioè regredito).Tuttavia, l'RMSE è a fisso valore basato sull'adattamento di regressione.Ovviamente c'è qualcosa che mi sfugge.In quali condizioni una SD potrebbe fluttuare nelle dimensioni nel tempo come nel grafico di esempio?(La trama proviene da una pagina Web USGS, quindi è stata senza dubbio completamente controllata.) https://www2.usgs.gov/water/southatlantic/ga/bacteria/#:~:text=coli)%20concentrations.-,E.,three%20stations%20within%20the%20CRNRA .
