similar to: indexing within panels in xyplot

Hello, I'm trying to build a frequency table for a vector, broken down by the combination of factors. For further analyses, I need to have the result arranged in a new data frame, with the upper limit of the histogram's breaks, the per bin count, and the factors identifying each record. My problem is including the latter: ---<---------------cut

Dear all, I wish to generate a lattice boxplot which skips an empty cell in a design. I have trawled r-help, scruitinized xyplot(lattice) help page, and merrily reproduced examples of using skip from a couple of previous r-help queries and the example given in Pinheiro & Bates. But I must be missing something... Here's an example (running R 1.9.1 on Win2k): # generate some data df1

Dear Latticers, I want to position two trellis plots of different forms side-by-side. The plot types are slightly different, aspect=1 required, but panels should look the same in both plots. Current workaround uses a guessed factor. Any way to improve this? Dieter library(lattice) n1 = 20 # I cannot rbind df1 and df2, because the x-dimensions are # different and must be scaled individually

I am still getting these in my maillog: Jul 13 07:39:38 pop5.cfu.net dovecot: IMAP(breu at cfu.net): Corrupted index cache file /var/spool/mail/filer/storage//cfu.net/b/r/breu/dovecot.index.cache: invalid record size Jul 13 09:35:03 pop5.cfu.net dovecot: IMAP(breu at cfu.net): Corrupted index cache file /var/spool/mail/filer/storage//cfu.net/b/r/breu/dovecot.index.cache: invalid record size

Hi, I have a question about how I should report the results for a linear mixed effects model where the model includes as predictors three factors (facA, facB and facC), one of which (facA) interacts with the other two. facA and facB have two levels and facC has 3 levels. There are also several other continuous predictors (e.g. varA, varB, varC). My mixed model is specified with the following

Hi people, I have a text file like this one posted: snp_id gene chromosome distance_from_gene_center position pop1 pop2 pop3 pop4 pop5 pop6 pop7 rs2129081 RAPT2 3 -129993 "upstream" 0.439009 1.169210 NA 0.233020 0.093042 NA -0.902596 rs1202698 RAPT2 3 -128695 "upstream" NA

Hi all i need to put on one graph 2 functions who's x axis is the same and y not. I mean on horizontal the time, and on vertical left: pressure, on vertical right: rpm of a motor, is R able to do that? i've found this that i could adapt maybe (i don't need time series really?) :/ : (http://tolstoy.newcastle.edu.au/R/help/04/03/1456.html) ## ## Description: A simple function which

I'm a newbie, so please forgive me: I have a data frame that includes two long columns of x and y coordinate values thusly: ... ... ... x y ... ... ... ... ... 115 127 ... ... ... ... ... 120 129 ... ... ... ... ... 0 0 ... ... ... ... ... 135 162 .... ... ect... i need to ask R tell me how many x-y pairs contain coords that are

#See David Williams' book "Weighing the odds", p286 y <- c(1.21, 0.51, 0.14, 1.62, -0.8, 0.72, -1.71, 0.84, 0.02, -0.12) ybar <- mean(y) ylength <- length(y) ybarv <- rep(ybar, ylength) x <- 1:ylength plot(x,y,asp=1,xlab="position",ylab="ybar",type="n",ylim=c(-1,1)) segments(x[1], ybar, x[ylength], ybar) segments(x,ybarv,x,y)

Dear all, Within a larger script I use the following code to produce barplots for all variables contained in a dataframe "type3m": for(xn in names(type3m)) barplot(table(type3m[[xn]]), main = paste ("Barplot of", xn), width=0.75,legend.text=barnamesm, xlim=c(0,4),ylim=c(0,400)) The plots are produced and printed into a *.ps file. However, R gives an error

Hi there, I found one reference to add a reactangle behind a plot using plot(...,add=T), I tried this but didn't obtain the desired result. If a I have the following code: plot(x,y) rect(xleft, ybottom, xright,ytop,col='green) The rectangle appear on top of the plot. Any help will be greatly appreciated, Judith

Good afternoon everybody, I'm quite new in functions optimization on R and, whereas I've read lot's of function descriptions, I'm not sure of the correct settings for function like "optimx" and "nlminb". I'd like to minimize my parameters and the loglikelihood result of the function. My parameters are a mean distance of dispersion and a proportion of

Dear all, I am struggling to modify the axis labels/ticks in a panel provided to xyplot. To begin with, I do not know the equivalent of the xaxt="n" directive for panels that would set the stage for no default x axis being drawn. My goal is to draw ticks and custom formatted labels at certain hours of the week. When I execute the code below, I get an error message in the plot window that

Dear all, is there a possibility to display pictures (jpegs, pngs, ...) in R? Maybe the following (hypothetical) code makes it clearer what I have in mind: plot(x=0:100, y=0:100, type="n") rect(xleft=25, xright=75, ybottom=30, ytop=95, fill="apicture.png") Thank you, Roland [[alternative HTML version deleted]]

Sun, 14 Mar 2004, Jan Verbesselt wrote: > Dear R specialists, > > I have two time series in a data.frame and want to plot them in the same > plot(), with the left axis scaled to time series 1 (-700,0) and the > right axis scaled to time series 2 (-0.2, 0.4). > > plot(timeserie1) > lines(timeserie2, col=c(2)) => this one should be scaled differently > with a new

The following gives two functions for producing distribution graphs: distribution-graph produces a single graph, and multiple.distribution.graph produces a number of graphs side by side. Regards, Tore Wentzel-Larsen statistician Centre for Clinical research Armauer Hansen house Haukeland University Hospital N-5021 Bergen tlf +47 55 97 55 39 (a) faks +47 55 97 60 88 (a) email

Still trying to get a handle on c, cat, paste, and lists. Here is my current problem. I have a good number of vectors that I want to selectively combine into a dataframe. The exact name of the vectors changes in a looping arrangement but can easily be grepped. The problem is: how can I use a list produced by a grep in a data.frame command: i.e. > ls() [1] "a"

Hi all, I am looking for a clever way to create the following graph using R: I got information on the shares of some subgroups over time (summing up to 1 in each year). The graph I want to create should display the development of the individual shares over time by shading rectangulars for each share in a different color. Is there a clever of doing this? Thanks for your help Stefan

I'm doing a complete upgrade of a server which I was using uw-imap for the imap/pop4 daemons.(This was an old RH 7.2 system and has been upgraded to CentOS 4.4 64bit) I notice there isn't an rpm at any of the repositories for this. I do pop before smpt for the remote users I have and also use procmail for calling spamc from SpamAssassin. Is there some place I can get uw's imap or

Hi, I'm trying to construct a legend which has four lines of text and associated symbols. The first two symbols need to be normal lines which vary only in colour. The second two symbols should have filled boxes. How do I suppress the fill boxes in the first two lines? J.