similar to: [OT] testing for synchronicity

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

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

Hi R-helpers, I am trying to build a package under FreeBSD 5.4-RELEASE #0 using R Version 2.1.1. I have constructed a package using package.skeleton(), when I try $ R CMD build foo * checking for file 'foo/DESCRIPTION' ... OK * preparing 'foo': * checking DESCRIPTION meta-information ... OK * cleaning src * removing junk files tar: Option -L is not permitted in mode -x Error:

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

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

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

Dear all, I'm curious about how to interpret the results of the following code. The first model is directly from the help page of lmer; the second is the same model but using the Gamma family with log link. The fixed effects make sense, because y = 251.40510 + 10.46729 * Days is about the same as log(y) = 5.53613298 + 0.03502057 * Days but the random effects seem quite

Dear Friends and Colleagues, I'm puzzling over how to interpret or use some bootstrap intervals. I think that I know what I should do, but I want to check with knowledgeable people first! I'm using a studentized non-parametric bootstrap to estimate 95% confidence intervals for three parameters. I estimate the variance of the bootstrap replicates using another bootstrap. The script

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