similar to: Help with STL function to decompose

Hi everybody. I'm using BFAST program and I have some data to be analized (named "META.csv" for example). When I call BFAST in the following way: meta <- read.csv("C:/META.csv") meta.ts<-ts(meta, frequency=12, start = c(2005,9)) dis<-frequency(meta.ts) / length(meta.ts) meta.sc<-bfast(meta.ts, h=dis, max.iter=1) meta.sc$output The last statement has the

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi I would like to add a line to the plot of the data panel of a stl decomposed time series. Example: plot(stl(nottem, "per")) plots a 4(or is t 8?) panel graph. I thought that I might be able to use mfg to plot in the top (left?) panel ("data"), but it is not working. Any help appreciated, Rainer - -- Rainer M. Krug,

Hello, Guys: I'm from China, my English is poor and I'm new to R. The first message I sent to R help meets some problems, so I send again. Hope that I can get useful suggestions from you warm-hearted guys. Thanks. I builded a multiplicative seasonal ARMA model to a series named "cDownRange". And the order is (1,1)*(0,1)45 The regular AR=1; regular MA=1; seasonal AR=0; seasonal

The decompose() function truncates the seasonal component unnecessarily. I've modified the function to fix this problem, and also added the original data to the object returned (to enable better plotting). I've also modified the plot.decomposed.ts() function so that it plots the original data in the top panel rather than the reconstructed data. The difference between the two is that the

Hi I'm just stepping through the decompose() function, in "stats". Does this contained line of code not work if you have a time series ending "unevenly" (i.e., middle of the year), or am I missing something? season <- na.omit(c(as.numeric(window(season, start(x) + c(1, 0), end(x))), as.numeric(window(season, start(x), start(x) + c(0, f))))) The line seems to

Hello , I have estimated the following model, a sarima: p=9 d=1 q=2 P=0 D=1 Q=1 S=12 In R 2.12.2 Call: arima(x = xdata, order = c(p, d, q), seasonal = list(order = c(P, D, Q), period = S), optim.control = list(reltol = tol)) Coefficients: ar1 ar2 ar3 ar4 ar5 ar6 ar7 ar8 ar9 0.3152 0.8762 -0.4413 0.0152 0.1500 0.0001 -0.0413 -0.1811

I'm trying to use the following command. arima (x, order = c(p,d,q), seasonal =list(order=c(P,D,Q), period=s) How can I write an estimated seasonal ARIMA model from the outputs. To be specifically, which sign to use? I know R uses a different signs from S plus. Is it correct that the model is: (1-ar1*B-ar2*B^2-...)(1-sar1*B^s-sar2*B^2s-....)(1-B)^d(1-B^s)^D

Drug US1 US2 Aptecha Celebrex 235.54 269.99 121.02 Detrol LA 157.99 190.99 55.3 Flomax 166.00 190.99 93.45 Lipitor 174.99 200.99 137.7 Novaldex 108.6 129.99 22.48 Norvasc 186.66 203.99 161.93 Plavix 107.99 106.99 64.53 Prevacid 117.39 134.99 59.83 Prilosec 115.99 126.99 57.75 Zyrtec 181.1 200.99 58.79 US1=c(235.54,157.99,166,174.99,108.60,186.66,107.99,117.39,115.99,181.10)

Hi R-help, I'd like to extract the seasonal component of a short timeseries, and was hoping to use stats::decompose. I don't want to decompose the 'trend' component so I thought I should call decompose(x,filter=0). I think I've either misunderstood the filter argument or come upon a bug/feature in decompose. # EXAMPLE

Greetings, For the following quarterly data I did a classical decomposition by hand in a spreadsheet and got reasonably similar results using Minitab 15. x 1 36 2 44 3 45 4 106 5 38 6 46 7 47 8 112 9 42 10 49 11 48 12 118 13 42 14 50 15

InDATA <-read.table("C:/Data/May 2004/season.txt",header=T) X <- decompose(InDATA) print(X) Period Connections Q1 67519 Q2 69713 Q3 68920 Q4 69452 Q1 70015 Q2 59273 Q3 57063 Q4 65596 Q1 73527 Q2 58586 Q3 69522 Q4 60091 Q1 51686 Q2 63490 Q3 55702 Q4 53200 Q1 51033 Q2 48175 Q3 52709 Q4 50106 Q1 50855 Q2 43466 Q3 48190 Q4 41702 Q1 48747 Q2 51441 Q3 42537

Dear R users, I have a very simple problem but I can't find the function in R to deal with it. I need to split (or decompose) one line into many lines using one field as a reference. I have a table with the following format: A B Frequency 23 3 2 24 2 5 25 1 3 And need to split each line into several lines according to the frequency to achieve something like this: A B

Hello, I'd like to know how R does calculate each component in the decompose() function? More precisely, how is calculated the final trend component in this function? Thanks for your answer Myriam -- View this message in context: http://www.nabble.com/Decompose-function-%3A-calculation-of-each-component-tp24362207p24362207.html Sent from the R help mailing list archive at Nabble.com.

Hi All, I am new to this forum. I have just started learning R. When i use plot(decompose(x)), i am getting the title " Additive time series decomposition". How to make this title off and change to some other title. Any help regarding this is highly appreciated. With sincerer regards, Upananda -- View this message in context:

Hi everyone! I would appreciate help with using decompose.graph(), community detection functions from igraph and lapply(). I have an igraph object G with vertex attribute "label" and edge attribute "weight". I want to calculate community memberships using different functions from igraph, for simplicity let it be "walktrap.community". This graph is not connected,

I just got this error: > library(igraph) > comp <- decompose.graph(gr) Error: protect(): protection stack overflow Error: protect(): protection stack overflow > what can I do? the digraph is, indeed, large (300,000 vertexes), but there are very many very small components (which I would rather not discard). PS. the doc for decompose.graph does not say which mode is the default. --

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

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 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,