similar to: [Re] add reference lines (or grid) in background

Hi I'm getting an error message with the roots() function in the vars package. Even the example in the help file comes up with an error: > data(Canada) > var.2c <- VAR(Canada, p = 2, type = "const") > roots(var.2c) Error in UseMethod("roots") : no applicable method for 'roots' applied to an object of class "varest" The error is odd, for

Hi I'm trying to measure the effect of a policy intervention (Box and Tiao, 1975). This query has to do with the coding of the model rather than with the particulars of my dataset, so I'm not providing the actual dataset (or a simulated one) in this case, apart from some general description. The time series are of length n=34 (annual observations between 1977 and 2010). The policy

Hi, Imagine I have the following data frame: > a <- c(1,NA,3) > b <- c(2,NA,NA) > c <- data.frame(cbind(a,b)) > c a b 1 1 2 2 NA NA 3 3 NA I want to delete the second row. If I use na.omit, that would also affect the third row. I tried to use a loop and an ifelse clause with is.na to get R identify that row in which all records are missing, as opposed to the first

Dear all, Is there a simple way to add reference lines in background? I am trying with abline() or grid() but the lines, since they are executed after the plot function, are draw on top. How can I draw such lines beneath the main plot? Here is an example: x<-rnorm(100) boxplot(x) abline(h=c(-1,0,1)) grid(NA, 4, lwd = 2) regards, Luigi Marongiu, MSc [[alternative HTML

Dear all, I am a total beginner in R, so sorry if this is the wrong place. I am using R 2.10.1 on a Mac (Mac OS 10.6.2). I have this small dataset : growth sugar 75 C 72 C 73 C 61 F 67 F 64 F 62 S 63 S I have no problem reading the table, or getting the summary, but if I try boxplot(growth~sugar, ylab="growth", xlab="sugar", data=Dataset), I have the following error :

Hi, i would like to plot a few hundred .csv files. Each file contains one curve with x,y values to plot. I have been searching for "gnu r read many files qplot" and similar words. I found for loops that use assign to generate one variable containing a dataframe. When i uesed the classic "plot' command i could add the curves with something like for... {

Hi everybody, I am quite new to data imputation, but I would like to use the R package ' Amelia II: A Program for Missing Data '. However, its unclear to me how the input for amelia should look like: I have a data frame consisting of numerous coulmns, which represent different experimental conditions, whereby each column has 3 replicates. I want amelia to perform an imputation across

Hello, Is there a function in R by which one can run AR model with Bayesian Information Criteria (BIC)? To my knowledge, functions ar and ar.ols could select the order only by AIC. Thanks, Miao [[alternative HTML version deleted]]

Hi all, I want to test whether the MEAN of two different variables, (and different number of observations) are the same. I am trying to use the anova test but it doesn't seem to like that the number of observations are different: a=c(1:5) b=c(1:3) aov_test=aov(a~b) >>Error in model.frame.default(formula = a ~ b, drop.unused.levels = TRUE) : variable lengths differ (found for

Dear colleagues. I wonder if anybody knows about a procedure in R to calculate the Berger-Parker index and the Brillouin Diversity index. I searched the net but did not find anything good about it. Thanks a lot for any help Mayo, Elvira. [[alternative HTML version deleted]]

Hello, I am totally new in the field of time series analysis and forecasting and R. I read that R is a powerful tool for time series. Could anyone give me navigation what models of time series are availiable in R etc? -- View this message in context: http://r.789695.n4.nabble.com/Time-Series-Analysis-and-Forecasting-tp4646637.html Sent from the R help mailing list archive at Nabble.com.

Hello! I would like to make a scatterplot of my data, but the problem is that several points have the same x and y values and are represented as only one point. I wonder if there is a way to represent the data points, but with point size representing the frequency of the depicted x-y value? Thank you! -- View this message in context:

Dear R users and developers, I have a dataset containing 34 variables measured in a survey, which has some missing items. I would like to conduct a factor analysis of this data. I tested mi, Amelia, and MissForest as alternative packages in order to impute the missing data. I now have 5 separate datasets with the variables I am interested in factor analysing. In my reading of the package

noob here trying to make boxplots of some data i would like to separate the boxplots according to conditons of various levels for example: i have group:1 and 2, each group performed tests consisting of condition A,B,C,D side: left and right time: 1 to 10 I would like separate boxplots of the results (x) of the tests (numeric) for each group under each condition on each side over time. so far i

Hello, I would like to plot my data in a fashion similar to a boxplot, but plot the true data points without a box, just overlay lines of the summary generated by the boxplot. I have less than 10 observations within each group, so I think showing the actual data would be more effective than the box of the boxplot. I have been unable to find a way to do this. Here is example data: >

I've got a Linux server with over 1.5 Million files in a single directory, which as you will probably agree is not a good idea. My problem is that the management refuses to accept that it's a bad idea to have that many files in a single directory. To make the problem worse, the folder in question is exported via Samba and NFS! Has anyone come across any references/documents that I could

No objections from Martin or elsewhere, so I have now committed this. Thanks, Uwe. On Tue, 24 Feb 2009, Uwe Ligges wrote: > [CCing Martin and Brian who had both done most svn commits of boxplot.R so > far] > > > A very minor wishlist item that I should have already reported years ago: > > All the time when I need presentation/publication quality boxplots, I add >

************************************************************ Important: We would draw your attention to the notices at the bottom of this e-mail, particularly before opening and reviewing any file attachment(s). ************************************************************ Martin, Thank you for letting us know about posdefify. It does do exactly what the Rebonato paper recommended and

************************************************************ Important: We would draw your attention to the notices at the bottom of this e-mail, particularly before opening and reviewing any file attachment(s). ************************************************************ Here is some code we have used. a<-array(c(1,.9,.7,.9,1,.3,.7,.3,1),dim=c(3,3)) a s<-eigen(a)$vectors

Dear All, I have attempted to install the R-Excel Interface with poor results. The version of R is 1.6.1, the version of R-Excel is 1.0, the version of COM is 0.99. All of these were downloaded Monday, 6 January 2003. R (version 1.6.1) is installed on the machine and appears to run correctly on its own. The computer is a Compaq Pentium 4 machine with 128 meg of RAM running Windows 2000 and