¿Cómo calculo +/- una desviación estándar alrededor de mi línea de regresión lineal?Hola, estoy tratando de hacer que jmp calcule una desviación estándar alrededor de mi línea de regresión lineal que estoy usando como modelo predictivo.Específicamente, para cada punto de muestra nuevo (futuro), me gustaría agregar + y menos una desviación estándar al resultado de la regresión para ese punto.Sé cómo calcular intervalos de confianza pero no una desviación estándar.Una vez que vea cómo hacer este cálculo, me gustaría colocarlo en un script para automatizar los resultados y trazar contra el tiempo.¡Gracias por adelantado!
Analíticasecuencias de comandosRe: ¿cómo calculo +/- una desviación estándar alrededor de mi línea de regresión lineal?Aquí hay un ejemplo de una manera de hacer esto
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: ¿cómo calculo +/- una desviación estándar alrededor de mi línea de regresión lineal?Gracias txnelson.Me doy cuenta de que su ejemplo muestra diagramas de línea recta.¿Funcionaría la secuencia de comandos para valores futuros predichos por mi fórmula de regresión?Mi fórmula de regresión predice el valor de ecoli (dependiente) en función de la turbidez (independiente), y estoy tratando de crear un modelo en el que uso futuros valores de turbidez en tiempo real para predecir ecoli... y también mostrar +/ un estándar desviación.La siguiente trama es mi objetivo.
Una cosa que no me queda clara es la parte "col std dev(:weight)" de su ejemplo.Necesito pensar qué desviación estándar se está calculando.En mi ejemplo, no sería turbidez, ni sería ecoli.En cambio, sería una desviación estándar del ajuste de regresión, ¿correcto?(¿Supongo que sería la desviación estándar de la pendiente?) ¡Agradezco cualquier consejo!
Re: ¿cómo calculo +/- una desviación estándar alrededor de mi línea de regresión lineal?Las líneas dibujadas se derivan de los valores predichos, por lo tanto, seguirán cualquier tiempo de regresión que ajuste. La elección de qué columna usar en la probabilidad de desarrollo estándar de columna sea la columna predicha
Re: ¿cómo calculo +/- una desviación estándar alrededor de mi línea de regresión lineal?No debe usar la función Col Std Dev() para calcular la desviación estándar sobre la línea. Esta función devuelve la SD estimada incondicional de la variación de la respuesta. Creo que quieres el SD condicional. Esta cantidad se informa como RMSE en la tabla Resumen de ajuste en el informe de regresión. Está condicionado por el predictor.
Re: ¿cómo calculo +/- una desviación estándar alrededor de mi línea de regresión lineal?gracias marcaEn mi gráfico de ejemplo (captura de pantalla) anterior, el rango de +/- una desviación estándar fluctúa en tamaño alrededor del valor ecoli predicho (es decir, retrocedido).Sin embargo, el RMSE es un reparado valor basado en el ajuste de regresión.Obviamente, hay algo que me estoy perdiendo.¿Bajo qué condiciones una SD fluctuaría en tamaño a través del tiempo como en la gráfica de ejemplo?(La trama es de una página web de USGS, por lo que sin duda ha sido examinada por completo). https://www2.usgs.gov/water/southatlantic/ga/bacteria/#:~:text=coli)%20concentrations.-,E.,tres%20estaciones%20dentro%20del%20CRNRA .
