similar to: plotting 1000 simulations, error message: plot.new has not been called yet

Hi, I want to make model selection with regsubsets. My code is: a<-regsubsets(Gesamt ~ CommunistSocialist + CountrySize + GNI + Lifeexp + Schoolyears + ExpMilitary + Mortality + PopPoverty + PopTotal + ExpEdu + ExpHealth, data=olympiadaten, nbest=2) summary(a) plot(a,scale="adjr2") (output attached) The problem is now, that I want to fit the best model again "manually"

Hi, I have a scatter plot of the variables GNI and Lifeexp (Gross National Income and Life Expectancy, both metric). So I plotted them and I want to add a regression line and a lowess line. I use lowess and not loess because I have missing values. My code: plot(GNI,Lifeexp) abline(lm(Lifeexp~GNI), col="red") y.loess<-loess(Lifeexp~GNI,na. action = na.exclude)

I was wondering if anyone out there has created a series of images (e.g. jpegs) using R, then animated them on a web page. I have looked around a bit and it seems fairly complicated. One approach that seems sensible is to use a java script that displays the list of images. I have only used html on web pages before so I don't know how to use java... Ideally the final result is easily ported to

Dear R-helpers, Is it possible to create a type of 'movie' in R based on the output of several figures (e.g., jpegs) via the plot() function. I obtained dynamic results with the plotting function and would like to save these as a movie (e.g., avi or other formats)? Regards, Jan _________________________________________________________________ Ir. Jan Verbesselt Research

Dear colleagues, Is there any way of saving an animated plot with R? For instance, any format that could be read by Microsoft windows media or whatever. Cheers, Jesus

Dear All, I'd like to know if it is possible to create animations with R. To be specific, I attach a code I am using for my research to plot some analytical results in 3D using the lattice package. It is not necessary to go through the code. Simply, it plots some 3D density profiles at two different times selected by the user. I wonder if it is possible to use the data generated for different

Hi there, I am running Monte Carlo Simulations in R using ordinary "while (condition)" loops. Since the number of iterations is something like 100.000 and within each iteration a given subsample is extended sequentially it takes hours to run the simulation. Does anyone know if there is either a way to avoid using loops in Monte Carlo Simulations or how to include possible faster

Hello all I am playing with R for to make a animated GIF. any suggestions, improvements are welcome :-) case somebody could help me, i thanks! Cleber N. Borges ( klebyn ) my objective: (steps TODO) ------------------- 1) to save PNG files; -----> i don't know the best way to make this; 2) transform the PNG files into GIF files (easy! no problem! ... i think ...) 3)

I am taking a basic statistics course this summer, and while the majority of the class is using a statistical package that came with the book, I am doing everything in R for practical reasons. Forgive me if there is documentation/instruction easily available on this topic, but Google and forums searches left me with nothing. We are learning about randomness at the moment and are required to

Hi - I am doing some monte carlo simulations comparing bayesian (using Plummer's jags) and maximum likelihood (using lmer from package lme4 by Bates et al). I would like to know if there is a way I can flag nonconvergence and exceptions. Currently the simulations just stop and the output reads things like: Error in optim(.Call("lmer_coef", x, 2, PACKAGE = "Matrix"), fn,

Greetings List, I am conducting some large simulations using R. As a result, I get many plots but I'm having some trouble with including some of them in a Microsoft Word document. Can any one tell me the easiest method of having copies of the R-graphs in the Word documents? Best regards Mahmoud [[alternative HTML version deleted]]

Hi to all the people (again), I'm doing some simulations with the memisc package of an own function, but I've the problem that I didn't know how to read the result of such simulations. My function is: > Torre<-function(a1,N1,a2,N2) + {Etorre<-(a1*N1)/(1+a1*N1) + Efuera<-(a2*N2)/(1+a2*N2) + if(Etorre>Efuera)Subir=TRUE + if (Etorre<Efuera)Subir=FALSE +

I want to create a 3d plot with densities. I use the function density to first create a 2d dimensional plot for specific x values, the function then creates the density and puts them into a y variable. Now I have a second set of x values and put it again into the density function and I get a second set of y variables and so on.... I want to put those sets into a 3d plot, I hope you know what I

I need help with these R simulations _______________________________________ The generation of pseudo-random numbers with uniform distribution [0.1] can be done through a method called Congruencial Committee, which consists the use of a recursive relationship of the type: xi +1 = (a.xi + c) mod m being a, c, and the positive integers m, c <m. Discuss in light of possible values of a, c, my

R-help, I'm using a function whose end result is a trellis plot. When I call the function I get sometimes the following message: "Note: The default device has been opened to honour attempt to modify trellis settings " leading up to any plot whatsoever. I call the package 'lattice' within the function and 'detach' after plotting. I have to close the graphics

Dear R-users, I am trying to perform some simulations from a model defined by ordinary differential equations. I would be grateful if someone could indicate some functions/packages/examples I could look at. Thank you in advance. Sebastien

Hello, I have estimated an ARIMA model and I would like to make simulations from this estimated model 1,5 and 10 steps ahead. Does anybody know how to do this? Thank You Felipe Parra [[alternative HTML version deleted]]

Hi all, I am a biostatistician and I have developed my own ranking system for clinical data. I would like to test the efficiency of it w.r.t. to other ranking systems. I would like to simulate the data and after assigning ranks to my observed scores(after neglecting dropouts), observe the type I error. If I want to do a Kruuskal Wallis type of test, what test statistic should I use to test for a

Dear R-help, I am trying to write a function to simulate datasets of size n which contain two time-to-event outcome variables with associated 'Event'/'Censored' indicator variables (flag1 and flag2 respectively). One of these indicator variables needs to be dependent on the other, so I am creating the first and trying to use this to create the second using an if/else statement.

Hello. I was exploring the R software package and received the error message "plot.new has not been called yet" when following "Appendix A A Sample Session" of R-intro.pdf. I searched the message archives and found no similiar report to mine. I am using R on CentOS and I am using the latest from R-project, version 2.2.0. Here is my .Rhistory file recorded when following the