Calcuation of STND residuals
[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Calcuation of STND residuals



Hi,
Could someone please tell me how the STND RES are calculated?
I can not get these residuals from the .yht file.
The attached is a simple example with a fixed effects only model.

I ran the analysis with and without the data point - S^2 is the variance
with it and S(i)^2 is the variance without.  Res and Hat are the
residuals and Hat value from the .yht file.

As per pg 376-377 of McCullagh and Nelder
Stand Res was calculated as Res /(sqrt(1-Hat))
Student Res was calculated as Stand Res/S
Del Stud Res was calculated as Student Res * S / S(i)

S^2         S(i)^2    Res     Hat      Stand Res        Student Res   
Del Stud Res
1.55357	0.758523	4	0.1295	4.287218357	3.439617822	7.044858309

The value from the .res file is 3.32.

Also, in at least this version of ASREML, potential outliers are those
with STND RES of >3.3 rather than the 3.5 said on page 115 of the March
21, 2000 manual.

Thanks 
Bruce
  Record        Yhat    Residual         Hat
       1      1.3750     -0.3750      0.1942    
       2      1.3750     -0.3750      0.1942    
       3      1.3750      0.6250      0.1942    
       4      1.3750      0.6250      0.1942    
       5      2.0000       0.000      0.1295    
       6      2.0000       0.000      0.1295    
       7      2.0000       1.000      0.1295    
       8      2.0000      -1.000      0.1295    
       9      2.0000       4.000      0.1295    
      10      1.3750     -0.3750      0.1942    
      11      1.3750     -0.3750      0.1942    
      12      1.3750     -0.3750      0.1942    
      13      2.0000      -1.000      0.1295    
      14      2.0000      -1.000      0.1295    
      15      2.0000       1.000      0.1295    
      16      1.7500       2.250      0.3884    
      17      1.7500     -0.7500      0.3884    
      18      1.7500     -0.7500      0.3884    
      19      1.3750      0.6250      0.1942    
      20      2.0000      -1.000      0.1295    
      21      2.0000      -1.000      0.1295    
      22      2.0000       0.000      0.1295    
      23      2.0000      -1.000      0.1295    
      24      1.7500     -0.7500      0.3884    

fig.mod

 SLOPES FOR LOG(ABS(RES)) on LOG(PV) for Section   1
   2.56
                  
 *                
 * *              
 * *              
 ***  *           
 ***  **          
 ***  **   *     *
   
 Min Mean Max -1.0000      0.0000      4.0000     omitting     3 zeros

 STND RES    9  3.32
 5 3
   0.000000       0.000000       1.000000    
RSTRUCTURE               1   1
 VARIANCE                1   1   0 1.00000    
 STRUCTURE              24   0   0
  parity                          1           0.000       0.000    
  parity                          2          0.6250      0.5689    
  parity                          3          0.3750      0.7633    
  mu                              1           1.375      0.4407    
   0.155357E+01   0.000000E+00   0.000000E+00   0.625000E+00   0.000000E+00
   0.323661E+00   0.375000E+00   0.000000E+00   0.194196E+00   0.582589E+00
   0.137500E+01   0.000000E+00  -0.194196E+00  -0.194196E+00   0.194196E+00
 Variance of Variance components     1
  0.229865    

fig.as