similar to: Lattice: merged strips?

Hi all, This is my first post to the developers list. As I understand it, aggregate() currently repeats a function across cells in a dataframe but is only able to handle functions with single value returns. Aggregate() also lacks the ability to retain the names given to the returned value. I've created an agg() function (pasted below) that is apparently backwards compatible (i.e.

Hi all, Consider the following function: #### my.func = function(x){ y=ifelse(x>-.5,0,ifelse(x< -.8,abs(x)/2,abs(x))) print(c(x,y)) #print what was tested and what the result is return(y) } curve(my.func,from=-1,1) #### When I attempt to find the maximum of this function, which should be -.8, I find that optimize gets stuck in the plateau area and doesn't bother testing the

Hi all, A suggestion derived from discussions amongst a number of R users in my research group: set the default column names produced by aggregate () equal to the names of the objects in the list passed to the 'by' object. ex. it is annoying to type with( my.data ,aggregate( my.dv ,list( one.iv = one.iv ,another.iv = another.iv ,yet.another.iv = yet.another.iv )

Hi all, I'm looking for an analytic method to obtain the PDF & CDF of the sampling distribution of a given correlation (rho) at a given sample size (N). I've attached code describing a monte carlo method of achieving this, and while it is relatively fast, an analytic solution would obviously be optimal. get.cors <- function(i, x, y, N){ end=i*N

I'm looking to analyze a large data set: a within-Ss 2*2*1500 design with 20 Ss. However, aov() gives me an error, reproducible as follows: id = factor(1:20) a = factor(1:2) b = factor(1:2) d = factor(1:1500) temp = expand.grid(id=id, a=a, b=b, d=d) temp$y = rnorm(length(temp[, 1])) #generate some random DV data this_aov = aov( y~a*b*d+Error(id/(a*b*d)) , data=temp ) While yields the

Hi all, Just thought I'd share something I discovered last night. I was interested in creating animations consisting of a series of plots and after finding very little in the usual sources regarding animation in R directly, and disliking the imagemagick method described here (http://tolstoy.newcastle.edu.au/R/help/05/10/13297.html), I discovered that if one exports the plots to a

Hi all, Hopefully this will be quick, I'm looking for pointers to packages/ functions that would allow me to calculate the power of a t.test when the DV has measurement error. That is, I understand that, ceteris paribus, experiments using measure with more error (lower reliability) will have lower power. Mike -- Mike Lawrence Graduate Student, Department of Psychology, Dalhousie

Hi all, Worked & looked around for a while on this to no avail. I'm trying to create a plotmath expression that achieves: ?i >> 0 and while: expression(Delta*i>0) comes close, I'd prefer to have the >> (denoting "very much greater than"). Maybe >> is a non-standard expression and therefore not supported? Mike -- Mike Lawrence Graduate

Hi all, I understand that it is simple to create data with a specific correlation (say, .5) using mvrnorm from the MASS library: > library(MASS) > set.seed(1) > > a=mvrnorm( + n=10 + ,mu=rep(0,2) + ,Sigma=matrix(c(1,.5,.5,1),2,2) + ,empirical=T + ) > a [,1] [,2] [1,] -1.0008380 -1.233467875 [2,] -0.1588633 -0.003410001 [3,] 1.2054727 -0.620558768

Hi all, I have collected response time data from 178 participants ('sub') for each combination of 4 within-Ss factors ('con','int','tone','cue'). Additionally, I have recorded the gender of each participant, so this forms a between-Ss factor ('gender'). Normally this would be analyzed using aov:

Has anyone successfully compiled pnmath (http://www.stat.uiowa.edu/~luke/R/experimental ) for an intel processor running mac OS 10.5? When I attempt to do so via the R package installer (choosing "Local Source Package" and pointing to the pnmath_0.0-2.tar.gz file), I get the following errors: * Installing *source* package 'pnmath' ... ** libs ** arch - i386 gcc -arch i386

Hi all, I am extremely newbie to R. Can anybody jump-start me with any clues as to how do I get a cumulative histogram from two independent variables, cumhist(X,Y) ? -jose [[alternative HTML version deleted]]

Hello, i'm trying to find a more modern system to reproduce the functionality that was available through the Histoscope program (from Fermilab). Namely, the capability of connecting to a running process and having plots update in realtime in response to new data. Is this possible with R? Thank you, Charles Cosse [[alternative HTML version deleted]]

#Quick one hopefully. Shouldn't dnorm be returning the pdf? Last time I checked, #a probability shouldn't be greater than 1 as produced by: curve(dnorm(x,0,.1),from=-1,to=t) #Shouldn't I be getting an axis more like that produced by: f=function(x,m,s){ y=rep(NA,length(x)) for(i in 1:length(x)){ y[i]=integrate( dnorm , upper=x[i]+sqrt(.Machine$double.eps) ,

I am using R 2.4.0 and lattice to produce some xyplots conditioned on a factor and a shingle. The shingle merely chops up the data along the x-axis, so it is easy to identify which part of the shingle a panel is in by looking at the x-axis markings. I only want to have a strip at the top for the factor. Is this possible? I looked into calculateGridLayout() and it seems to me that there

Dear all, back in 2002 Martin Henry H. Stevens wrote (https://stat.ethz.ch/pipermail/r-help/2002-May/019851.html) > How do I control the text in strips? Specifically, I want to put in the > ranges generated in shingle(x) where x is continuous. with an answer from Deepyan Sarkar (see strip.new towards the end of this message). I assume that the answer worked back then, but I've tried

Hello everybody My collegue and I noticed a strange behaviour of R on different platforms. It's a simple computation, but results are rather different. On Windows XP: > floor(log(8,2)) [1] 3 which is what one should expect. Here's instead the result with Mac OS X (same version, 2.5.0 (2007-04-23)) > floor(log(8,2)) [1] 2 Is it a "bug" in R or in the operating

Hello (R-)Experts I hope someone can help with this problem concerning axis annotation of a lattice shingle plot. I want a plot with three shingles to display some laboratory value over time. In the first panel over the first few days, then in the next panel some months, and in the last panel some years. In the following minimal example the axis annotation will be in days, but I'd like to

I'm interested in a version of tapply() that operates with shingles instead of factors. For instance: x <- c(1,1,2,2,3,3) y <- c(1,1,1,0,0,0) s <- shingle(x,intervals=cbind(c(0.5,1.5),c(2.5,3.5))) # the following function should exist! tapply.shingle(x,s,mean) # returns the vector c(0.75,0.25) I've written such a function as follows: tapply.shingle <-

Hi, My problem is simple: since having updated the lattice package, I cannot load lattice anymore. If I type in the command 'library(lattice)' the loading fails with the following message: --- cut here --- Error in library.dynam(lib, package, package.lib) : shared library 'lattice' not found In addition: Warning messages: 1: In loadNamespace(package, c(which.lib.loc,