similar to: sorting data

Dear useRs, I'm trying to specify the colour of a factor with ggplot2. The example below gets me close to what I want, but it's missing a legend. Any ideas? Thanks, Thierry library(ggplot2) dataset <- data.frame(x = rnorm(40), y = runif(40), z = gl(4, 10, labels = LETTERS[1:4])) ggplot(data = dataset, aes(x = x, y = y, group = z)) + geom_point(colour = c("red",

Hi I have the following code which I would like to simplify. Id does linear regressions and returns the r-squares, and the coefficients. It runs slow, as it is doing the regressions for each - is it possible to get the values in a dataframe which looks as follow: expert | xx | seeds | r.squared | slope | intercept Thanks in advance, Rainer library(reshape) rsqs <- as.data.frame(

Hi, I am trying to perform Maximum Likelihood estimation of a Multivariate model (2 independent variables + intercept) with autocorrelated errors of 1st order (ar(1)). Does R have a function for that? I could only find an univariate option (ar.mle function) and when writing my own I find that it is pretty memory-consuming (and sometimes wrong) so there must be a better way. Thanks, KB

If you tell me how to update R itself automatically, I will go for your advice. I am not aware of any method to do it... Bye Rick "ONKELINX, Thierry"

Se ha borrado un texto insertado con un juego de caracteres sin especificar... Nombre: no disponible Url: https://stat.ethz.ch/pipermail/r-help/attachments/20080212/16ba8b97/attachment.pl

Dear R list members I am trying to improve a boxplot with 2 data sets. I run somethinkg like boxplot(data1 ~ month, add=F, col = "red", ...) boxplot(data2 ~ month, add=T, col = "blue", ...) The problem is that the data from February are missing for data2, so R think that must take little more space between the data classes in data 2 and then both data gropus are not aligned.

Dear useRs, Lately I noticed a strange warning in the summary of a lm-object. Any idea what this warning is about? I'm using R 2.5.1 on Win XP pro. > x <- rnorm(100) > y <- rnorm(100) > summary(lm(y~x)) Call: lm(formula = y ~ x) Residuals: Min 1Q Median 3Q Max -1,77809 -0,68438 -0,04409 0,63891 2,30863 Coefficients: Estimate Std.

Hi, I'm scratching my head as to why I can't use the subset() command to remove one line of data from a data frame. There is just one row (out of 45840) that I'd like to remove and it can be identified using.... > dim(raw.all.clean) [1] 45840 10 > subset(raw.all.clean, Height.1 == 0 & Height.2 == 0) Sample.Name Well SNP Allele.1 Allele.2 Size.1 Size.2

Hello My colleagues and I have recently established a large database (40 tables each with greater than 15 variables) in Microsoft's SQL Server 2000. Currently we are accessing this database via SQL client running an Windows XP. Our objectives are many fold including running SQL applications, outputting results to ARC/INFO IMS, production of summarizing tables - graphs and web interfaces for

I need the most straightforward way to build semivariograms within R i am currently using s-gems software but i would like to even test R Thanks D [[alternative HTML version deleted]]

Hello, I use newly R! I'd like to plot several data set together in one output window! How can I do that? Best regards Hadi --------------------------------- Never miss an email again! [[alternative HTML version deleted]]

Hi all, I am wondering if there is a function available in R for simulating spatial distribution of objects (plants in this case) in order to simulate sampling of a population . Specifically, I would like to be able to change the spatial correlation of individuals. I don't want to reinvent the wheel if it already exists. Thanks, Wade

Hi, my problem belongs to the basic ones. I want to get cumulated sum over the matrix columns by one command (if such exists). Ordinary R's cumsum(x) when x is: [,1] [,2] [1,] 1 5 [2,] 2 6 [3,] 3 7 [4,] 4 8 yields: 1 3 6 10 15 21 28 36 I want: [,1] [,2] [1,] 1 5 [2,] 3 11 [3,] 6 18 [4,] 10 26 Is there any command to do so?? best,

Hi R-users, I have a list a <- list(one=matrix(rnorm(20), 5, 4), two=matrix(rnorm(20, 3, 0.5),5,4)) How to add rowSums (calculated using lapply) to corresponding matrix in this list lapply(a, function(x) rowSums(x)) ?? -Lauri

Hello, I have a question about "qplot": How can I add another line to the same plot ? (like function "lines" in "plot"). Thank You, Sigalit. [[alternative HTML version deleted]]

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible Url : https://stat.ethz.ch/pipermail/r-help/attachments/20080305/f2c52eea/attachment.pl

I am confused about why the following occurs: > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor 4.0

Hello, I'm wondering if there are any packages/functions that can perform conditional sequential gaussian simulation. I'm following an article written by Grunwald, Reddy, Prenger and Fisher 2007. Modeling of the spatial variability of biogeochemical soil properties in a freshwater ecosystem. Ecological Modelling 201: 521 - 535, and would like to explore this methodology.

Hello again With the poylgon function it's possible to plot shaded areas under a curve. But somehow it connects the start and the endpoint of a line and fills whats between them. I would actually like to set the boarders of the shading by two min and max curves, but I failed. i tried stl like: polygon(min, max, border=0), where min and max are my curves plotet by plot(max) lines(min)

Hi, dear R developers, I've got a vector of monthly volatilities and i would like to get the position of the highest volatility of the vector without computing a loop. Is there a function that could give me such a result ? a<-c(1,2,4,100,3) the highest value is the fourth of the vector. how can i get "4" without a loop going through the vector ? Thanks ! Benoit.