search for: darrenleeweb

Hi, if we have a file called Rscript.R that contains the following, for example: x <- 1:100 outfile = "Rscript.Rout" sink(outfile) print(x) and then we run >> source("Rscript.R") we get an output file called Rscript.Rout - great! Is there an internal variable, something like .Platform, that holds the script name when it is being executed? I would like to use

Dear Prof. Ripley et al., I have a data frame with some variable names that contain the underscore character ('_'). It seems that RODBC is not converting them to a suitable variable name for MySQL. > sqlSave(channel, zz, addPK = TRUE, verbose = TRUE) Query: CREATE TABLE zz (rownames varchar(255) PRIMARY KEY, NAME varchar(255), ID varchar(255), GROUP varchar(255), OAC_RT double,

Hi, I can't find a simple command to switch the plot colors from paper to presentation formats. Has anyone defined styles that can be easily applied in one command? It would be nice to have a command for papers, using ps or eps outputs, that has a white background and colors that will print nicely on that. On the other hand, it would be nice to have exactly the same plot for a presentation

With respect to calculating the epsilon index of sphericity for ANOVA, discussed on pp. 45-47 of: http://www.psych.upenn.edu/~baron/rpsych.pdf It notes that epsilon is not required for a repeated measures design with only k=2 levels, as the minimum value of epsilon (e) is given by: e = 1/(k-1) so for k=2, we have e = 1 (ie, no correction of the F test df; see p. 46). These notes apply to a

Hi, are there any good charting and analysis tools for use with currencies, stocks, etc. in R? I have some tools to download currency data from the NYFRB using python and XML. Can we get and parse an XML download using R? Can we have interaction in R plots? Does anyone use R for back-testing trading strategies? Are there any forums for discussion of using R for this specific purpose (apart

Hi, I'm stuck on creating a plot with x tick labels for every Nth tick mark - how is that done? I don't see a simple solution to this in help(plot) or help(par) and what I've tried is not working, eg, the following does not work, although it seems intuitive to me that it should work: x <- seq(-100,1000,25) y <- x * x % find all the x values that are multiples of 100 tmp <-

Hi, where is the latest version of the RSPython library? Is it compatible with the current stable release of R? Thanks, Darren

Is there a function in R for doing Greenhouse-Geiser correction in ANOVA models? Is it already available in the aov function? How do we use it? Best, Darren [[alternative HTML version deleted]]

Dear Prof. Baron, thankyou for your pdf paper on using R, at http://www.psych.upenn.edu/~baron/rpsych.pdf I am currently interested in the sphericity notes on pp. 45-47. The equation on p. 46 suggests the following R calculation: D <- k^2 * ( mean(diag(S)) - mean(S) )^2 while the notes in the pdf indicate the following: D <- k^2 * ( mean(S) - mean(diag(S)) )^2 Can you please consider

For example, this works: x = seq(-100, 1000, 25) y = x * x plot(x,y, xaxt="n") axis(1,x,FALSE,tcl=-0.3) axis(1,x[x %% 100 ==0]) It creates two axis objects and the values of the x-axis are the labels. The following scenario is more difficult, because it uses 'image' to plot a grid of values: a = matrix(rnorm(100),ncol=10) image(1:ncol(a), 1:nrow(a), a, xaxt="n",

Dear Erich and John, thankyou for providing RColorBrewer! Are you working on options for CMYK or RGB return values? Thanks again, Darren [[alternative HTML version deleted]]

Why do I see this error? > library(stats) > require(stats) [1] TRUE > > tHSD <- TukeyHSD(aov) Error in TukeyHSD(aov) : no applicable method for "TukeyHSD" In case it helps: > aov Call: aov(formula = roi ~ (Cue * Hemisphere) + Error(Subject/(Cue * Hemisphere)), data = roiDataframe) Grand Mean: 8.195069 Stratum 1: Subject Terms: Residuals Sum

Hi, I find a lot of the R-help email traffic overloads my inbox. My IT managers are not really happy for me to be subscribed to several high-traffic email lists. I don't want to lose my contact with the R-help emails, so I'm having to consider various ways of handling the traffic. Anyhow, I'm wondering how many people on the R-help email list would prefer that most of the traffic

Hi, we need some help to define planned comparisons. I've based my understanding of the problem on reading Tabachnick and Fidell (2006), ie: http://www.ablongman.com/catalog/academic/product/0,1144,0205459382,00.html I don't understand how to specify planned comparisons in R. I've not found explanations for this in MASS or elsewhere. There is only discussion of the contrast

How do we specify a contrast interaction matrix for an ANOVA model? We have a two-factor, repeated measures design, with Cue Direction (2) x Brain Hemisphere(2) Each of these has 2 levels, 'left' and 'right', so it's a simple 2x2 design matrix. We have 8 subjects in each cell (a balanced design) and we want to specify the interaction contrast so that: CueLeft>CueRght

Hi, I have some general questions about statistical analysis for a research dataset and a request for advice on using R and associated packages for a valid analysis of this data. I can only pose the problem as how to run multiple ANOVA tests on time series data, with reasonable controls of the family-wise error rate. If we run analysis at many small sections of a long time-series, the Type-I