similar to: summary on predict with arima0

I've got a data frame that looks like this: subject foo bar 2 1.7 3.2 2 2.3 4.1 3 7.6 2.3 3 7.1 3.3 3 7.3 2.3 3 7.4 1.3 5 6.2 6.1 5 3.4 6.9 ... That is, I've got multiple rows per subject. I need to compute summaries within categories where the subject has the same number

I am generating some graphs (on GNU/Linux using R 1.5.x) for a PowerPoint presentation my boss is giving next week. He just tried cut/pasted a plot off of a PDF file I had created, but he complained that the fonts were too small and fuzzy, and that he couldn't change the background. I have been playing with png(bg=transparent), but I am not sure what appropriate height,width parameters and

I am writing some R code that generates some figures. I am primarily developing on GNU/Linux, but will ultimately run the code on Windows in order to get WMF files I can integrate in an existing document. The win.metafile() function, of course, doesn't exist in the GNU/Linux version of R. I am contemplating writing an abstraction layer that runs the right thing, but I need to know how the

I am getting ready to load a bunch of data into R. The data is all numeric, but some of the numbers are integer codes representing non-numeric semantics. What is the best way to "fix" the data frames so that these compenents are recognized as "factors" or "ordered", as appropriate? Can I "assign" to some attribute of the data frame component, like the

Dear R People: I'm trying to use the predict command on an arima0 object. I do the following: xm.arma <- arima0(xm2,order=c(1,0,1)) predict(xm.arma,n.ahead=2) and I get the message: Error in round(x, digits) : Non-numeric argument to mathematical function Any ideas what the problem might be, please? R version 1 4 1 on Windows. Thanks in advance! Sincerely, Erin Hodgess Associate

Hi, I'm trying to go through the examples for function arima0() in ts package, i.e, >data(lh) >arima0(lh, order = c(1,0,0)) each time the call to arima0() causes a segmentation fault. I checked the earlier version (1.1.1) of R, the function arima0 works fine. Tracing the call indicates that the function "setup_starma" (in pacf.c under ts) interprets the addresses of the

In the 'RWeaveLatexSetup' function in the 'tools' package in R 1.6.0 , the code below does not work (pdflatex can't find 'Sweave.sty'): if (stylepath) styfile <- file.path(R.home(), "share", "texmf", "Sweave") else styfile <- "Sweave" In R 1.5.1, that was: if(stylepath) styfile <-

Hi, I work on order estimation for autoregressive processes and after some inconsistencies cropped up I implemented the AIC criterion myself. Its results do not match the implementation in R and there are a few things I can not understand even after reading the source code of R. I used the function called arima0 (with empty "ma" coeficient vector), and I do not understand how some of

[accidentally sent to owner-r-help -- please do NOT! it's "r-help" !] ------- start of forwarded message ------- From: Marcos_Sanches at gallup.com To: owner-r-help at stat.math.ethz.ch Subject: How to get residuasl with arima0? Date: Wed, 19 Sep 2001 15:19:07 -0300 I know this is a basic question, but I've never used the 'ts' package and I'm having some

Hi, I just had a very quick look at the StatDataML proposal --- nice work! At the risk of showing my ignorance, I want to mention my first impressions. My first impression is that defining datasets in terms of arrays and list is a bit too high a level. What about simpler vectors, scalars? (I know that R/S don't have scalars, but other systems/applications do.) Can we think of a core

Dear all, this may be a stupid question but... The underlying model in procedure arima0 is X[t] = a[1]X[t-1] + ... + a[p]X[t-p] + e[t] + b[1]e[t-1] + ... +b[q]e[t-q] Is it possible to get an estimate of e for every point t, t-1 etc. or at least an estimate of the variance of e? Thanks a lot in advance for any hints Kai Arzheimer

What is the difference between arima and arima0?

Full_Name: Arto Luoma Version: 1.1.0 OS: Windows 98 Submission from: (NULL) (153.1.53.119) In arima0 it is possible to specify whether the mean of the original series is included in the model or not. However, it is not possible to specify whether the mean of the differenced series is included. It seems that it is not included. However, if differencing is used to eliminate trend, the mean of the

Dear R People: Suppose I want to estimate the parameters of the following AR model: (1 - phi_1 B - phi_2 B^2 - phi_9 B^9) x_t = a_t and I want to use the arima0 command from the ts library. How would I use the order subcommand, please? R Version 1.4.1 for Windows. Thanks! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston -

does the arima/arima0 function use the state space form of the model equation even when fitting with the "CSS"-method? regards Christoph [[alternative HTML version deleted]]

I'd like to generate some plots like you'd see on the old "normal probability graph paper", like the first plot in: <http://www.itl.nist.gov/div898/handbook/eda/section3/normprpl.htm> except the horizontal scale would have 1%, 5%, 25%, 50%, 75%, 95%, 99%, or similar quantiles, with associated tick/grid lines. [still hunting around for a good example...] something like

Full_Name: Ahmad Abu Hammour Version: rw0651 OS: windows 95 Submission from: (NULL) (63.23.128.44) Although I know that "ts package" is preliminary, I wanted to compare the results from R and SPSS. I ran ARIMA(2,1,2) in both softwares. I got NaN in standard errors of coefficients from R and real figures from SPSS. I changed "delta" in R to match that used by SPSS, I received

Dear R People: I have a data frame, xm1, which has 12 rows and 4 columns. If I put is xm1[,-4], I get all rows, and columns 1 - 3, which is as it should be. Now, is there a way to use the names of the columns to omit them, please? Thanks so much in advance! Sincerely, Erin -- Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown

Hello! Those of you, who use Sweave a lot, will probably find my shell script usable. You can get it at: http://www.bfro.uni-lj.si/MR/ggorjan/programs/shell/Sweave.sh No warranty, however don't hesitate to contact me if you find an error or have a patch! -- Lep pozdrav / With regards, Gregor GORJANC --------------------------------------------------------------- University of

Hello. I was wondering how one can find the exact values of the range of an axis within a plot. In xlispstat it was (send plot :range 0). Thanks much, Jason