similar to: R News, volume 8, issue 2 is now available

Dear useRs, I very much like the effect display of the proportional odds model on page 29 (Figure 8) of the following paper by John Fox: http://socserv.mcmaster.ca/jfox/Papers/logit-effect-displays.pdf It really gives a very concise overview of the model. I would like to use it to illustrate the proportional odds mixed models we fit here for a project on Diabetes but I can't seem to reproduce

I've prepared preliminary versions of two packages that I plan eventually to contribute to CRAN: car (for "Companion to Applied Regression") is a package that provides a variety of functions in support of linear and generalized linear models, including regression diagnostics (e.g., studentized residuals, hat-values, Cook's distances, dfbeta, dfbetas, added-variable plots,

I've prepared preliminary versions of two packages that I plan eventually to contribute to CRAN: car (for "Companion to Applied Regression") is a package that provides a variety of functions in support of linear and generalized linear models, including regression diagnostics (e.g., studentized residuals, hat-values, Cook's distances, dfbeta, dfbetas, added-variable plots,

I've prepared preliminary versions of two packages that I plan eventually to contribute to CRAN: car (for "Companion to Applied Regression") is a package that provides a variety of functions in support of linear and generalized linear models, including regression diagnostics (e.g., studentized residuals, hat-values, Cook's distances, dfbeta, dfbetas, added-variable plots,

Dear list members, I've attached some email correspondence with Jaro Lajovic (with his permission), detailing a problem with the Slovenian translation file for the Rcmdr package. In brief, while certain UTF-8 characters used in Slovenian used to appear properly in older versions of R, some characters do not display properly in the Rcmdr menus and output window under R 2.7.x. I've

Dear list members, The polr() function in the MASS package takes an optional weights argument for case weights. Is there any way to retrieve the case weights from the fitted "polr" object? Examining both the object and the code, I don't see how this can be done, but perhaps I've missed something. Any help would be appreciated. John ------------------------------ John Fox,

Dear John Fox, I am using Snowleopard with a Intel Core 2 Duo processor. I have installed the R console and followed your steps at: http://socserv.mcmaster.ca/jfox/Misc/Rcmdr/installation-notes.html but every time I try "library(Rcmdr)" it loads Tcl/tk infinitely. I have tried loading "library(tcltk)" by itself and it does the same thing. I have tried with X11 running to

Dear Roman, I've attached the latex file for the copy-edited version of your manuscript -- correcting English grammar and style, and making the paper conform more closely to the R News style. Can you check this over to confirm that I've not altered the meaning of your text and to approve the copy-edits? Thanks, John -------------------------------- John Fox Department of Sociology

Dear r-help members, Can lmer() in the lme4 package fit models that have a specified within-group correlation structure, as provided, for example, by the correlation argument to lme() in the nlme package? Thanks, John -------------------------------- John Fox Department of Sociology McMaster University Hamilton, Ontario Canada L8S 4M4 905-525-9140x23604 http://socserv.mcmaster.ca/jfox

Dear list members, I think that it's a pity that there's no Windows binary for the rgl package on CRAN. I understand that rgl doesn't compile with the standard Windows tools and that there's a Windows binary on Daniel Adler's web site. But the lack of an rgl binary on CRAN provides presents a potential obstacle to at least some Windows users, such as many students in

Dear r-devel list members, By placing a call to the new (with R 2.15.1) globalVariables() in the development version of the Rcmdr package, I've succeeded in reducing the "notes" produced by R CMD check from dozens to two. But I still get the following notes, even though '.commander.done' is in the call to globalVariables(): ------------- snip ------------ * checking R code

Dear list members, The release of R 2.9.0, reminds me of a long-standing nit that I have to pick. I prefer not to transfer and update all of the packages in my old library, because I see a new release of R as an opportunity to start with a clean slate. I'd rather install packages as I need them. My personal solution is to use the following function in place of library(): package <-

Dear r-devel list members, I find that adding --sdi to the Target field in the short-cut for R is ineffective in R 2.3.0 beta for Windows (version/platform details below). Setting MDI = no in the Rconsole file works, however. It's probably a bit late to request this for version 2.3.0, but it would be nice to be able to select either the SDI or MDI in the R Windows installer. Regards, John

I am trying to run a confirmatory factor analysis using the SEM package. My data are ordinal. I have read http://socserv.mcmaster.ca/jfox/Misc/sem/SEM-paper.pdf. When I apply the hetcor function, I receive the following error: Error in checkmvArgs(lower = lower, upper = upper, mean = mean, corr = corr, : at least one element of 'lower' is larger than 'upper' Example:

Dear r-devel members, My Rcmdr package "depends" on several other packages (tcltk, grDevices, utils, and car) and "suggests" a number of others (abind, aplpack, colorspace, effects, foreign, grid, Hmisc, lattice, leaps, lmtest, MASS, mgcv, multcomp, nlme, nnet, relimp, rgl, and RODBC). The reason for the distinction is that I don't want all of these packages to load when

Dear list members, In sem, std.coef() will give me standardized coefficients from a sem model. But is there a trick so that path.diagram can use these coefficients rather than unstandardized ones? Thanks Steve Powell From: John Fox <jfox_at_mcmaster.ca> Date: Wed 28 Feb 2007 - 14:37:22 GMT Dear Tim, See ?standardized.coefficients (after loading the sem package). Regards, John John

Dear list members (esp. Duncan), I've run into the following curious problem with Version 2.2.0 alpha for Windows: Options in the Profile file in R's etc directory don't appear to be set. (I habitually uncomment options(chmhelp=TRUE), to no effect in this case.) As far as I can see, the only thing non-standard about my installation is that I put R 2.2.0 alpha under c:\R rather than

Hi, In this newsletter (Vol 7, 1),the article on AMMI by Onofri and Ciriofolo presented a AMMI function. One of arguments for this function AMMI (Page 17) is biplot. There is a biplot fucntion from {stats} package. I guess they are not the same. Could the authors clarify that? Thanks, Yuandan [[alternative HTML version deleted]]

Dear Ashim, I?ll address your questions briefly but they?re really not appropriate for this list, which is for questions about using R, not general statistical questions. (1) The relevant distribution is within cells of the wool x tension cross-classification because it?s the deviations from the cell means that are supposed to be normally distributed with equal variance. In the warpbreaks data

Dear r-helpers, I'm looking for a way of sorting variable names in a "natural" order, when the names are composed of digits and other characters. I know that this is a vague idea, and that sorting character strings is a complex topic, but perhaps a couple of examples will clarify what I mean: > s <- c("x1b", "x1a", "x02b", "x02a",