similar to: the predict.lda function

Hi, I've noted that not all `predict' methods are hidden in the namespace: > methods("predict") [1] predict.ar* predict.Arima* [3] predict.arima0* predict.glm [5] predict.HoltWinters* predict.lm [7] predict.loess* predict.mlm [9] predict.nls* predict.poly [11] predict.ppr* predict.prcomp* [13]

Hi Every Body, I don't know why some regression functions have no related prediction function. For example lowess, loess, lokerns, lpridge, lpepa, and ksmooth. What could help? Is there any global or wrapper function so that can help? Regards, Amir Safari --------------------------------- [[alternative HTML version deleted]]

Dear All I am a biginner of R. I have difficulty with reading the code of a method. I am using the vars package to estimate a VAR model and I want to view the code of "predict" method for objects with class attribute "varest". I thougt I could just type the name "predict" without anything to display the code of the method as I often do with generic function.

Oh yes. The "load package" under the "packages menu" in the Windows version does that. To check I typed "library(assist)" after starting R. Same behavior, ssr is found, but others like predict.ssr, and plot.ssr, give a "not found" message. Thanks for the suggestion. Mike

Good Morning! May you help me? I need to understand the function predict. I need to understand the algorithm implemented, the calculations associated. Where can I find this information? Thank You!

Hello, I'm getting the following error message when I try to run 'tsdiag' on what seems to be a valid time series: > tsdiag(small) returns: [Error in tsdiag(small) : no applicable method for "tsdiag"] where small is a little test series where I have isolated this problem (the original has 30-years worth of daily data) When I print it (small), it looks like the

Hello, I just ran the predict.StructTS function using the AirPassengers data and got a ridiculous result. Here's what I ended up with: http://24.210.155.111/PredictWhat!.pdf Who wrote this? Am I seriously supposed to think this function would accurately predict a time series? -AnalogKid

I have a set of data (in a matrix). I spliced a column out and parsed it as.ts (time series). I then plotted the time series but I found that it was very noisy. I wanted to smooth it out. However, I am having some problems smoothing and plotting the smoothed version. > A <- as.ts(read.table(choose.files())) > x <- as.ts(A[,10]) > plot(x) > > plot(smooth(x))

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 there, Does anyone know how to predict ARMA? It doesn?t have either predict or forecast methods. I found couple of packages called fbasic and fseries at http://www.itp.phys.ethz.ch/econophysics/R/, which has ?arma.predict? in it, but it doesn?t seem to be working. Any help in this regard would be appreciated. Thanks in advance. Regards Skanda Kallur "Prediction is very difficult,

Using > methods("plot") [1] plot.Date* plot.HoltWinters* plot.POSIXct* [4] plot.POSIXlt* plot.TukeyHSD plot.acf* [7] plot.data.frame* plot.decomposed.ts* plot.default [10] plot.dendrogram* plot.density plot.ecdf [13] plot.factor* plot.formula* plot.hclust* [16] plot.histogram* plot.isoreg* plot.lm [19] plot.medpolish*

Dear All, R 1.9.1, Windows When copying and pasting a few lines from the 'predict.Arima' help, I get an error message: > data(lh) > predict(arima(lh, order = c(3,0,0)), n.ahead = 12) Error in eval(expr, envir, enclos) : Object "xreg" not found On the other hand, the following is OK: > data(lh) > predict(arima0(lh, order = c(3,0,0)), n.ahead = 12) $pred Time

I apologize in advance if I picked the wrong list to post this to. I have made an effort to find the answers to these questions on CRAN, but if they are there, I couldn't find them, and I was going to email the developer of StructTS directly but could not find who that is. I have 2 interrelated questions about StructTS 1. Where can I obtain the source code for StructTS if I wanted to

I am attempting to configure a dovecot 2.0 + exim + fetchmail (demon) combination on openSUSE 11.4, initially with the standard dovecot package (2.0.9-3.4-x86_64) and now with the later server:mail build service package (2.0.11-1.1-x86_64). I followed the instructions in http://wiki2.dovecot.org/HowTo/EximAndDovecotSASL but despite the changes in conf.d/10-master.conf the permissions on

Dear all, I'd like to discuss about a possible bug in function StructTS of stats package. It seems that the function returns wrong value of the log-likelihood, as the added constant to the relevant part of the log-likelihood is misspecified. Here is an simple example: > data(Nile) > fit <- StructTS(Nile, type = "level") > fit$loglik [1] -367.5194 When computing the

Dear all, I am applying the StructTS function in ts-package. For some time series the program terminates and the following error appears: Error in optim(init[mask], getLike, method = "L-BFGS-B", lower = rep(0, : L-BFGS-B needs finite values of fn Do someone know what do I have to adjust in the original time series to avoid this error? It works fine for some subsets of the

Hi all, I've experienced this problem using StructTS(x) when the *first* element of x is a NA (R:R1.7.0, os: w2ksp4). Please look at the following code: a=rep(1:7,10) library(ts) #this works StructTS(a) #this works x=a x[2]=NA StructTS(x) #this doesn't work x=a x[1]=NA StructTS(x) The last command returns this error "Error in optim(init[mask], getLike,

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?

Dear All, I am trying to fit a structural time series model using the StructTS function (package stats) with only 2 seasons (summer and winter). More than 2 seasons work fine but with 2 seasons I get this error: > fit <- StructTS(y.ts, type="BSM") Error in T[cbind(ind + 1L, ind)] <- 1 : subscript out of bounds I have looked at Prof. Ripley's 2002 RNews article but cannot