comment calculer +/- un écart type autour de ma droite de régression linéaire ?Salut - J'essaie de faire en sorte que jmp calcule un écart type autour de ma ligne de régression linéaire que j'utilise comme modèle prédictif.Plus précisément, pour chaque nouveau point d'échantillonnage (futur), je voudrais ajouter + et moins un écart type au résultat de la régression pour ce point.Je sais calculer des intervalles de confiance mais pas un écart type.Une fois que j'ai vu comment faire ce calcul, je voudrais le placer dans un script pour automatiser les résultats et tracer en fonction du temps.Merci d'avance!
AnalytiqueScriptRe : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Voici un exemple d'une façon de procéder
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 : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Merci txnelson.Je remarque que votre exemple montre des tracés linéaires.Le script fonctionnerait-il pour les valeurs futures prédites par ma formule de régression ?Ma formule de régression prédit la valeur d'ecoli (dépendante) en fonction de la turbidité (indépendante), et j'essaie de créer un modèle dans lequel j'utilise les futures valeurs de turbidité en temps réel pour prédire l'ecoli ..... et également afficher +/ une norme déviation.L'intrigue ci-dessous est mon objectif.
Une chose sur laquelle je ne suis pas clair est la partie "col std dev( :weight)" de votre exemple.J'ai besoin de réfléchir à ce que std dev est calculé.Dans mon exemple, ce ne serait pas la turbidité, ni l'ecoli.Au lieu de cela, ce serait un écart type de l'ajustement de régression, n'est-ce pas ?(Je suppose que ce serait le développement standard de la pente ?) J'apprécie tout conseil !
Re : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Les lignes tracées sont dérivées des valeurs prédites, elles suivront donc toute courbe de régression que vous ajusterez. Le choix de la colonne à utiliser dans la probabilité col std dev soit la colonne prédite
Re : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Vous ne devez pas utiliser la fonction Col Std Dev() pour calculer l'écart type autour de la droite. Cette fonction renvoie l'écart-type estimé inconditionnel de la variation de la réponse. Je pense que vous voulez le SD conditionnel. Cette quantité est rapportée en tant que RMSE dans le tableau Summary of Fit du rapport de régression. Il est conditionné par le prédicteur.
Re : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Merci Marc.Dans mon exemple de tracé (capture d'écran) ci-dessus, la plage d'écart type +/- un fluctue en taille autour de la valeur d'ecoli prédite (c'est-à-dire régressée).Cependant, le RMSE est un fixé valeur basée sur l'ajustement de régression.De toute évidence, il y a quelque chose qui me manque.Dans quelles conditions la taille d'un SD fluctuerait-elle dans le temps, comme dans l'exemple de graphique ?(L'intrigue provient d'une page Web de l'USGS, elle a donc sans aucun doute été entièrement vérifiée.) https://www2.usgs.gov/water/southatlantic/ga/bacteria/#:~:text=coli)%20concentrations.-,E.,three%20stations%20within%20the%20CRNRA .
