similar to: Where's iris?

Hi: I need to plot two time series in the same plot and they cover the same time range and have the same frequency. With RSiteSearch("multiple series plot") i found this post by Gabor Grothendieck: http://finzi.psych.upenn.edu/R/Rhelp02a/archive/42281.html Exactly what i need except for one detail. I want one series to be made of points and the other by a line. When I simply try: xts

De: Rubén Roa Enviado el: jueves, 17 de noviembre de 2011 9:53 Para: 'Users@admb-project.org' Asunto: Reporting a conflict between ADMB and Rtools on Windows systems Hi, I have to work under Windows, it's a company policy. I've just found that there is a conflict between tools used to build R packages (Rtools) and ADMB due to the need to put Rtools compiler's

Does anyone know of a package to make a map from GIS data, and/or would it be easier in one of the free GIS programs. I would like to make a map of the savannah river area with our sampling locations. thanks stephen -- Let's not spend our time and resources thinking about things that are so little or so large that all they really do for us is puff us up and make us feel like gods. We are

DeaR ComRades, This is a quote from a News article in Science's 11-February issue, about competitions to model data: "For Chris Raimondi, a search-engine expert based in Baltimore, Maryland, and winner of the HIV-treatment competition, the Kaggle contest motivated him to hone his skills in a newly learned computer language called R, which he used to encode the winning data model.

Hi R-listers, I have developed 47 GLM models with different combinations of interactions from 1 variable to 5 variables. I have manually made each model separately and put them into individual tables (organized by the number of variables) showing the AIC score. I want to compare all of these models. 1) What is the best way to compare various models with unique combinations and different number

Dear R Helpers, I have attempted "optim" function to solve a two-dimensional optimization problem. It took around 25 second to complete the procedure. However, I want to reduce the computation time: less than 7 second. Is there any optimization function in R which is very rapid? Best Regards, Wonsang ----- Wonsang You Leibniz Institute for Neurobiology -- View this message in

It's in the Venables & Ripley MASS (ed 3) book in the section on principal components. The context is as follows > ir <- rbind(iris[,,1], iris[,,2], iris[,,3]) > ir.species <- factor(c(rep("s",50),rep("c",50),rep("v",50))) (then they use brush(ir) which I guess is not an R function) and then > princomp(log(ir[1:4]),cor=T) (there is no [1:4]

Hi, I'm trying to figure out how to use lme() for analyzing a split-plot experiment. I've been looking at the examples from the 'R Book', those are nested but with only one factor at the whole-plot level, my test is 2^2 at the whole-plot level, with a single many level factor at the sub-plot level. My question is about properly specifying the random effects part of the model,

Dear all What is the easy way to drop a variable by using its name (and not its number)? Example: > data(iris) > head(iris) Sepal.Length Sepal.Width Petal.Length Petal.Width Species 1 5.1 3.5 1.4 0.2 setosa 2 4.9 3.0 1.4 0.2 setosa 3 4.7 3.2 1.3 0.2 setosa 4 4.6 3.1

I tried it using file and it seems to work for saving: > data(iris) > con <- file("clipboard","w") > save(iris,ascii=T,file=con) > close(con) > readLines("clipboard") ... lengthy output follows which seems correct ... but not for loading: > con <- file("clipboard","r") > load(con) Error in load(con) : loading from

Hi there, Suppose mu is constant, and error is normally distributed with mean 0 and fixed variance s. I need to find a statistics that: Y_i = mu + beta1* I1_i beta2*I2_i + beta3*I1_i*I2_i + +error, where I_i is 1 Y_i is from group A, and 0 if Y_i is from group B. It is large when beta1=beta2=0 It is small when beta1 and/or beta2 is not equal to 0 How can I find it by R? Thank you very much

DeaR ComRades, require(lattice) data <- data.frame(SP=sort(rep(as.factor(c('A','B','C','D','E')),12)), x=rpois(60,10), y=rep(c(rep(0,4),rep(10,4),rep(20,4)),5), z=rep(1:4,15)) xyplot(x~y|SP,data=data,groups=z,layout=c(2,3),pch=1:4,lty=1:4,col='black',type='b') How do I put a legend

Dear R list, I'm trying to calculate Mahalanobis distances for 'Species' of 'iris' data as obtained below: Squared Distance to Species From Species: Setosa Versicolor Virginica Setosa 0 89.86419 179.38471 Versicolor 89.86419 0 17.20107 Virginica 179.38471 17.20107 0 This distances above were obtained with proc

I am refering to a function call like this: >data(iris) >x <- svmlight(Species ~ ., data = iris) I tried to see the content of it by typing: > Species ~ . but it gives nothing. How can I see it's content ? - P.Dubois

DeaR ComRades, This is a quote from a News article in Science's 11-February issue, about competitions to model data: "For Chris Raimondi, a search-engine expert based in Baltimore, Maryland, and winner of the HIV-treatment competition, the Kaggle contest motivated him to hone his skills in a newly learned computer language called R, which he used to encode the winning data model.

Hola, ¿qué tal? Me han pedido que genere unos gráficos para insertarlos en un docuemento de Word. Yo suelo generar gráficas en formato PNG, pero me han dicho que "quedan mal". Y yo apenas sé nada de Word. ¿Cuál sería el formato gráfico más adecuado en este caso? Un saludo y muchas gracias, Carlos J. Gil Bellosta http://www.datanalytics.com

Is this a bug? I thought that "clipboard" could always be substituted for a filename when dealing with ASCII files. > data(iris) > save(iris,ascii=TRUE,file="clipboard") Error in file(file, "wb") : `mode' for the clipboard must be `r' or `w' Also this (where gclip is a utility found at unxutils.sourceforge.net that copies its standard input to the

Hi all! Does anyone knows if it exists a "symbolic math" package in R, that allows to compute derivatives, integrals, etc.? Does exist a freeware version of Maple? Cheers, Marco [[alternative HTML version deleted]]

I did the same transition, facing the same issues, so I settled for installation of RKWard. This pretty much reproduces my previous experience with Tinn-R and the Rgui. HTH Ruben ________________________________ Ruben H. Roa-Ureta, Ph. D. Senior Scientist, Center for Environment and Water, Marine Studies Section, King Fahd University of Petroleum and Minerals, KFUPM Box 1927, Dhahran 31261,

On 02/02/2010 6:20 AM, Dimitri Shvorob wrote: > Ruben Roa has kindly suggested using 'scipen' option - cf. > >> fixed notation will be preferred unless it is more than ???scipen??? digits >> wider. > > However, > > options(scipen = 50) > x = c(1e7, 2e7) > barplot(x) > > still does not produce the desired result. This is strange. I see what