Comment calculer +/- un écart type autour de ma droite de régression linéaire ?Bonjour, j'essaie d'utiliser JMP pour calculer l'écart type autour de ma droite 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 souhaite ajouter + et - un écart type au résultat de la régression pour ce point. Je sais calculer les intervalles de confiance, mais pas l'écart type. Une fois ce calcul effectué, j'aimerais l'intégrer à un script pour automatiser le calcul des résultats et les représenter graphiquement en fonction du temps. Merci d'avance !
Exploration et visualisation des donnéesPartage et communication des résultatsRe : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Merci beaucoup txnelson.Ce jsl m'aidera une fois que j'aurai établi les statistiques sous-jacentes.Je vous en suis reconnaissant!
Re : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Merci Marc.J'apprécie votre aide et vos idées !
Re : comment calculer +/- un écart type autour de ma ligne de régression linéaire ?Je pense que le tracé USGS pourrait être une moyenne d'échantillon et un écart type, et non une moyenne et un écart type prédits à partir d'un modèle.
Tout modèle doit supposer quelque chose sur la variance de la réponse. Par exemple, les modèles de régression linéaire supposent qu'il est constant (indépendant des conditions). De nombreux modèles courants font cette hypothèse. Vous avez déclaré que vous utilisiez la régression linéaire, j'ai donc supposé qu'elle est constante et estimée par le RMSE.
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 tracé ?(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%20CRRNA .
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 ?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 ?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 ?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 ) )
)
);
