similar to: Coding categorical -> dummy

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

Many thanks for your very useful comments and suggestions. Renaud 2006/5/30, Prof Brian Ripley <ripley at stats.ox.ac.uk>: > On Tue, 30 May 2006, Prof Brian Ripley wrote: > > > This is not really a bug. See > > > > http://developer.r-project.org/model-fitting-functions.txt > > > > for how this is handled in other packages. All model-fitting in R used =

I met today a computer crash and our maintenance officer had to reinstall some components of the OS (MS Windows XP Pro) as well as the Internet browser (among other things). Now, I cannot install packages from a distant repository: > utils:::menuInstallPkgs() Error in .readRDS(pfile) : unknown input format > traceback() 5: .readRDS(pfile) 4: .packages(all.available = TRUE) 3:

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 all, I've just downloaded R.12.0 from CRAN (http://cran.r-project.org/bin/windows/base/release.htm). I am unable to find RGUI.exe as well as Rterm.exe. Conversely to what is explained in readme.R.2.12.0: There are two versions of the R executable in R-2.12.0\bin\i386 (32-bit) or R-2.12.0\bin\x64 (64-bit). I can't find directories i386 or x64 in R-2.12.0\bin -- Renaud Lancelot EDEN

Dear all, I have installed the latest version of R 2.12.0 available on CRAN (http://cran.r-project.org). When I try to load the recommended package lattice: > library(lattice) Error: package 'lattice' is not installed for 'arch=i386' I am running Rgui using C:\R\R-2.12.0\bin\i386\Rgui.exe --vanilla > sessionInfo() R version 2.12.0 (2010-10-15) Platform:

"Writing R Extensions" says {quotes} If you have a crash which gives a core dump you can use something like gdb /path/to/R/bin/exec/R core.12345 to examine the core dump. If core dumps are disabled... {unquotes} sadly it doesn't go on to say how to enable if core dumps are disabled. I understand that in bash I need to do $ ulimit -c unlimited but this doesn't seem to

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

> 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

Hello, All: The vignette with the sos package used "upquote.sty", required for R Journal when it was published in 2009. Current CRAN policy disallows "upquote.sty", and I've so far not found a way to pass "R CMD check" with sos without upquote.sty. I changed sos.Rnw per an email exchange with Prof. Ripley without solving the problem; see below. The

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

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, 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, When I use the cacheSweave package together with the Sweave option keep.source = TRUE, all the LaTeX code before the Sweave code chunk is included in the TeX file to be compiled. For example, with the following Sweave file \documentclass[12pt]{article} \usepackage[nogin]{Sweave} % Sweave options \SweaveOpts{keep.source=TRUE, strip.white=TRUE, eps=FALSE, pdf=TRUE} \begin{document}

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