224723_fr-FR

liliconcom

Résidus studentisés en externe

Salut à tous! 

 

J'ai un client qui souhaite accéder aux résidus étudiés en externe utilisés pour le graphique des résidus étudiés dans la régression par la méthode des moindres carrés.  Quelqu'un a-t-il un moyen de calculer cela en dehors de la plateforme ou de savoir comment atteindre ces valeurs?

 

Meilleur,

 

M

Re: Résidus externe de l'étudiant

Il me semble que les colonnes Enregistrer> Résidus studentisés sont les valeurs tracées sur ce graphique. 

Ledi

Re: Résidus externe de l'étudiantJe le pensais aussi, mais dans la doc, ils sont différents ( https://www.jmp.com/support/help/14-2/row-diagnostics.shtml ).

M
Re: Résidus externe de l'étudiant

Je n'ai pas regardé les chiffres, mais vous pouvez toujours obtenir les valeurs directement à partir de la position du marqueur dans le rapport:

NamesDefaultToHere(1);

dt = Open("$SAMPLE_DATA/Big Class.jmp");

// Fit Model
fm = dt << Fit Model(
     Y( :height ),
     Effects( :weight ),
     Personality( "Standard Least Squares" ),
     Emphasis( "Effect Leverage" ),
     Run(
      :height << {Summary of Fit( 0 ), Analysis of Variance( 0 ),
      Parameter Estimates( 0 ), Effect Tests( 0 ), Effect Details( 0 ),
      Lack of Fit( 0 ), Scaled Estimates( 0 ), Plot Actual by Predicted( 1 ),
      Plot Residual by Predicted( 1 ), Plot Studentized Residuals( 1 ),
      Plot Effect Leverage( 1 ), Plot Residual by Normal Quantiles( 0 ),
      Box Cox Y Transformation( 0 )}
     ),
     SendToReport(
      Dispatch(
       {"Response height", "Whole Model"},
       "Regression Plot",
       OutlineBox,
       {Close( 1 )}
      ),
      Dispatch(
       {"Response height", "Whole Model"},
       "Actual by Predicted Plot",
       OutlineBox,
       {Close( 1 )}
      )
     )
    );

// Get values directly from the marker seg . . .    
fmRep = Report(fm);
frame = fmRep[FrameBox(4)];
mySeg = (frame << Find Seg( Marker Seg( 1 ) ));
xVals = mySeg << Get X Values;
yVals = mySeg << Get Y Values;
dt2 = asTable(xVals||yVals, << columnNames({"Row Number", "Studentized Residual"}));
dt2 << setName("Residuals");

Re: Résidus externe de l'étudiant

C'est intéressant. J'ai essayé de cliquer avec le bouton droit de la souris et de sélectionner Personnaliser dans l'intrigue pour voir ce qui pourrait être disponible. Il n'y a rien que je puisse utiliser, donc je n'ai pas pensé aux segments.

Re: Résidus externe de l'étudiant

FWIW, cette méthode ne m’a jamais manqué (pour le moment. . . ).

Re: Résidus externe de l'étudiant

Impressionnant! Merci, @ian_jmp !

 

Cela les a!  

 

De plus, j'étais perplexe au sujet de cette nuit et je suis tombé sur une équation qui convertira les résidus étudiés en résidus étudiés à l'extérieur.  Vous devez simplement fournir le nombre de facteurs dans le modèle.  L'équation est extraite de la régression linéaire appliquée de Weisburg. c'est à la page 216 de la 4ème édition.  J'ai ajouté un script à la fin de votre script pour que tout soit au même endroit. 

 

Meilleur, 

 

M

</CODE></P> <PRE><CODE class=" language-jsl">Names Default To Here( 1 ); Clear Log(); dt = Open( "$SAMPLE_DATA/Big Class.jmp" ); // Fit Model fm = dt << Fit Model( Y( :height ), Effects( :weight ), Personality( "Standard Least Squares" ), Emphasis( "Effect Leverage" ), Run( :height << {Summary of Fit( 0 ), Analysis of Variance( 0 ), Parameter Estimates( 0 ), Effect Tests( 0 ), Effect Details( 0 ), Lack of Fit( 0 ), Scaled Estimates( 0 ), Plot Actual by Predicted( 1 ), Plot Residual by Predicted( 1 ), Plot Studentized Residuals( 1 ), Plot Effect Leverage( 1 ), Plot Residual by Normal Quantiles( 0 ), Box Cox Y Transformation( 0 )} ), SendToReport( Dispatch( {"Response height", "Whole Model"}, "Regression Plot", OutlineBox, {Close( 1 )} ), Dispatch( {"Response height", "Whole Model"}, "Actual by Predicted Plot", OutlineBox, {Close( 1 )} ) ) ); /*METHOD 1 Pull Externally Studentized Residuals from the plot */ // Grab the report object fmRep = Report( fm ); // Grab the ESR plot frame box frame = fmRep[FrameBox( 4 )]; // Grab the marker segment in the frame box and pull the values mySeg = (frame << Find Seg( Marker Seg( 1 ) )); yVals = mySeg << Get Y Values; /*METHOD 2 Convert Studentized Residuals to Externally Studentized Residuals (Applied Linear Regression, 4th Ed, 2014 - pg 216) note the denominator removes one df for the intercept so p should be the df from the model ANOVA*/ // define the total number of terms in the model from the model ANOVA p = 1; //export the studentized residuals to the data table and gather table information fm << Studentized Residuals; r = Column( N Cols() ); n = N Rows( dt ); // create a blank matrix for the values ESR = []; // loop through the residuals matrix (last column in the dt) and calculate the ESR. For( i = 1, i <= n, i++, // calculation per reference above k = r[i] * Sqrt( (n - p - 2) / (n - p - 1 - (r[i] ^ 2)) ); // insert the ESR value into the matrix ESR = V Concat( ESR, k ); ); // create a new column in the data table with the values (from both methods and method 2 as a formula column) dt << New Column( "Externally Studentized Residuals (From Report)", Set Values( yVals ) ); dt << New Column( "Externally Studentized Residuals", Set Values( ESR ) ); dt << New Column( "Externally Studentized Residuals (Formula)", Formula( Parameter( {p = 1}, :Studentized Resid height * Sqrt( (N Rows() - p - 2) / (N Rows() - p - 1 - :Studentized Resid height ^ 2) ) ) ) );