similar to: Enhanced version of plot.lm()

Dear all, I noticed the following behaviour of plot.lm: > fm1 <- lm(time~dist, data=hills, weights=c(0,0,rep(1,33))) > par(mfrow=c(2,2)) > plot(fm1) Warning messages: 1: longer object length is not a multiple of shorter object length in: res/(sd * (1 - hat)) 2: longer object length is not a multiple of shorter object length in: (res/(sd * (1 - hat)))^2 * hat which seems to be

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

As some of you R-devel readers may know, the plot() method for "lm" objects is based in large parts on contributions by John Maindonald, subsequently "massaged" by me and other R-core members. In the statistics litterature on applied regression, people have had diverse oppinions on what (and how many!) plots should be used for goodness-of-fit / residual diagnostics, and to my

# Your mailer is set to "none" (default on Windows), # hence we cannot send the bug report directly from R. # Please copy the bug report (after finishing it) to # your favorite email program and send it to # # r-bugs@biostat.ku.dk # ###################################################### The error occurs following the third Hit <Return> below: > y_rnorm(40) >

The following code demonstrates an annoyance with plot.lm(): library(DAAGxtras) x11(width=3.75, height=4) nihills.lm <- lm(log(time) ~ log(dist) + log(climb), data = nihills) plot(nihills.lm, which=5) OR try the following xy <- data.frame(x=c(3,1:5), y=c(-2, 1:5)) plot(lm(y ~ x, data=xy), which=5) The "Cook's distance" text overplots the label for the point with the

Confirmation from John Maindonald. -- Brian D. Ripley, ripley@stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272860 (secr) Oxford OX1 3TG, UK Fax: +44 1865 272595 ---------- Forwarded message ---------- Date: Fri, 3

Hello, While using the "plot" function against a model with a long formula, I get a garbled label in the bottom of the plot. It looks like there are two lines printed one over another, here's my example: http://pico.magnum2.pl/maciej/NodalInvolvement/MisclassificationRate?action=AttachFile&do=get&target=TwoFactorResiduals.png This was from using

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

Hi, I just try to draw multiple plots in one page using R, I used par command. For example I have 7 plots, but instead of arranging them in the default way plot1 plot2 plot3 plot4 plot5 plot6 plot7 I want them in this order plot1 plot2 plot3 plot4 plot5 plot6 plot7 Could somebody tell me how to do this, please? Thanks so many. Suyan

Hi, Figure 13 of Emmanuel Paradis's "R for Beginners" was produced by termplot working on an aov object. The lower right-hand plot is labelled "Cook's distance plot", and I'd really like to produce a similar type of figure, but in a totally different context. (I'm not even sure what this kind of figure is called, perhaps an "impulse plot", where

(1) read.table(), with sep="\t", identifies 13 our of 1400 records, in a file with 1400 records of 3 fields each, as having only 2 fields. This happens under version 2.3.1 for Windows as well as with R 2.3.1 for Mac OS X, and with R-devel under Mac OS X. [R version 2.4.0 Under development (unstable) (2006-07-03 r38478)] (2) Using read.table() with sep="\t", the first 1569

Does anyone know what algorithm R uses in glm, family=binomial (i.e. a logit model)? I assume that it's in the source somewhere, but I wasn't able to find it. I'd like to know what file it's in (in a unix distribution of R). Thanks for your help. --------------------------- Barnet Wagman wagman at enteract.com 1361 N. Hoyne, 2nd floor Chicago, IL 60622 773-645-8369

If anyone has a list of application areas where there is extensive use of R, I'd like to hear of it. My current short list is: Bioinformatics Epidemiology Geophysics Agriculture and crop science John Maindonald Mathematical Sciences Institute, Australian National University. john.maindonald at anu.edu.au

Hi there, I need to save a series of lattice plots as a PDF, this is my code so far: windows(height=8,width=6) plot.new() library('grid') lattice.options(layout.heights=list(top.padding=list(x=0.15, units="inches"))) print(plot1, split=c(1,1,2,3), more=TRUE) print(plot2, split=c(1,2,2,3), more=TRUE) print(plot3, split=c(1,3,2,3), more=TRUE) print(plot4, split=c(2,1,2,3),

Can anyone suggest an equivalent, for lattice graphs, of the base graphics argument bty="l"? NB that I am leaving off the box around the strip, with a strip function: stripfun <- function(which.given,which.panel, factor.levels=as.expression(levlist), ...){ panel.text(x=0, y=0.5, lab = as.expression(levlist[which.panel[which.given]]),

Hi Deepayan and everyone, I need to add a common legend to a group of latice graphs, I have tried different ways using viewport and grid.draw without success. Here is what I have: plot.new() library(grid) library('IDPmisc') print(plot1, split=c(1,1,2,4), more=TRUE) print(plot4, split=c(2,1,2,4), more=TRUE) print(plot2, split=c(1,2,2,4), more=TRUE) print(plot5, split=c(2,2,2,4),

Does anyone know any more than is in the following press release about REvolution Computing and their commercialization of R? http://www.intel.com/capital/news/releases/080122.htm "Intel Capital, the global investment arm of Intel Corporation, today announced that it has invested in the Series A financing of REvolution Computing, creator of parallel computing software for computational

I am looking for information on experimentation with the use of R in the teaching of statistics and science in schools. Any leads would be very welcome. I am certain that there is such experimentation. I've made this inquiry on r-sig-teaching, with no response. John. John Maindonald email: john.maindonald at anu.edu.au phone : +61 2 (6125)3473 fax : +61 2(6125)5549 Centre

help(letters, htmlhelp=TRUE) fails. Under the Mac OSX gui, the message is 'Help for the topic "a" was not found.' Under the version documented below, and under Windows, the message is "No documentation for 'a' in specified packages and libraries:" repeated for all the elements of letters, then followed by "you could try

In Section 1.4.2 of "Writing R Extensions" %\VignetteEngine{knitr::knitr} should be %\VignetteEngine{knitr::knit} > sessionInfo() R version 3.0.2 (2013-09-25) Platform: x86_64-apple-darwin10.8.0 (64-bit) Is this sort of thing best reported here, or is a huge report in order? John Maindonald email: john.maindonald at anu.edu.au phone : +61 2 (6125)3473 fax