similar to: new lattice/grid

Hello all, I have a very simple set of data and I would like to analyze them with probit analysis. The data are: X Event Trial 100 8 8 75 8 8 50 6 8 25 4 8 10 2 8 0 0 8 I want to estimate the value of X that will give a 95% hit rate (Event/Trial) and the corresponding 95% CI. Anyone can offer some help? Thanks!! -

I just ran into this strange behavior. y <- rnorm(1000) x <- seq(0, length=length(y)) plot(x, y, type='l', lty=1) - instantaneous plot(x, y, type='l', lty=2) - 18s plotting lines + 15s plotting axes = 33s plot(x, y, type='l', lty=3) - 76s ,, + 75s ,, = 151s plot(x, y, type='l', lty=4) -

Is there any way to process the "vignette" documantation files? I have been looking at the bioconductor packages and most of them have *.Rnw files in doc directories. They look like regular TeX but have sections showing R code that seem to choke TeX. Thanks, Andy __________________________________ Andy Jaworski Engineering Systems Technology Center 3M Center, 518-1-01 St. Paul, MN

Hi, I am not sure if this is a bug, so I am posting it here. I also searched the FAQ file for the readline string and did not find anything about the behavior below. I just successfully rsync'ed, compiled and tested R-1.4.0 on my RH 7.2 Linux machine. I encountered one problem though. When running ./configure I noticed that I was getting errors from the readline library. Although

Hi, I am having the following problem and I am not sure if this is bug or I am missing something obvious. I have a bunch of R functions that automatically process large data sets. The processing produces several numerical summaries and several graphs which are subsequently put in reports. I usually do my processing in R BATCH mode on a Win2000 machine. All the graphical functions in my package

Hi, I have a problem of estimating a mixture of two normal distributions. I need to find the starting points automatically, since this is a part of a larger piece of image processing code. I found the mix2normal1 function in VGAM package that mentions a method of finding starting values for mu1 and mu2 but refers the reader to a book by Everitt and Hand. Unfortunately, I do not have an easy

Hi, This seems a very basic problem, but I cannot seem to find the solution to it. I am trying to use Rcmdr on a data frame I created in my current R session. More specifically, I did x <- data.frame(matrix(0, ncol=3, nrow=5)) library(Rcmdr) Now I would like to be able to edit, view and perhaps analyze x, but I cannot seem to figure out how to make it active. Any hint will be welcome.

Hi, I am trying to compile R-1.3.1 under RH-7.2 Linux. I used the old gcc compiler (kgcc) and the compilation went fine. I can start R and run some rudimentary tasks with no problem. However, when I did make check I got the following error in base-Ex.Rout > fa <- file.access(dir(".")) Error in file.access(dir(".")) : file.access is not implemented on

I am not sure if this is already known. I checked BUGS and found some references to similar behavior when R_HOME is set? Here is what I get: (1) Fresh build of 1.4.1 under Linux/RH-7.2. All tests pass. > R CMD check ctest * checking for working latex ... OK * using log directory `/home/andyj/stat/R-1.4.1/library/ctest.Rcheck' Installing *source* package `ctest'

I think that there is a bug in the as.POSIXct function on Windows. Here is what I get on Win2000, Pentium III machine in R 1.8.1. > dd1 <- ISOdatetime(2003, 10, 26, 0, 59, 59) > dd2 <- ISOdatetime(2003, 10, 26, 1, 0, 0) > dd2 - dd1 Time difference of 1.000278 hours Now, the 26th of October was the day that change to the standard time occurred, so I suspect that this has

Hi, I just installed the odesolve package and ran the check command on it. It failed trying to execute library(gnlm). Sure enough, there is no gnlm library on my system and I could not find it on the CRAN archive either. Am I missing anything obvious or is gnlm some private library that somehow found its way into the example section of odesolve? Thanks in advance, Andy

Hi, I have the following question about creating data frames. I want to create a data frame with 2 components: a vector and a matrix. Let me use a simple example: y <- rnorm(10) x <- matrix(rnorm(150), nrow=10) Now if I do dd <- data.frame(x=x, y=y) I get a data frame with 16 colums, but if, according to the documentation, I do dd <- data.frame(x=I(x), y=y) then str(dd)

I also downloaded the rgl package (after reading recent Uwe Ligges' post) and installed in on my Win2k Pro sp2 machine. It seems to work with no problems. I ran most of the examples, including the hist3d one, and all of them worked fine. A couple of thoughts. I have administrative priviliges on my machine. Could this possible make a difference. The rgl library uses OpenGL. Is it a

I wonder how one in R can fit a 3rd degree polynomial to some data? Say the data is: y <- c(15.51, 12.44, 31.5, 21.5, 17.89, 27.09, 15.02, 13.43, 18.18, 11.32) x <- seq(3.75, 6, 0.25) And resulting degrees of polynomial are: 5.8007 -91.6339 472.1726 -774.2584 THanks in advance! -- Jonas Malmros Stockholm University Stockholm, Sweden

I am about to get a new machine at work - an IBM Intellistation with the Xeon 2.8 GB processor. It will run Windows 2000. I would like to install the proper ATLAS dll for this machine, but I am not sure if Xeon is P4? Does anybody have any experience with Xeon? Thanks in advance, Andy __________________________________ Andy Jaworski 518-1-01 Process Laboratory 3M Corporate Research Laboratory

This is (hopefully) a small request. I routinely build (and test) both R versions from sources on my system (WinXP). I would like to be sure that I have the current version of Rtools. I always go to Robert Murdoch's site "Building R for Windows" and check the latest news. Everything I need is explained there, except I never know if I have the most recent version of Rtools.exe.

Hi, This happens in 2.8 and 2.9. I am posting this here because I am not sure if one should use bug.report for reporting bugs in unreleased versions. I was trying to run the example from the help page of plot.rpart. I got > library(rpart) > fit <- rpart(Price ~ Mileage + Type + Country, cu.summary) > plot(fit, compress=TRUE) Error in get(getOption("device")) : invalid

A strange thing happened. After installing the most recent RC version of 2.7 (2008-04-20 r45403) I compiled it as usual with no errors. After installing it I ran "make fullcheck". I got as far as testing recommended packages and I got the following: ... -------- Testing package foreign -------- Running examples in 'foreign-Ex.R' ... Comparing `foreign-Ex.Rout' to

Hello, I need help with the details of loess prediction algorithm. I would like to get it implemented as a part of a measurement system programmed in LabView. My job is provide a detailed description of the algorithm. This is a simple one-dimensional problem - smoothing an (x, y) data set. I found quite a detailed description of the fitting procedure in the "white book". It is also

Hello, I faced something I can't understand. When I use boxplot.stats(1:10) and quantiles(1:10) the results are different for 25% and 75%: > boxplot.stats(1:10) $stats [1] 1.0 3.0 5.5 8.0 10.0 > quantile(1:10) 0% 25% 50% 75% 100% 1.00 3.25 5.50 7.75 10.00 Actually, I expected the value 3 for 25% and 8 for 75% as results of quantile(1:10). Can you please explain me