similar to: rigde regression and influence measures

Dear list, There is, I believe, a minor typo in the example section of influence.measures.Rd. In the final example the word `does` appears where I suspect `dose` is required: I couldn't remember exactly what format patches should be in, so here is one as diff would produce: Index: devel/src/library/stats/man/influence.measures.Rd

Hi dear all, I'm wondering about the question that; Does the influence.measures(model) for linear models valid for general linear models such as logistic regression models? That is; If I fit the model like model <- glm( y~X1+X2, family=binomial) Then, if i apply the function "influence.measures(model), i will get the result of influence measures. These result are valid for

influence.measures gives several measures of influence for each observation (Cook's Distance, etc) and actually flags observations that it determines are influential by any of the measures. Looks good! But how does it discriminate between the influential and non- influential observations by each of the measures? Like does it do a Bonferroni-corrected t on the residuals identified by

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

A while back I asked about getting a list of points that R considers influential after fitting a linear model, and very quickly got a helpful pointer to influence.measures(). But "it has happened again." The trouble I am having is that points marked on plots are not flagged in the output from influence.measures(), and I can't read them on the plots. I tried some successive

Dear list, I've noticed that influence.measures and cooks.distance gives different results for non-gaussian GLMs. For example, using R-1.9.0 alpha (2003-03-17) under Windows: > ## Dobson (1990) Page 93: Randomized Controlled Trial : > counts <- c(18,17,15,20,10,20,25,13,12) > outcome <- gl(3,1,9) > treatment <- gl(3,3) > glm.D93 <- glm(counts ~ outcome +

I'm following up on a question I posted 8/10/2010, but my newsreader has lost this thread. > Barrett & Ling, JASA, 1992, v.87(417), pp184-191 define general > classes of influence measures for multivariate > regression models, including analogs of Cook's D, Andrews & Pregibon > COVRATIO, etc. As in univariate > response models, these are based on leverage and

Many thanks to those who replied to my message regarding installation of the lattice package. I have now upgraded my version of R. I have a quick question: I am running mixed effects models using the lme command. I am wanting to assess the fit of these models using the usual plots including Cooks Distance plots. I have noticed from the help files that it is possible to run influence measures

Barrett & Ling, JASA, 1992, v.87(417), pp184-191 define general classes of influence measures for multivariate regression models, including analogs of Cook's D, Andrews & Pregibon COVRATIO, etc. As in univariate response models, these are based on leverage and residuals based on omitting one (or more) observations at a time and refitting, although, in the univariate case, the

I?ve a data set and i want fit a kernel density estimate to the data. but using the k-nearest neighbour method. How i do this with R. thanks -- bertola at fastmail.fm --

Is there a command in R that make the same regression like l1fit in S-plus? -- bertola at fastmail.fm --

Hi, I've used the Mathematica to produce 3D graphics, contour plots of a bivariate normal distribution Now I want make these graphics in R, but i do not know how. I would like to: - Plot a 3D graph for some different variance matrix - Plot the contour plots - Find and try to plot (in the 3d graph ou contour plot) the (1-a)% confidence region based in a chi-square(a) with the degrees of

You'll find that it is a lot easier to do it in R: # lets first simulate a bivariate normal sample library(MASS) bivn <- mvrnorm(1000, mu = c(0, 0), Sigma = matrix(c(1, .5, .5, 1), 2)) # now we do a kernel density estimate bivn.kde <- kde2d(bivn[,1], bivn[,2], n = 50) # now plot your results contour(bivn.kde) image(bivn.kde) persp(bivn.kde, phi = 45, theta = 30) # fancy contour with

I'm looking for data to use with nonparametrics exploration technics. I'm a undergraduate student in statistics at the Unicamp (Brazil), and i've to make a project (with orientation of a professor). I choose nonparametrics statistics to make the project. But in this project i must have to use a databank and present analisys. So, if somebody give me indications where i can find

I am a novice R user, and I am having difficulty understanding R's influence plots. I am trying to remove outliers from a particular variable, "sib." I am able to generate influence plots and further outlier information such as below (which is a shortened example). For my analyses, I end up excluding the points R refers to, 7, 18, 26, and 105. However, my question is, how can I

Hello: I am trying to understand the method 'hatvalues(...)', which returns something similar to the diagonals of the plain vanilla hat matrix [X(X'X)^(-1)X'], but not quite.  A Fortran programmer I am not, but tracing through the code it looks like perhaps some sort of correction based on the notion of 'leave-one-out' variance is being applied. Whatever the

Hi all, Does lme has function to compute the cook's distance or influence measure like lm? I can't find them. Thanks. Yen Lin [[alternative HTML version deleted]]

Hi all I'm getting a NaN returned on using dffits, as explained below. To me, there seems no obvious (or non-obvious reason for that matter) reason why a NaN appears. Before I start digging further, can anyone see why dffits might be failing? Is there a problem with the data? Consider: # Load data dep <-

I am hoping that some one might be to tell me whether there are any functions that produce influence measures for lme /nlme objects (i.e those suggested by Lesaffre and Verbeke or Langford and Lewis for multilevel models). Thanks in advance. ----------------------------------------------- Anthony A. Pezzola apezzola@uc.cl (02) 354-7823 Profesor de Ciencia Política Instituto de Ciencia

Greetings! Recently when I tried to use lm.influence I get the following error: Error in .Fortran("lminfl", model$qr$qr, n, n, k, as.integer(do.coef), : Fortran symbol name "lminfl" not in DLL for package "base" This occurs on both Linux and Windows platforms (details below). Searching the mail lists and other sources indicates that the fortran code for