similar to: Bug in plot.lm function (stats package): positioning of labels for extreme points.

I propose the following enhancements and changes to plot.lm(), the most important of which is the addition of a Residuals vs Leverage plot. (1) A residual versus leverage plot has been added, available by specifying which = 5, and not included as one of the default plots. Contours of Cook's distance are included, by default at values of 0.5 and 1.0. The labeled points, if any, are those

Hi, I'd like to request that the 'write.dta' command also save the "var.labels" attributes of the data frame. 'read.dta' can read them, but 'write.dta' doesn't write them back. Thanks.... Edward McNeil -- Epidemiology Unit Faculty of Medicine Prince of Songkla University Hat Yai THAILAND (\_ _/) (='.'=) ('')_('') -- This

Dear r-helpers, I estimated a generalized additive model (GAM) using Hastie's package GAM. Example: gam1 <- gam(vegetation ~ s(slope), family = binomial, data=aufnahmen_0708, trace=TRUE) pred <- predict(gam1, type = "response") vegetation is a categorial, slope a numerical variable. Now I want to assess the accurancy of the model using k-fold cross validation. I found the

Dear All, I am not able to run the random forest with my dataset.. X<- 280 records with satellite data(28 columns) - B1min, b1max, b1std etc.. y<- 280 records with 3 columns - TotBasal Area, Stem density and Volume yref <- y[1:230,] #Keeping 1st 230 records as reference records want to set 0 to y values for records 231 to 280.. yimp <- y[231:280,] #records for which we want

I'm using randomForest in yaImpute to create a yai-type object which associates "L" with landscape features. Then I use the sp() package to impute L to a landscape consisting of four ascii files) I keep getting the message "NA's generated due to illegal level(s)" when I do the imputation. It's probably because one of the landscape features ("as", for

Dear fellow R-users I am trying to do some imputation using K-NN with the yaImpute library. All seems to be going well until I try to use AsciiGridImpute. All my data are correctly formatted and I am able to run and view the results of yai. Below is my code: **************************************************************************************************** library("yaImpute") data

Dear All, I would like to plot text with a box around it. I used strwidth and strheight to compute the size of the box which is plotted with rect: z <- rnorm(10) # horizontal text works plot(rnorm(10)) x1 <- 5 y1 <- 0 label <- "Label" cha <- paste(" ", label, " ", sep = "") xh <- strwidth(cha, cex = par("cex")) yh <-

My mistake, but also a bug. I did not open a device and called strheight("X"). ==> R crashed ...\rgui --vanilla [...] > strwidth("X") Error in strwidth("X") : plot.new has not been called yet (OK.) > strheight("X") R crashed. (Happens also on Solaris.) Regards, Uwe Ligges --please do not edit the information below-- Version: platform =

What version of GCC? We user GCC 3.4.0. -G > -----Original Message----- > From: Marc Schwartz [mailto:MSchwartz@MedAnalytics.com] > Sent: Monday, December 20, 2004 4:30 PM > To: Warnes, Gregory R > Cc: Paul Murrell; Jain, Nitin; R-Devel; Frank E Harrell Jr > Subject: RE: [Rd] RE: [R] SAS or R software > > > On Mon, 2004-12-20 at 16:13 -0500, Warnes, Gregory R wrote:

I just completed YAI. The install was done via HTTP from an HTTP server (duh!). Is there some file I could copy into the directory structure (anaconda-ks.cfg) so that once I point the install to that server, it will do the second install just like the first?

If these things bother you, could you please supply fixes? For some reason you omitted the part of the help page dealing with one of these (which is as documented), and I also think you are confusing `plot' and `device' sizes. The latter will be the displayed device surface, I believe. Nothing on the help page suggested that the fine details of the R read-only parameters will be as you

Yes, par(family="mono") would work, except that I get R segfaults from this sequence: > > plot.new() > par(family="mono") > par(cex=8) > strheight("foo") Process R segmentation fault (core dumped) at Mon Dec 20 16:07:56 2004 on R 2.0.1 (2004-11-15), Red Hat Enterprise Linux AS release 3 In my code I call strheight and strwidth several times in

############################################################################### Before reporting 4 problems with windows(rescale=) I want to congrat on R1.2 and to thank r-developers for quickly adding the rescale workaround to the windows version. Happy New Year Jens Oehlschlaegel ###############################################################################

If I have drawn a string with ``ltext( x, y, labels="first string" )'' how can a draw a second string just one line (or strheight("X") below the first string regardless of the size and scales of the panel? Thanks Wolfram

Hi I am trying to perform a bootstrap estimate of classification accuracy of a logistic regression using the 'Daim' package in r using the code at the bottom of this post, this all works great and I get the .632+ misclassification accuracy, specificity, sensitivity, AUC etc etc but what I would like is to access the list of AUC for each of the bootstrap samples as I need calculate the 95%

Dear all, the help of par() claims that cin and cra are c(width, height) but it appears to be rather c(height, width) Best regards Jens Oehlschl?gel > plot.new() > strheight("W", unit="inches") [1] 0.1354167 > par("cin") [1] 0.1354167 0.1875000 > version _ platform i386-pc-mingw32 arch i386 os mingw32

Dear HelpeRs, I would like to do sg similar to: plot(c(0,100), c(0,100), xaxs="i", yaxs="i", type="n") txt <- paste("a =", a, "%\n b =", b, " km2", sep="") text(95, 95, txt, adj=c(1,1)) just with the "km2" formatted with 2 in a superscript. I thought txt <- substitute("a = "*a*" %\n b

I would suggest that you make a more thorough search of the R-Archives. (http://finzi.psych.upenn.edu/search.html) If you do you will find this discussion has been had several times and that the type of machine you are running will have an impact upon what you can do. My feeling is that you are going have to knuckle down with the documentation and understand how R works and then when you have

Dear All, I am Ricky Jacob, a project Student from India who is working on Forest Inventories. Input data: Plot(area = .1 ha) data having the following information: 1) Basal Area 2)Tree Density 3)Volume So I am applying this information to the corresponding pixels in the satellite imagery of the study area. I also have given the same values to a 3x3 window around that pixel. Inventory is taken

> > > > I've tested values of cex. Anything above 2.1249 seems to > segfault. > > > The X11 device works pretty hard to find a font to use and if > it really > can't find anything it should throw an error, not segfault. Are you > able to debug this and see where it is failing? I haven't had time to get gdb correctly pointed at a source tree