Hi:

A standardized residual in (non-mixed) linear models has the form

                                      e_i / (s * sqrt(1 - h_i)),

where s is the (residual) standard deviation and h_i is the i-th element of
the 'hat matrix'. h_i is also called the leverage that the i-th observation
exerts on the model. You could look into the influence.ME package on CRAN to
see if something like it has already been implemented, but if not, all you
need is the hat matrix analogue of a mixed model and you can write it up
yourself. You can get the residuals and residual standard deviation from a
lme(r) object directly, and should be able to access the hat matrix, too,
perhaps indirectly.

HTH,
Dennis

On Tue, Nov 9, 2010 at 9:05 AM, Michael Larkin <mlarkin@rsmas.miami.edu>wrote:

> Anyone use the standardized residuals function in R?  It is called stdres
>
> I can't get it to work.  It tells me there is no documentation for stdres.
> I am assuming there is a package that I need to load.  I looked at the list
> of packages and I can't find one for stdres.  It must be bundled in some
> package.
>
> Any advice on how I can get stdres to work would be greatly appreciated.
>
> Mike
>
> _______________________________________________
> R-sig-mixed-models@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>

	[[alternative HTML version deleted]]


