similar to: surfaces and digital terrain model

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

> 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 problems to install some of the packages dowloaded from CRAN (whatever the

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

Hi, lme() users, Can some one tell me how to do this. I model Orthodont with the same G for random variables, but different R{i}'s for boys and girls, so that I can get sigma1_square_hat for boys and sigma2_square_hat for girls. The model is Y{i}=X{i}beta + Z{i}b + e{i} b ~ iid N(0,G) and e{i} ~ iid N(0,R{i}) i=1,2 orth.lme <- lme(distance ~ Sex * age, data=Orthodont, random=~age|Subject,

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

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

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

thanks a lot Renaud. but i was interested in Finney's fiducial confidence intervals of LD50 so to obtain comparable results with SPSS. But your reply leads me to the next question: does anybody know what is the best method (asymptotic, bootstrap etc.) for calculating confidence intervals of LD50? i could "get rid" of Finney's fiducial confidence intervals but

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

Dear all, platform i386-pc-mingw32 arch x86 os Win32 system x86, Win32 status major 1 minor 2.3 year 2001 month 04 day 26 language R I wrote the following to import tables in different database formats (Access, Visual FoxPro, etc.): getTable <-

Dear all, I am trying to implement goodness-of-fit tests for logistic regression models described in : Hosmer, D.W., Lemeshow, S., 2000. Applied logistic regression. New-York, John Wiley & Sons, Inc., 373 p. (pp. 152-154) Namely I would like to reproduce Osius & Rojek's test (Osius, G., Rojek, D., 1992. Normal goodness-of-fit tests for multinomial models with large degrees of

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

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

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

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