come calcolo +/- una deviazione standard attorno alla mia linea di regressione lineare?Ciao, sto cercando di far calcolare a JMP una deviazione standard attorno alla mia retta di regressione lineare che utilizzo come modello predittivo. Nello specifico, per ogni nuovo punto campione (futuro), vorrei aggiungere + e - una deviazione standard al risultato della regressione per quel punto. So come calcolare gli intervalli di confidenza, ma non la deviazione standard. Una volta capito come eseguire questo calcolo, vorrei inserirlo in uno script per automatizzare i risultati e rappresentarli graficamente in funzione del tempo. Grazie in anticipo!
Esplorazione e visualizzazione dei datiCondivisione e comunicazione dei risultatiRe: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Grazie mille txnelson.Questo jsl aiuterà una volta che avrò elaborato le statistiche sottostanti.Lo apprezzo!
Re: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Grazie Marco.Apprezzo il tuo aiuto e le tue intuizioni!
Re: come faccio a calcolare +/- una deviazione standard attorno alla mia retta di regressione lineare?Penso che il grafico USGS potrebbe essere una media campionaria e una deviazione standard, non una media prevista e una deviazione standard da un modello.
Qualsiasi modello deve presumere qualcosa sulla varianza della risposta. Ad esempio, i modelli di regressione lineare presuppongono che sia costante (indipendente dalle condizioni). Molti modelli comuni fanno questa ipotesi. Hai affermato che stai usando la regressione lineare, quindi ho assunto che fosse costante e stimato dall'RMSE.
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 .
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?Le linee tracciate 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?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?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 ) )
)
);
