similar to: Kalman-filtering

Dear all, I am pleased to announce the CRAN release of a new package called 'KFAS' - Kalman filter and smoother. The package KFAS contains functions of multivariate Kalman filter, smoother, simulation smoother and forecasting. It uses univariate approach algorithm (aka sequential processing), which is faster than normal method, and it also allows mean square prediction error matrix Ft to

Dear all, I am pleased to announce the CRAN release of a new package called 'KFAS' - Kalman filter and smoother. The package KFAS contains functions of multivariate Kalman filter, smoother, simulation smoother and forecasting. It uses univariate approach algorithm (aka sequential processing), which is faster than normal method, and it also allows mean square prediction error matrix Ft to

Dear all, I am currently trying to use the "dlm" package for Kalman filtering. My model is very simple: Y_t = F'_t Theta_t + v_t Theta_t = G_t Theta_t-1 + w_t v_t ~ N(0,V_t) = N(0,V) w_t ~ N(0,W_t) = N(0,W) Y_ t is a univariate time series (1x1) F_t is a vector of factor returns (Kx1) Theta_t is the state vector (Kx1) G_t is the identity matrix My first

Hi Does anyone know of any R code for estimating a *multivariate* state space model using a Kalman filter where the output matrix H(t) is time-varying but predictable (i.e. measurable w.r.t information at time t-1) in the observation equation y(t) = H(t) z(t) + R w(t)? [Here y(t) are the observations, z(t) is the state variable, w(t) the observation error and R R' the observation error

I am looking for a Kalman filter that can handle a control input. I thought that l.SS was suitable however, I can't get it to work, and wonder if I am not using the right function. What I want is a Kalman filter that accepts exogenous inputs where the input is found using the algebraic Ricatti equation solution to a penalty function. If K is the gain matrix then the exogenous input

Since notation for state-space models vary, I'll use the following convention: x(t) indicates the state vector, y(t) indicates the vector of observed quantities. State Transition Equation: x(t+1) = Fx(t) + v(t) Observation Equation: y(t) = Gx(t) + w(t) Cov[v(t)] = V Cov[w(t)] = W I've found myself in a situation where I will have V = s%*%t(s)*k^2, with s a vector the same length as the

Hello again, A different but related question to my last one: Does anyone know if one can easily estimate state-space models using ML and the kalman filter using R? I would be especially interested in a relatively flexible function that would allow for estimation of hyperparameters, or could be made to do so. Thanks Michael J. Roberts Resource Economics Division, PMT USDA-ERS 202-654-5557

Hi My name is Vladimir Samaj. I am a student of Univerzity of Zilina. I am trying to implement Kalman Filter into my school work. I have some problems with understanding of R version of Kalman Filter in package stats( functions KalmanLike, KalmanRun, KalmanSmooth,KalmanForecast). 1) Can you tell me how are you seting the initial values of state vector in Kalman Filter? Are you using some method?

Hi all! I'm sure this must have been asked many times before but here goes anyway. I'm looking for a kalman filter in R for ar(i)ma time series. I'm sure there must be one around but it does not seem to be in either ts or tseries packages? Any suggestions welcome. Thanks Gerard Keogh The information in this email, and any attachments transmitted with it, are confidential and are

Dear all: Greetings! I am trying to replicate a simple state space model in R, using the package 'dlm'. This model has two observation equations and three state equations. Each observation equation represents structural shocks based on SVAR for country i, where i=1 to 2. The structural shocks (S1 and S2) are to be decomposed into common (sv1) and country-specific (sv2, sv3) shocks. We

Hi, Following convention below: y(t) = Ax(t)+Bu(t)+eps(t) # observation eq x(t) = Cx(t-1)+Du(t)+eta(t) # state eq I modified the following routine (which I copied from: http://www.stat.pitt.edu/stoffer/tsa2/Rcode/Kall.R) to accommodate u(t), an exogenous input to the system. for (i in 2:N){ xp[[i]]=C%*%xf[[i-1]] Pp[[i]]=C%*%Pf[[i-1]]%*%t(C)+Q siginv=A[[i]]%*%Pp[[i]]%*%t(A[[i]])+R

help.search("Kalman") says to look at help(KalmanLike, package=ts). Andy -----Original Message----- From: Mohamed A. Kerasha [mailto:mohamed at engr.uconn.edu] Sent: Tuesday, November 19, 2002 9:27 AM To: r-help at stat.math.ethz.ch Subject: [R] Kalman Filter Hi all, Does any one know if there is Kalman Filter code or library in R. Thanks, Mohamed.

Hi All, Currently I'm using DSE package for Kalman Filtering. I have a dataset of one dependent variable and seven other independent variables. I'm confused at one point. How to declare the input-output series using TSdata command. Because the given example at page 37 showing some error. rain <- matrix(rnorm(86*17), 86,17) radar <- matrix(rnorm(86*5), 86,5) mydata <-

Hello, I would like use Kalman filter for estimating parameters of a stochastic model. I have developed the state space model but I don’t know the correct way use Kalman filter for parameter estimation. Has anybody experience in work with Kalman filter in R. I don’t know the correct function. Maybe it is - KalmanLike; but what is the correct Input? - tsmooth? -

Dear R Users, I am new to state-space modeling. I am using SSPIR package for Kalman Filter. I have a data set containing one dependent variable and 7 independent variables with 250 data points. I want to use Kalman Filter for forecast the future values of the dependent variable using a multiple regression framework. I have used ssm function to produce the state space (SS)

There is now a preliminary version of a time series package in the R-devel snapshots, and we would welcome feedback on it. It is based in part on the packages bats (Martyn Plummer) and tseries (Adrian Trapletti) and in part on code I had or have written. (Thanks for the contributions, Martyn and Adrian!) Some of the existing ts code has been changed, for example to plot multiple time series, so

I have built a Kalman Filter model for flu forecasting as shown below. Y - Target Variable X1 - Predictor1 X2 - Predictor2 While forecasting into the future, I will NOT have data for all three variables. So, I am predicting X1 and X2 using two Kalman filters. The code is below x1.model <- dlmModSeas(52) + dlmModPoly(1, dV=5, dW=10) x2.model <- dlmModSeas(52) + dlmModPoly(1, dV=10,

1. The function "KalmanLike" seems to change its inputs AND PREVIOUSLY MADE copies of the inputs. Consider the following (using R 2.1.0 patched under Windows XP): > Fig2.1 <- StructTS(x=Nile, type="level") > unlist(Fig2.1$model0[2:3]) a P 1120 286379470 > tst2 <- tst <- Fig2.1$model0 > tst23 <- tst[2:3] > tst23u <-

Dear R Helpers, Currently I'm working with the ts package of R and created a TimeSerie from pixels extracted from satellite imagery(S10 NDVI data, 10 daily composites). I'm trying to decompose this signal in different signals (seasonal and trend). When testing out the STL method is says => Only univariate timeseries are allowed, but the current Timeserie I'm using is univariate!

Hi, I am studying Kalman Filter and it seems to be difficult for me to apply the filter on a simple ARMA. It is easy to construct the state-space model, for instance: dlmModARMA(ar=c(0.4,-0.2),ma=c(0.2,-0.1, sigma2=1) but applying the dlmFilter on it, it doesn't work... I don't know if my problem is clear but if anyone has already worked on Kalman filter, it could be great to advise me!