similar to: problem with package binaries

Dear all, Has anyone written functions to compute (and ideally project / plot onto a plane) confidence ellipsoids for q linear model parameters, q > 2 ? Greetings and best wishes, Renaud -- Dr Renaud Lancelot, v?t?rinaire CIRAD, D?partement Elevage et M?decine V?t?rinaire (CIRAD-Emvt) Programme Productions Animales http://www.cirad.fr/presentation/programmes/prod-ani.shtml ISRA-LNERV

Dear all, I have computed a digital terrain model from a set of points (x, y, z) using the function interp() in package akima. I want to predict flooded surfaces given target values of z. I can display the flooded surfaces with contour() or image(), but I don't know how to get the polygons delimiting the surfaces. Did anyone write a function for this purpose ? Best regards, Renaud --

Dear all, platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status Patched major 1 minor 5.0 year 2002 month 05 day 15 language R Is there a way to include a graph in a help file (Rd file processed with make pkg-mypkg) ? Thanks in advance, Renaud -- Dr

Dear R-Helpers, Is there a way to draw shaded (filled) polygons in panel graphs produced by xyplot ? I made different trials (see below) with polygon() and grid.polygon() but polygons are not drawn (probably for the same reason that prevents the use of lines, points, etc. in xyplot panels). Thanks in advance, Renaud > Water id label longitude latitude 515 21 Zone inondable

Dear all, 1. How can I compute the intersecting area between 2 polygons ? 2. I have polygons with coordinates in decimal degrees (i.e. 13 deg 30 min = 13.5 decimal degrees). I want to compute their area and get the results in square meters or square kiometers. Can anyone give me a conversion coefficient or a pointer where I can find this information (sorry for this off topic question) ? Thanks

Dear all, > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 5.0 year 2002 month 04 day 29 language R OS: Win98 I have a strange problem with a particular data.frame built with

I've put up a build of the current R-patched (and updated recommended packages) for Windows at http://www.stats.ox.ac.uk/pub/bdr/RWin/Rpre/SetupRpatched.exe This should be able to rename files on Win 9x (unlike the mid-May build) and can build packages accessing Rblas.dll. As ever, this is less tested than an official release, in particular only on Windows XP, but ought to have fewer

Please reply to: lancelot at sentoo.sn because my address: renaud.lancelot at cirad.fr is broken (should be fixed soon...). platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status Patched major 1 minor 5.0 year 2002 month 05 day 15 language R

To Windows R users: The binary build of R that I uploaded on Tuesday contained errors in the compiled HTML help files and some of the PDF versions of the manuals. In particular, the index was missing from the CHM files, and the table of contents was missing from at least one of the manuals. I've since found the error in my setup, and have rebuilt the distribution. I uploaded it to CRAN

To Windows R users: The binary build of R that I uploaded on Tuesday contained errors in the compiled HTML help files and some of the PDF versions of the manuals. In particular, the index was missing from the CHM files, and the table of contents was missing from at least one of the manuals. I've since found the error in my setup, and have rebuilt the distribution. I uploaded it to CRAN

Dear all, I'd like to draw a color key beside a graph drawn with xyplot (lattice library). I am aware of the draw.colorkey function (grid library) but don't know how to handle it. Any hint would be appreciated. Renaud -- Dr Renaud Lancelot, v?t?rinaire CIRAD, D?partement Elevage et M?decine V?t?rinaire (CIRAD-Emvt) Programme Productions Animales

Sorry if this is obvious, but my S skills aren't great and I haven't been able to find it documented anywhere. I want to write a new function for use with lme objects; the function will simply calculate an ICC (aka "rho") for each level of a mixed-effects model. What I need for this is pretty simple: (c(var1..varn, residual)) / sum(c(var1..varn, residual)) where var1..varn

I have a number of plots that I'm trying to do using the lattice package, unfortunately, they involve Posix values. A small sample of the data is as follows: Time TC.1 <Snipped> 1 2003-04-24 13:29:15 26.61 2 2003-04-24 13:29:30 26.48 3 2003-04-24 13:30:00 26.53 4 2003-04-24 13:30:30 27.85 <Snipped> where str(data) `data.frame': 415 obs. of

Please, forgive my ignorance on statistics, but I have a rather simple question concerning AIC. Small values of AIC mean good fit of the model. How about negative values of AIC? Is a model with AIC=0.5 considered to have better fit than a model with AIC=-500? I couldn't find anything explaining this in the documentation or any elementary statistics textbook that I have at my disposal. Thanks

Dear all, It looks like the parameters var-cov matrix returned by gee() is not adjusted for the scale parameter: > fm1 <- gee(nbtrp ~ strate * saison + offset(log(surf)), family = poisson, data = Eff2001, + id = loc, tol = 1e-10, corstr = "exchangeable") [1] "Beginning Cgee S-function, @(#) geeformula.q 4.13 98/01/27" [1] "running glm to get initial

I've been using par() to check the graphics parameters associated with both plot(<fitted linear model>) and plot(<grouped data object>). AFAIK the only differences are in the $cxy, $usr, $xaxp, and $yaxp parameters but the background color for the grouped data plot is grey while the linear model plot has a white background. When I've tried par(bg = "white") prior to

Hello R-users, I'm using the boot's library. I have two questions: 1. I execute for example data<-rnorm(100) boot(data,mean,500) and I obtain 500 equal bootstrap samples of data. Why? 2. I need to calculate the bootstrap confidence intervals for a statistic, but when I use the function boot.ci it shows the error message: Error in if

hi, during a for(i in 1:xx]) loop I always newly calculate a vector (e.g. tmp<-c(a,b,c,d,e) ) now i need that vector to be attached at the bootom of a matrix (or data.frame). e.g. m<-matrix() for(i in 1:5]){ #some calculations for a,b,c,d,e a<- ... b<- ... c<- .. .... tmp<-c(a,b,c,d,e) ??? now I need to attach this tmp to the matrix m ??? } i couldn't get rbind

Dear all, I am trying to avoid the warnings produced by: > x <- -2:2 > log(x) [1] NaN NaN -Inf 0.0000000 0.6931472 Warning message: production de NaN in: log(x) I thought that using ifelse would be a solution, but it is not the case: > ifelse(test = x < 0, yes = NaN, no = log(x)) [1] NaN NaN -Inf 0.0000000 0.6931472 Warning message: production

In Baron and Li's "Notes on the use of R for psychology experiments and questionnaires" http://cran.r-project.org/doc/contrib/rpsych.htm they describe a balanced data set for a drug experiment: "... a test of drug treatment effect by one between-subject factor: group (two groups of 8 subjects each) and two within-subject factors: drug (2 levels) and dose (3 levels). "