similar to: I'm trying to use package ts (decompose). How do you set up the data/ See attached. thanks

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

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

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

Hi everyone. I'm having some troubles with STL function to decompose some data. My issue is that I have monthly data from September 2005 up to August 2010 i. e. 60 observations. I define it in the following way: *u<-read.csv("C:/CELEBREX.csv",header = TRUE) u.ts<-ts(u, start=c(2005,9), frequency=12) * The issue is that when I try to use stl(u.ts, 'per') Then the

Is there any way to give a "decent" title after I plot something generated by decompose? For example: # generate something with period 12 x <- rnorm(600) + sin(2 * pi * (1:600) / 12) # transform to a monthy time series y <- ts(x, frequency=12, start=c(1950,1)) # decompose z <- decompose(y) # plot plot(z) Now, the title is the ugly "Decomposition of additive time

Is there a simple way to decompose the upper triangle of a correlation matrix to a linear list; For example: X Y Z X 1 2 3 Y 2 1 4 Z 3 4 1 so you get a list like: xy 2 XZ 3 YZ 4 I suspect you can do it with a matrix transformation, but that beyond me at present. Many thanks Mark _________________________ Department of Molecular and Human Genetics, Baylor College of Medicine,

Consider the following analysis of a class experiment done as a Latin Square: > spinner <- gl(4,4,16,label=c("Murray","Angela","Shasha","Stephen")) > order <- gl(4,1,16) > treat <- scan() 1: 1 2 4 3 5: 4 3 1 2 9: 3 4 2 1 13: 2 1 3 4 17: Read 16 items > coin <-

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

Dear all, I have a vector of length 1 but within this vector there are 9 elements > vectorNOcorrection [1] "(O, O, O, H, L, O, O, O, O)" can someone help me with how I can have it as with 9 element separated such that it will have length 9 as "O", "O"," O", "H","L", "O", "O", "O", "O" Kind

I did composition of butterfly8 and butterfly16 and I found, that these functions are well decomposed - decomposition doesn't lower computional speed. On the other hand the same can be done with butterfly8 - decomposition to butterfly4 (further decomposition is not possible) but there's no reason to do this. I think little improvement can be done by inlining them. Compiler and processor

Hi Duncan, On Wednesday 02 May 2012 10:15:17 Duncan Sands wrote: > yup, front-ends have to take care of more complicated ABI details. For > example the front-end should currently use "byval" for any (parts of) > structs that need to be passed on the stack, and explicit scalars for > struct bits that should go in registers. > > If this interface could be improved, I

Hi. I'd like to find some way for Samba to convert the filenames created by Macintosh OS X into something windows can use. According to their docs: http://developer.apple.com/techpubs/macosx/Essentials/SystemOverview/FileSystem/chapter_9_section_13.html filenames are stored as canonically decomposed Unicode 2.1 in UTF-16 format (a sequence of 16-bit codes) Samba running on a

Greetings! For some reason I am not managing to work out how to do this (in principle) simple task! As a result of applying strsplit() to a vector of character strings, I have a long list L (N elements), where each element is a vector of two character strings, like: L[1] = c("A1","B1") L[2] = c("A2","B2") L[3] = c("A3","B3")