[R] Influence.measures

Peter B. Mandeville mandevip at uaslp.mx
Wed Jun 23 11:00:19 CEST 1999


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 at deimos.tc.uaslp.mx
Jefe del Depto. de Informática y Bioestadística rpe1531 at 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 at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list