Message-Id: <3.0.5.16.19990623110019.37d7a712@deimos.tc.uaslp.mx>
Date: Wed, 23 Jun 1999 11:00:19
To: r-help@stat.math.ethz.ch
From: "Peter B. Mandeville" <mandevip@uaslp.mx>
Subject: [R] Influence.measures
I am using rw0641 with Windows 98. To list just the influential
repetitiones that result from "influence.measures", I am using the input
result <- lm(y~x)
and the code from the example in the help for "influence.measures"
INFLM <- function(result){
inflm <- influence.measures(result)
which(apply(inflm$is.inf,1,any))
}
It works fine up to now with the exception of Huber's data published in
Atkinson 1985:8 (Plots, Transformations and Regression)
Observation x y
1 -4 2.48
2 -3 0.73
3 -2 -0.04
4 -1 -1.44
5 0 -1.32
6 10 0.00
which is not very well excepted. The command
> inflm <- influence.measures(result)
reports
Warning: NaNs produced in function "qf"
and
> which(apply(inflm$is.inf,1,any))
reports
Error: names attribute must be the same length as the vector
What should I do?
Thank you very much.
Peter B.
-- Peter B. Mandeville mandevip@deimos.tc.uaslp.mx Jefe del Depto. de Informática y Bioestadística rpe1531@pasteur.fmed.uaslp.mx Facultad de Medicine Tel: 48 26-23-45 ext. 232 Universidad Autónoma de San Luis Potosí Fax: 48 28-23-52 Av. V. Carranza 2405 Col. Los Filtros Apartado Postal 145 San Luis Potosí, S.L.P. 78210 México-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request@stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._