similar to: Initial covariance matrix in StructTS

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 was playing around with factor contrasts, and found the argument nmax on function factor. When using nmax=1, R froze completely, and I had to close it from task manager. After some debugging, I found that the problem is actually in unique-function, where the internal unique function is called: .Internal(unique(x, incomparables, fromLast, nmax)) More generally, it looks like

Dear all, I recently uploaded a new package Rlibeemd to CRAN (http://cran.r-project.org/web/packages/Rlibeemd/index.html, and in github: https://github.com/helske/Rlibeemd), and simple example in help files seems to throw an error on Solaris: Version: 1.3.1 Check: examples Result: ERROR Running examples in ?Rlibeemd-Ex.R? failed The error most likely occurred in: > ### Name:

Dear all: I want to modify the 'StructTS' function from the 'stats' package. First, I am writing a working copy of the original version and got some problems. I have two versions of the function plus the original one. The first version is the same code as the 'StructTS' function: StructTS.v1 <- function (x, type = c("level", "trend",

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,

I am having problems with forward prediction using the output of the Basic Structural Model from StructTS. The following snippet illustrates the problem: t_end <- 139 nahead <- 20 data(AirPassengers) ap <- log10(AirPassengers)-2 fit <- StructTS(ts(ap[1:t_end], freq=12), type="BSM") p <- stats:::predict.StructTS(fit, n.ahead=nahead) plot(1:t_end, ap[1:t_end],

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

Dear R-bugs I have been studying the StructTS function (in package 'stats') and functions supplied with it. I think I have found a few minor bugs in the documentation. I am referring to the version of StructTS supplied with the release R 2.0.0. Output from 'version' platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status

Dear List, I am thinking about ways in which I might analyse some stratigraphic data. The nature of the data series I have generates a number of issues: 1) The data I have in mind come from a sediment core sequence taken from the bottom of a lake. The sequence is sliced into a priori defined slices, in this case 0.2cm per slice. in this way a sequence of 0.2cm slices is produced for the

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

In the examples for StructTS -- ($RHOME)/library/stats/man/StructTS.Rd -- could par(mfrow = c(4, 1)) plot(log10(UKgas)) plot(cbind(fitted(fit), resids=resid(fit)), main = "UK gas consumption") become plot(log10(UKgas)) par(mfrow = c(4, 1)) plot(cbind(fitted(fit), resids=resid(fit)), main = "UK gas consumption") ## Note that par was moved down This makes the plot of UKgas

Dear R-list, I try to use StructTS for some univariate time series, with a type="BSM". I have two questions: 1-in the "fitted" value of the object, it seems to me that for getting back fitted values, one have to add level, slope and seasonal (or "level" contains already the two other ones) 2-i am really circumspect about the results of "KalmanForecast"

I have seen a couple of posts about this, but no solutions. The problem is fitting the Basic Structural Model (BSM) to the AirPassengers time series using StructTS. ?For this particular time series, if the length is reduced below 140 months, the BSM fits are bad.? The following illustrates the problem ap0 <- log10(AirPassengers) - 2ap<-ts(ap0[1:139], freq=12) ?# bad fits for length <

Hi, I used function StructTS some time ago to fit a structural model to a time series. Now with R 1.6.2-1 I repeated the analysis with the same series and I get the following error: Error in KalmanLike2(y, Z, -1) : invalid argument type I tried with other series and I get the same error; I checked the examples in the documentation and they work fine. I suspect I am missing something

Hello, Does anybody know if (and how) it is possible to obtain standard errors of estimated variances from StructTS? (R 2.10.0). Thank you in advance, Giovanni

Dear All I would like to use the StructTS function in the ts library to fit the ' BSM ' model. I have some, probably basics, questions about the model and about the function(s): 1) How can I check the statistical significance of the estimated parameters(variances)? 2) Is there some way to find what component "dominate" the series? 3) Is there a function to produce

Dear All I would like to use the StructTS function in the ts library to fit the ' BSM ' model. I have some, probably basics, questions about the model and about the function(s): 1) How can I check the statistical significance of the estimated parameters(variances)? 2) Is there some way to find what component "dominate" the series? 3) Is there a function to produce

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

Hello everybody, I'm investigating several time series with StructTS, and for one series I get zero variance for the measurement errors (local level model). I've read in Brian Dipleys article on time series in the R-Newsletter June/2002, that the airline passengers data set yields the same result. Having the variance of epsilon as zero, I would expect the residuals (y - Z a) to be

To whomever it may concern, I'm a young Industrial Engineer working on Senior Design at Georgia Tech and have found the StructTS method to be excellent for the training set for my forecasting project. There's only one problem: I don't actually understand what a Structural Time Series IS. I've looked up resources on it, and get that essentially you're dividing the Time