similar to: influence.measures()

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

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 <-

R-users E-mail: r-help@r-project.org Hi! R-users. I am just wondering what the definition of "dffits" in R language is. Let me show you an simple example. function() { library(MASS) xx <- c(1,2,3,4,5) yy <- c(1,3,4,2,4) data1 <- data.frame(x=xx, y=yy) lm.out <- lm(y~., data=data1, x=T) lev1 <- lm.influence(lm.out)$hat sig1 <-

Hey guys, I've search a few threads about deleting a value from a vector, but no one has addressed this question so far. I want to delete the last value from a string of values I have: r = [ 1, 2, 3, 4, 5 ], and i want to make r2 = to [ 1, 2, 3, 4] So that r2 is just like r, except that it missing the final value. Thanks, -- View this message in context:

All - My question is a bit involved, so bear with me. I have some data that looks like: Lake LL LW 81 2.176091259 1.342422681 81 2.176091259 1.414973348 81 2.176091259 1.447158031 81 2.181843588 1.414973348 81 2.181843588 1.447158031 81 2.184691431 1.462397998 81 2.187520721 1.447158031 81 2.187520721 1.477121255 81 2.187520721 1.505149978 ... [truncated] I'm trying to: 1) fit a simple

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

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

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 use R on my Ubuntu 9.04 laptop, which was installed by "aptitude install" way Now i'm learning a book of R which needs "MASS" and "DAAG" installed. So i followed the instructions: >install.packages("MASS") >library("MASS") MASS works fine. But DAAG doesn't. Anyone who could help would be appreciated a lot ! Below are my error

I use R on my Ubuntu 9.04 laptop, which was installed by "aptitude install" way Now i'm learning a book of R which needs "MASS" and "DAAG" installed. So i followed the instructions: >install.packages("MASS") >library("MASS") MASS works fine. But DAAG doesn't. Anyone who could help would be appreciated a lot ! Below are my error

Hi Everyone, I am building a script for a client and my imagecreatefromjpeg() is throwing a fatal exception, and dying. It seems as my server with CentOS 4, PHP 4.3.9 doesn't have >= GD-1.8 support or better. I did a phpinfo() and I can not find any reference to GD..............Any help on this would be well appricated -- Best Regards, Jon McCauley

Hi everyone I'm interested in using the Telappliant/voip-talk offering as an alternative to my DDI analog problem. (see [Asterisk-Users] Multiple DDI & Hunting on Analog Lines (UK) for details) Does anyone on the list have any recent comments on reliability etc? I would really appricated some positive and negative comments. Cheers Matt

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]]

Dear all, When I print out princomp's loading outputs, there is alwasy a section for "SS loading", "Proportional Var" and "Cumulative Var". Anybody can tell what they are for? Or anyone can direct me to some reference to read about? Any help will be highly appricated. Hongsheng [[alternative HTML version deleted]]

Hey, I have rather new installation of xen, dom0 is running ntpd and is perfectly in time. All domU''s are ~23s in future. Running ntpd on the domU''s seem to have no effect, /proc/sys/xen/independent_wallclock appears to have no effect. I''m running xen-unstable and 2.6.16-rc4 in AMD64 box. Any tip would be appricated. Thanks, -- ++ytti

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

Hi there I have created a multivariate logistic regression model looking at the presence/absence of disease on farms. I would like to plot the diagnostic plots recommended by Hosmer & Lemeshow to look particularly for any points of high influence. In order to do this I need to extract values for delta-beta. The function dfbeta gives a value for change in each coefficient but I am looking

hi, I have strange problems with HTB and several hundred classes flat structure i.e. root |--50kbps |--30kbps |--50kbps |--80kbps |--100kbps .... several hundred classes like this Ceil is the same as rate. The machine get no more than 2-3% average cpu(2.4Ghz pentium). What happens is that from time to time the traffic got "stalled". I tried numerous things to solve

What I am trying to do is the following: A call is sent to the * box via a SIP invite. The * box answers via an IVR menu system with " Enter the extension you want to dial" so I enter in my 5 digit extension and get the below message. Jan 18 10:10:03 WARNING[-1380238416]: channel.c:1860 ast_request: No channel type registered for 'SIP)' Jan 18 10:10:03 NOTICE[-1380238416]:

Hello, I recently installed the ruby on rails framework from this tutorial: http://www.hivelogic.com/narrative/articles/ruby_rails_lighttpd_mysql_tiger So, I installed: Ruby on Rails 1.1 Ruby 1.8.4 LightTPD 1.4.13 FastCGI 2.4.0 RubyGems 0.9.0 Readline 5.1 PCRE 6.6 FastCGI and MySQL bindings I figured out how to keep my ruby on rails up to date by using the command "gem update rails