similar to: Cook's Distance in GLM (PR#9316)

Bug PR#9316 noted an inconsistency between the Cook's distance contours on plot.lm(x, which = 5) and the values given by cooks.distance(x) -- as shown in plot.lm(x, which = 4) -- for glms: http://bugs.r-project.org/cgi-bin/R/Analyses-fixed?id=9316;user=guest;selectid=9316 The suggested fix was to modify the contour levels by a dispersion factor, implemented as follows: dispersion <-

In previous versions of the gnm package, the terms component of "gnm" objects had a "classID" attribute. This caused problems when used with str as the following simple example illustrates: > x <- 1 > attr(x, "classID") <- "type1" > str(x) Class 'type1' Class 'type1' Class 'type1' Class 'type1' Class

Hi, I was wondering if I could use the matprod function from array.c in my own C routine. I tried the following as a test /* my_matprod.c */ # include <Rinternals.h> /* for REAL, SEXP etc */ # include <R_ext/Applic.h> /* array.c says need for dgemm */ /* following copied from array.c */ static void matprod(double *x, int nrx, int ncx, double *y, int nry, int ncy, double *z)

Hi, I am trying to speed up part of an algorithm in which certain columns of a large matrix (X) are replaced by the element-wise product of a matrix (M) and a vector (v). In R, the code might be X[, ind] <- M * v I have written a small C routine to do this for me, but the timing depends on how I define the storage.mode of the objects in R and the data types in C, in a way which I don't

Last week Giovanni Parrinello posted a message asking why various packages were loaded when he loaded an .Rdata file. Brian Ripley replied saying he thought it was because the saved workspace contained a reference to the namespace of ipred. (Correspondence copied below). This begs the question: how did the reference to the namespace of ipred come to be in the .Rdata file? Brian did say it is

Last week Giovanni Parrinello posted a message asking why various packages were loaded when he loaded an .Rdata file. Brian Ripley replied saying he thought it was because the saved workspace contained a reference to the namespace of ipred. (Correspondence copied below). This begs the question: how did the reference to the namespace of ipred come to be in the .Rdata file? Brian did say it is

This bug is not quite fixed - the example from my original report now = works using R-2.2.1, but plot(Uniform, 6) does not. The bug is due to if (show[6]) { ymx <- max(cook, na.rm =3D TRUE) * 1.025 g <- hatval/(1 - hatval) # Potential division by zero here # plot(g, cook, xlim =3D c(0, max(g)), ylim =3D c(0, ymx),=20 main =3D main, xlab =3D

I have a puzzle as to how R is computing Cook's distance in weighted linear regression. In this case cook's distance should be given not as in OLS case by h_ii*r_i^2/(1-hii)^2 divided by k*s^2 (1) (where r is plain unadjusted residual, k is number of parameters in model, etc. ) but rather by w_ii*h_ii*r_i^2/(1-hii)^2 divided by k*s^2,

-----Original Message----- From: Robert Gentleman [mailto:rgentlem@fhcrc.org] Sent: Wed 10/12/2008 18:47 To: Legaie, Roxane Subject: Re: Kegg.db with GOstats Hi Roxanne, Can you redirect your question to the mailing list. And, you can find the answer in the mailing list archives... best wishes Robert Legaie, Roxane wrote: > Dear Robert Gentleman, > I am currently working on

Is there a reliable way to link to a vignette or other overview document, which has been included in the "inst/doc" directory of a source package, within a package help file in Rd format? or maybe link to an index of such documents? David Professor David Firth Dept of Statistics University of Warwick Coventry CV4 7AL United Kingdom Voice: +44 (0)247 657 2581 Fax: +44 (0)247 652

On 29/01/2023 18:24, Orion Poplawski wrote: > It seems that thunderbird-102.7.1-1.el8 (at least on CentOS Stream) > broke OAuth authentication with outlook.office365.com.? Downgrading to > 102.4.0-1.el8 resolved the issue. > > Error console reports: > > XHR POST https://login.microsoftonline.com/common/oauth2/v2.0/token > [HTTP/1.1 400 Bad Request 293ms] > > Is

Hello everybody, this might be a trivial question, but I have been unable to find this using Google. I am trying to replace double backslashes with single backslashes using gsub. There seems to be some unexpected behaviour with regards to the replacement string "\\". The following example uses the string C:\\ which should be converted to C:\ . > gsub("\\\\",

No, blas/veclib is tested, so aprt this extreme case you should report some other more commonly used cases in which something fails on OS X. This will help us to work it out. As said, I'll try some tests without using veclib and let you know. I've fowarded this mail to r-devel, which seems to be the right place, so for future msg on the subject please use r-devel. stefano On

Hi all, I have a problem with cpu affinity in my R-2.11.1 installation compiled against ATLAS running on a Linux (Ubuntu 10.04) cluster under GridEngine. I wish to use Grid Engine's core binding feature to bind user processes into the number of cores they request on the cluster, thus preventing badly behaved multi-threaded libraries from consuming more cores than requested. An example of

Hi Community, I'm trying to reconcile Cook's Distances computed in glm. The following snippet of code shows that the Cook's Distances contours on the plot of Residuals v Leverage do not seem to be the same as the values produced by cooks.distance() or in the Cook's Distance against observation number plot. counts <- c(18,17,15,20,10,20,25,13,12) outcome <- gl(3,1,9)

I would like to follow up on another one of the documentation issues raised in the discussion on function hints. Duncan mentioned that the R core were working on preprocessing directives for .Rd files, which could possibly include some sort of include directive. I was wondering if a "includeexamples" directive might also be considered. It often makes sense to use the same example to

Hi all, I'm not familiar with R programming and I'm trying to reproduce a result from a paper. Basically, I have a dataset which I would like to model in terms of successive increments, i.e. (y denote empirical values of y) y_1 = y1, y_2 = y1 + delta1, y_3 = y1 + delta1 + delta2. ... y_m = y1 + sum_2^m delta j where delta_j donote successive increments in the y-values, i.e. delta

Full_Name: Heather Turner Version: 2.2.0 OS: Windows XP Submission from: (NULL) (137.205.240.44) Standardized residuals as calculated by rstandard.lm, rstandard.glm and plot.lm are Inf/NaN rather than zero when the un-standardized residuals are zero. This causes plot.lm to break when calculating 'ylim' for any of the plots of standardized residuals. Example:

Hi, Searching this alias I can find a number of guides and scripts that describe the configuration of Solaris to boot from a ZFS rootpool. However, these guides appear to be Solaris 10 x86 specific. Is the creation of a ZFS boot disk available for Solaris SPARC ? If so, could you point me in the direction of where I can obtain details of this new feature from. Thanks and Regards, Paul. PS:

I have a server that has one interface onto a private subnet, and another connected to our LAN. I found that the IP address in WINS was the IP address of the private interface, so things weren't working. By using the "interfaces" setting in smb.conf I was able to force the use of the correct IP address. Now, I notice that the manual says, "By default Samba will query the