similar to: R CMD build produces tar error under FreeBSD 5.4

Hi everyone, I was wondering if anyone has any code they could share for creating thumbnail plots in Sweave. For example, I'd like a plot like the following: y <- c(40, 46, 39, 44, 23, 36, 70, 39, 30, 73, 53, 74) x <- c(6, 4, 3, 6, 1, 5, 6, 2, 1, 8, 4, 6) opar <- par(mar=c(3,3,0,0)) plot(x, y, xlab="", ylab="") abline(h=mean(y), col="red") par(opar)

Dear community, can anyone provide a snippet of code to force the lm() to fit a model with terms in the formula in an arbitrary order? I am interested in something like: lm(y ~ A * B + C, data=data) where the interaction of A and B should be in the formula before C. My goal is to simplify my presentation of models using the anova() statement. I have found that this should be possible using

Dear R community, I'm trying to figure out if there is any way to split the printed output of some commands, for example summary.lme, so that I can intersperse comments in Sweave. I don't mind running the command numerous times and masking various portions of the output, or saving the output as an object and printing it, but I can't figure out how to do either. Does anyone have any

Hi Developers, The alpha, compiles successfully, but it is failing make check-all (on two seperate machines, both FreeBSD 6.1). Here is the version string: platform i386-unknown-freebsd6.1 arch i386 os freebsd6.1 system i386, freebsd6.1 status alpha major 2 minor 3.0 year 2006 month 03 day 27 svn rev

Dear R community, I'm fitting a complex mixed-effects model that requires numerous iterations and function evaluations. I note that nlminb accepts a list of control parameters, including eval.max. Is there a way to change the default eval.max value for nlminb when it is being called from lme? Thanks for any thoughts, Andrew -- Andrew Robinson Department of Mathematics and Statistics

Dear R developers, is there a preferred format or strategy for making a patch to contribute to a package that is maintained by R-core? Berwin Turlach and I have written a very minor extension to lmeControl to allow it to pass an argument to nlminb for the maximum number of evaluations of the objective function. I've edited the nlme/R/lme.R and nlme/man/lmeControl.Rd files. I can diff the

Dear Community, I'm interested in developing a package that could ease the command-line learning curve for new users. It would provide more detailed syntax checking and commentary as feedback. It would try to anticipate common new-user errors, and provide feedback to help correct them. As a trivial example, instead of > mean(c(1,2,NA)) [1] NA we might have > mean(c(1,2,NA)) [1]

The Rassist package has been loaded to CRAN. This package is designed to make R easier for new users, by providing extra checks and feedback. Presently the package functionality includes: * offers an alternative help facility, eg(.), with examples first, with additional examples included. eg() offers a start help menu, and eg(.) incorporates help.search(.) automatically. It also

The Rassist package has been loaded to CRAN. This package is designed to make R easier for new users, by providing extra checks and feedback. Presently the package functionality includes: * offers an alternative help facility, eg(.), with examples first, with additional examples included. eg() offers a start help menu, and eg(.) incorporates help.search(.) automatically. It also

Greetings, friends in the R community, this is an OT question about statistics. Given four time series of events, what possibilities do I have to test for synchronicity? e.g. times <- data.frame(year= c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10), event.1=c(1, 0, 0, 1, 2, 4, 1, 0, 0, 0), event.2=c(0, 0, 0, 1, 0, 2, 1, 0, 0, 0), event.3=c(1, 0, 0, 0, 1, 2, 4, 1, 0, 0),

Dear community, I'm trying to track down a quote, but can't recall the source or the exact structure - not very helpful, I know - something along the lines that: 80% of [applied] statistics is linear regression ... ? Does this ring a bell for anyone? Thanks, Andrew -- Andrew Robinson Department of Mathematics and Statistics Tel: +61-3-8344-9763 University of

Hi everyone, I'm witnessing an odd modelling phenomenon that I can't explain. If anyone has seen this before, or can explain what's going on would let me know, I'd be very grateful! Especially if I'm just being dim. I'm fitting a pair of continuous variates and their interaction to some residuals from another model. The sequential anova statement changes with the term

Hi dear R community, I have up to 12 measures of a protein for each of 6 patients, taken every two or three days. The pattern of the protein looks periodic, but the height of the peaks is highly variable. I'm testing for periodicity using a Monte Carlo simulation envelope approach applied to a cumulative periodogram. Now I want to predict the location of the peaks in time. Of course, the

Hi dear R community, I have up to 12 measures of a protein for each of 6 patients, taken every two or three days. The pattern of the protein looks periodic, but the height of the peaks is highly variable. It's something like this: patient <- data.frame( day = c(1, 3, 5, 8, 10, 12, 15, 17, 19, 22, 24, 26), protein = c(5, 3, 10, 7, 2, 8, 25, 12, 7, 20, 10, 5) ) plot(patient$day,

Hi R-help citizens, I'm having trouble making version 0.99-6 of Matrix on FreeBSD 6.0. The error message is: * Installing *source* package 'Matrix' ... ** libs gcc -I/usr/local/lib/R/include -I/usr/local/include -D__NO_MATH_INLINES -fPIC -g -O2 -c Csparse.c -o Csparse.o ... numerous lines deleted ... gcc -I/usr/local/lib/R/include -I/usr/local/include -D__NO_MATH_INLINES

Full_Name: Byron Ellis Version: R 1.9.0 (and 2.0.0) OS: Linux (Redhat Fedora Core) Submission from: (NULL) (140.247.241.197) It appears that `R CMD build` cannot handle spaces in the path when building packages for distribution. For instance: [ellis@net-78815 ~/Bayesian Networks]$ R CMD build bnsl * checking for file 'bnsl/DESCRIPTION' ... OK * preparing 'bnsl': * cleaning src *

Hi Community, I'm trying to reconcile Cook's Distances computed in glm. The following snippet of code shows that the Cook's Distances contours on the plot of Residuals v Leverage do not seem to be the same as the values produced by cooks.distance() or in the Cook's Distance against observation number plot. counts <- c(18,17,15,20,10,20,25,13,12) outcome <- gl(3,1,9)

Hi everyone, I'm involved in a discussion with a colleague. He suggested a sample design for a finite-sized process that (to all intents and purposes) involves tossing a coin and examining the unit if the coin shows Heads. I should emphasize that we're both approaching the problem from a design-based sampling theory point of view. So I have no argument about the appropriateness of the

Dear R community, I'm curious to know how people go about estimating standard errors for parameter estimates after model selection by ridge regression and the lasso. Do you have any practical or theoretical advice? Warmly, Andrew -- Andrew Robinson Department of Mathematics and Statistics Tel: +61-3-8344-9763 University of Melbourne, VIC 3010 Australia Fax:

Dear R-helpers, when in the browser, is it possible to change the environment, so as to be able to easily access (print, manipulate) objects in the parent, or elsehwere? I know that it is possible to evaluate expressions in different environments, using eval(), but I would prefer to avoid that if possible. Thanks, Andrew -- Andrew Robinson Department of Mathematics and Statistics