similar to: suggested changes to summary.glm and summary.lm (PR#2776)

Dear list members, I was reminded of what I see as a problem with summary.glm this week when some of my students fell into a trap in a homework exercise, defining a logit model in which coefficients were aliased. When this happens in lm, summary.lm prints a message ["Coefficients: (1 not defined because of singularities)"], but summary.glm is silent. In both instances, the print

Hi! In R 3.0.0 from current SVN, ?summary.lm says: > Value [...] > df degrees of freedom, a 3-vector (p, n-p, p*), the last > being the number of non-aliased coefficients. ?summary.glm says: > df a 3-vector of the rank of the model and the number of residual > degrees of freedom, plus number of non-aliased coefficients. It seems to me that the description is reversed: p is

Dear R and S list members: I've just published a book entitled An R and S-PLUS Companion to Applied Regression. The book introduces R and S-PLUS in the context of linear and generalized-linear models, and assumes familiarity with -- or a simultaneous course or text on -- these topics. I thought that the book might be of interest to the members of these lists or to their students. More

I've uploaded to CRAN a new package called effects. The package contains functions for tabular or graphical display of terms in a linear or generalized linear model, and is particularly suitable for models that contain terms -- such as main effects and interactions, or polynomial regressors -- related by marginality (hierarchy). A draft paper describing the package is located at

Dear R-announce list members, I've submitted a new package to CRAN, Rcmdr, which provides a cross-platform basic-statistics GUI for R. Rcmdr is based on the tcltk package and uses only standard Tk widgets. It therefore runs under Windows if Tcl/Tk support is installed with R, although it requires the single-document interface (SDI). I've tested Rcmdr with Windows and Red Hat Linux;

Dear R-announce list members, I've submitted a new package to CRAN, Rcmdr, which provides a cross-platform basic-statistics GUI for R. Rcmdr is based on the tcltk package and uses only standard Tk widgets. It therefore runs under Windows if Tcl/Tk support is installed with R, although it requires the single-document interface (SDI). I've tested Rcmdr with Windows and Red Hat Linux;

I've uploaded to CRAN a new package called effects. The package contains functions for tabular or graphical display of terms in a linear or generalized linear model, and is particularly suitable for models that contain terms -- such as main effects and interactions, or polynomial regressors -- related by marginality (hierarchy). A draft paper describing the package is located at

Dear R-help list members, I'm encountering problems using expand.model.frame. To illustrate, consider the function > test <- function(model){ + expand.model.frame(model, "income") + } > The data frame Prestige (from the car library) has several variables, including prestige, income and education. I've attached this data frame >

Dear list members, I wonder whether anyone has succeeded in building either the RPgSQL or the Rdbi.PGSQL and Rdbi packages for Windows. If so, would you be willing to share either the binary package(s) or instructions about how to set up the configure scripts? My preliminary efforts to get things going under Windows 2000 haven't met with much success. Thanks, John

Dear list members, I'm experiencing a variety of intermittent problems with GUI-related commands (e.g., for dialog boxes and graphics-device windows) in R for Windows running under XEmacs/ESS. I'm using version 1.6.1 of R, version 21.4.9 of XEmacs, version 1.02 of ESS, and various versions of Windows, including 2000, 95, and 98SE. The problems don't occur in every circumstance,

Dear r-devel list members, I've prepared a version of the Rcmdr package that uses the new tools for internationalization and localization in R 2.1.x. For the present, this is the development version of the Rcmdr package (1.1-0), available at <http://socserv.socsci.mcmaster.ca/jfox/Misc/Rcmdr/index.html>. There is currently a French translation provided with the package, courtesy of

Dear R-help list members, Yesterday, I described a problem that I encountered with the use of expand.model.frame when the model object passed to this function was created without an explicit data argument in the call. I've had no responses, but I've been able to work around the problem by making a change to expand.model.frame, replacing the line rval <-

Dear R-announce list members, I've uploaded to CRAN a new version of my Rcmdr package, which provides a cross-platform basic-statistics GUI for R. This version has a few new features, but mostly addresses some problems that arose with the original version. The relevant portion of the CHANGES file for the package is appended. As before, suggestions and reports of problems are appreciated.

Dear list members, I've just uploaded version 1.0-0 of the Rcmdr package to CRAN. For people who haven't seen the package before, the "R Commander" provides a basic-statistics graphical user interface to R, based on the tcltk package. The new version incorporates a number of improvements to the R Commander interface, as documented in the CHANGES file distributed with the

Dear list members, I've just uploaded version 1.0-0 of the Rcmdr package to CRAN. For people who haven't seen the package before, the "R Commander" provides a basic-statistics graphical user interface to R, based on the tcltk package. The new version incorporates a number of improvements to the R Commander interface, as documented in the CHANGES file distributed with the

Dear list members, I and my students have encountered an intermittent problem using the Rcmdr package (version 0.9-0) under Windows (with the SDI). The problem occurs both in versions 1.7.1 and 1.8.0 of R. The problem seems to occur only in the following context: The "Edit data set" button is pressed in the Rcmdr GUI. This executes the fix() function on the active data frame,

Dear Steve, Since this relates as well to the message I posted a couple of minutes before yours, I agree that it?s possible to phrase ?best practices? too categorically. In the current case, I believe that it?s reasonable to say that specifying the data argument is ?generally? or ?usually? the best option. That doesn?t rule out exceptions. Best, John

>>>>> Fox, John <jfox at mcmaster.ca> >>>>> on Thu, 14 Sep 2017 13:46:44 +0000 writes: > Dear Martin, I made three points which likely got lost > because of the way I presented them: > (1) Singularity is an unusual situation and should be made > more prominent. It typically reflects a problem with the > data or the

Dear list members, In re-running with GLMM() from the lme4 package a generalized-linear mixed model that I had previously fit with glmmPQL() from MASS, I'm getting a warning of a convergence failure, even when I set the method argument of GLMM() to "PQL": > bang.mod.1 <- glmmPQL(contraception ~ as.factor(children) + cage + urban, + random=~as.factor(children) + cage +

Dear Martin, I think that everyone agrees that it?s generally preferable to use the data argument to lm() and I have nothing significant to add to the substance of the discussion, but I think that it?s a mistake not to add to the current examples, for the following reasons: (1) Relegating examples using the data argument to ?see also? doesn?t suggest that using the argument is a best practice.