search for: chabotd

Hi, Without trying to print 1000000 points (see <http:// finzi.psych.upenn.edu/R/Rhelp02a/archive/42105.html>), I often print maps for which I do not want to loose too much of coastline detail, and/or plots with 1000-5000 points (yes, some are on top of each other, but using transparency (i.e. rgb colors with alpha information) this actually comes through as useful information.

Hi, Can we attach a more descriptive "label" (I may use the wrong terminology, which would explain why I found nothing on the FAQ) to variable names, and later have an easy way to switch to these labels in plots? I fear this is not possible and one must enter this by hand as ylab and xlab when making plots. Thanks in advance, Denis Chabot

Hi, I am fairly new to GAM and started using package mgcv. I like the fact that optimal smoothing is automatically used (i.e. df are not determined a priori but calculated by the gam procedure). But the mgcv manual warns that p-level for the smooth can be underestimated when df are estimated by the model. Most of the time my p-levels are so small that even doubling them would not result

Hi, I don't understand the following. When I create a small artificial set of date information in class POSIXct, I can calculate the mean and the median: a = as.POSIXct(Sys.time()) a = a + 60*0:10; a [1] "2009-03-10 11:30:16 EDT" "2009-03-10 11:31:16 EDT" "2009-03-10 11:32:16 EDT" [4] "2009-03-10 11:33:16 EDT" "2009-03-10 11:34:16

Hi, I'm wondering why 48 %/% 2 gives 24 but 4.8 %/% 0.2 gives 23... I'm not trying to round up here, but to find out how many times something fits into something else, and the answer should have been the same for both examples, no? On a different topic, I like the behavior of NAs better in R than in SAS (at least they are not considered the smallest value for a variable), but at the

Hi, I have used package quantreg to estimate a non-linear fit to the lowest part of my data points. It works great, by the way. But I'd like to extract the predicted values. The help for predict.qss1 indicates this: predict.qss1(object, newdata, ...) and states that newdata is a data frame describing the observations at which prediction is to be made. I used the same technique I used

Hi, Sometimes data series (not necessarily time series) suffer breaks where data were expected, but not collected. Often the regular "lines" command to add such data to a plot is what I want, but other times I'd like the line to break where the data series is interrupted, instead of the line jumping to the next point in the series uninterrupted. Usually my data file

Dear list, I've tried to search the archives but found nothing, although I may use the wrong wording in my searches. I've also double-checked the upData function in Hmisc, but it does something else. I'm wondering if one can update a dataframe by "forcing into" it a shorter dataframe containing the corrections, like the "update" provided in SAS data steps.

Hi, As part of the legend to a plot, I need to have the "n" in italics because it is a requirement of the journal I aim to publish in: "This study, n = 3293" Presently I have: legend(20, 105, "This study, n = 3293", pch=1, col=rgb(0,0,0,0.5), pt.cex=0.3, cex=0.8, bty="n") I suppose I could leave a blank in place of the "n",

Hi, I've looked at the parameters available for the legend function and cannot find a way to change the distance between the top of the box surrounding a legend and the legend's title. I have a math expression that raises the height of my title. If you don't mind the non-sensical title I give to the legend for this plot (Figure 3.20 in R Graphics): with(iris,

Dear list members, I have to combine dataframes together. However they contain different numbers of variables. It is possible that all the variables in the dataframe with fewer variables are contained in the dataframe with more variables, though it is not always the case. There are key variables identifying observations. These could be used in a merge statement, although this won't

Hi, I have data on stomach contents. Possible prey species are in the hundreds, so a list of prey codes has been in used in many labs doing this kind of work. When comes time to do analyses on these data one often wants to regroup prey in broader categories, especially for rare prey. In SAS you can nest a large number of "if-else", or do this more cleanly with "select"

Hello R users, I have to import a file with one column containing dates written in French short format, such as: 7-d?c-07 11-d?c-07 14-d?c-07 18-d?c-07 21-d?c-07 24-d?c-07 26-d?c-07 28-d?c-07 31-d?c-07 2-janv-08 4-janv-08 7-janv-08 9-janv-08 11-janv-08 14-janv-08 16-janv-08 18-janv-08 There are other columns for other (numeric) variables in the

Hi, I'm trying to move to R the last few data handling routines I was performing in SAS. I'm working on stomach content data. In the simplified example I provide below, there are variables describing the origin of each prey item (nbpc is a ship number, each ship may have been used on different trips, each trip has stations, and individual fish (tagno) can be caught at each

...at stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide! Message: 71 Date: Wed, 8 Mar 2006 10:49:28 +0000 (GMT) From: Simon Wood <sw283 at maths.bath.ac.uk> Subject: Re: [R] predicted values in mgcv gam To: Denis Chabot <chabotd at globetrotter.net> Cc: R list <r-help at stat.math.ethz.ch> Message-ID: Content-Type: TEXT/PLAIN; charset=US-ASCII Hi Denis, Your first plot is of f(x) against x where \sum_i f(x_i) = 0 (x_i are observed x's). Your second plot is of \exp(\alpha + f(x)) against x where \sum_i f(x_i...

Hi, I'd like to compute the weighted mean of some variables, all using the same weight variable, for each combination of 3 factor variables. I found how I could use "summarize" (from Hmisc) to do normal means for combinations of 3 factors, but I cannot find a way of doing weighted means. Is it possible in R? Thanks in advance, Denis Chabot

Hi, I do have it installed on 2 Macs as well (OS X 10.2.8 and 10.3.7) and what I need does work, however if you do the command check routine some problems will likely be revealed. At least there were problems for me. Denis Le 08 févr. 2005, à 12:23, r-help-request@stat.math.ethz.ch a écrit : > De: Don MacQueen <macq@llnl.gov> > Date: 07 février 2005 16:05:14 GMT+01:00 > À:

> De: Depiereux Constant <constant.depiereux at aqte.be> > Date: 12 mars 2005 09:44:08 GMT-05:00 > ?: r-help at stat.math.ethz.ch > Objet: [R] New user of R on Mac OS X - Please help > > > Brand new Mac OS X user, I am transfering my R stuffs from my windows > machine. > > When porting some of my functions, I got messages such as : > > 2005-03-12

Hi, I would like to create a new variable that would be the rank of a continuous variable within each level of a grouping variable. Say the first level of the grouping variable has 5 observations, I'd like them ranked from one to five and a new variable would hold the rank value (one to five). So forth for each level of the grouping variable. I'm quite new with R and cannot

Hi, My problem does not happen all the time, nor with all files I save to csv format. I can send my test file (format rda or csv) to whoever would like to replicate (and hopefully explain) my problem. In short, I have a dataset with mostly numerical variables. One of my variable is called pfi2 and is definitely numerical, as shown by this: > summary(test$pfi2) Min. 1st Qu. Median