similar to: scatterplot with wrong line offset

Dear r-users! I have another question regarding the dlm package and I would be very happy if someone could give me a hint! I am using the dlm package to get estimates for an endogenous rate of capacity utilization over time. The general form of a state space model is (1) b_t = G * b_t-1 + w_t w_t ~ N(0,W) (2) y_t= A' * x_t + H' * b_t + v_t v_t ~ N(0,V) (Hamilton 1984: 372) The

Hi, I am trying to figure how to use dlm with constant terms (possibly time-dependent) added to both equations y_t = c_t + F_t\theta_t + v_t \theta_t = d_t + G_t\theta_{t-1} + w_t, in the way that S-PLUS Finmetrics does? Is there any straightforward way to transform the above to the default setup? Thanks, Tsvetan -------------------------------------------------------- NOTICE: If received in

Dear all, I have a serious problem to solve my model. I study over exploitation of fish in the bay of biscay (france). I know only the level of catch and the fishing effort (see data below) by year. My model is composed by the following equations: * the growth function Gt(St) = r*St*(1-St/sbar) with Gt the growth of each period t r intrinsec growth of the stock sbar carriyng capacity of the

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)

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

Dear R-users, I am trying to estimate a state space model of the form (1) b_t = G * b_t-1 + w_t w_t ~ N(0,W) (2) y_t= A' * x_t + H' * b_t + v_t v_t ~ N(0,V) (Hamilton 1984: 372) In particular my estimation in state space form looks like (3) a3_t = 1 * a3_t-1 + w_t w_t ~ N(0,W) (4) g_t = (a1, a2) * (1, P_t)' + u_t * a3_t + v_t v_t ~ N(0,V) where g_t is the

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 all, I'm working on a term structure estimation using state-space modeling for 1, 2 and 3 factor models. When I started to read the functions on R, I got to the function ss on the library sspir. From what I understood this function is similar to SsfFit from S-PLUS. But for my models purpose there is something left to be desired. Its formulation follow these equations: *Y_t = F_t^T *

Dear R listers, I have a very annoying problem using nlm(). I want to find the minimizer of my target function, if written in \LaTeX is f(\mu1,\mu2,\sigma1,\sigma2) = \sum_i^n( w_ig_t(z_i) ), where g_t(z) is a pdf of bivariate normal distribution and z_i is my samples. I cannot get the estimation result generated by nlm(), and I got the following errors " Error in

I am having a problem with the package SSPIR. The code below illustrates it. I keep getting the message: "Error in y - f : non-conformable arrays." I tried to tweak the code below in many different ways, for example, substituting rbind for cbind, and sometimes I get a different error message, but I could not find a variation of this code that would work. Any help will be greatly

i am working on the project to analyze hedge fund performance, i would appreciate that if you guys could spare some time helping me out with the R code. Thanks. The senario is: i applied BOXPLOT() to plot the performance of all hedge funds with 7 strategies. And right now in this boxplot I need to plot the points of 30 individual hedge funds from my portfolio. And I applied POINTS() and

Hi There, I'm attempting to plot 10 values on a three-dimensional PCA with text labels next to each point. While i have no trouble doing this on 2D plots using the 'text' or 'textxy' function, I cannot find a function to do this on a 3D plot. I am using princomp for my PCA: >PCA<-princomp(eucdata, cor=TRUE) >PCA$scores [,1:3] # the three principal components i

Hi R folks, In my previous post I forgot to mention that I was new to R. I was really grateful for your quick help. I have two further questions: 1) In the graph of a regression line I would like to show one specific residual yi obs - yi pred (let's take the person whose residual is 76). How do I add a bracket to this vertical distance and name it? I'am getting stuck after the

Hi guys, I am doing a project to evaluate the 7 individual fund performance from a portfolio compared with all fund performance utilizing the same strategy. Lets say in total there are 10 strategies and in the portfolio there exist 5 strategies. First, i use the BOXPLOT() and SUBSET() to produce the box plot of all the 5 funds performance individually: #identify funds in strategies within

Hi R users, I would like to add year information to each point in a scatter plot. The following example shows what i mean: df=data.frame(year=c(2001,2002,2003),a=c(8,9,7),b=c(100,90,95)) df plot(b~a,data=df) Now, I would like to see which point belongs to 2001, 2002 and 2003 in the above graphic. Thank you very much, ozan [[alternative HTML version deleted]]

Hi everyone, i?m trying to graphically display distributions with r and i?m working with makrodata from the WVS. the command i?m using is > plot (Makrodata$v11, Makrodata$v12, xlab="Democracy Score Economist", > ylab= share religious people") i?m having an additional variable that identifies respectively labels the