similar to: Outlier Detection for timeseries

Hi R users, Is there an easy way in R to generate the results table below using table 1 and the formula (simplified version of the real problem)? It would be easy if I knew the R equivalent of SAS's retain function, but could not find one. Thanks in Advance for any help! table1: ID X2 X3 1.00 1.00 0 2.00 0.00 3.00 1.00 4.00 3058 5.00 0.00 6.00 6.00 Formula: X3 = x2 + (.24 *

Hi R users, I am doing cross correlation analysis on 2 time series (call them y-series and x-series) where I need the use the model developed on the x-series to prewhiten the yseries.. Can someone point me to a function/filter in R that would allow me to do that? Thanks in advance for any help! -- View this message in context:

Hi R users, I have a data frame that contains 10K obs and 200 variables where I am trying to format the numeric columns to look like the output table below (format to 2 decimal places) but I am having no luck.. Can someone tell me the best way to accomplist this? Thanks in advance for any help! str(ad.test) 'data.frame': 10,000 obs. of 200 variables: $ ID : Factor

Hello R users, I have been trying to output all my results (text, plots, etc) into the same postscript file as one document, but have been unable to...Can anyone help me improve my code below so that I can accomplish this? Currently I have to output them separately then piece them back together into one document.. Thanks in Advance for any help! options (scipen=999, digits=7)

Working - among other things- in the field of (short & long term) electricity forecast, we are now using too many & too expensive pieces of licensed software: SAS, SPSS, EViews. This "sedimentation" is due to the fact that my predecessors in the past used different consultant companies to manage each procedure. Having attended the useR2004! Conference with the aim of assessing

Dear R People: When looking at the previous postings regarding TRAMO-SEATS, I am somewhat puzzled. Is it true that we CANNOT replicate TRAMO-SEATS because of licensing or ownership issues, please? If not, would anyone be interested in an R version of it, please? Thanks, Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston -

I have a Mac OS X system. To deal with a long monthly electricity demand time-series I use the procedures TRAMO/SEATS with the MS-windows only Demetra programme and X12 under R resorting to the awkward - as far as the output is concerned - x12 R package running the relating Fortran code. I wonder if someone out there has attempted to translate TRAMO/SEATS and X12 into R native language? Ciao

Dear Colleagues would someone know a suitable online-source for information regarding the TRAMO SEATS method for time series ? (an alternative to X12 ARIMA and earlier ARIMAs used by the US census bureau) Cheers -- ----------------------------- Soren Wilkening Principal Consultant phone: +49-30-96080121 wilkening at censix.com CENSIX Consulting Statistics, Surveys, Censuses

Are any seasonal adjustment programs, like Tramo/Seats, Census X12 ARIMA or Berliner Verfahren implemented in R? I am doing a simulation study and I don't know how to adjust the series in R. The possibility to access external the exe.files of the seasonal adjustment programs seems to be quite difficult. Can anyone help me? Thanks, Ingo

Hi R users, I have a time series variable that is only available at a monthly level for 1 years that I need to decompose to a weekly time series level - can anyone recommend a R function that I can use to decompose this series? eg. if month1 = 1200 I would to decompose so that the sum of the weeks for month1 equals 1200, etc.. Many thanks in advance for any help. -- View this message in

Hi R users, I have an R object with the following attributes: > str(sales.bykey1) 'by' int [1:3, 1:2, 1:52] 268 79 118 359 87 147 453 130 81 483 ... - attr(*, "dimnames")=List of 3 ..$ GROUP: chr [1:3] "III" "II" "I" ..$ year : chr [1:2] "2006" "2007" ..$ week : chr [1:52] "1"

Hello, i am a newby on R and i am trying to make a backward selection on a binomial-logit glm on a large dataset (69000 lines for 145 predictors). After 3 days working, the stepAIC function did not terminate. I do not know if that is normal but i would like to try computing a "homemade" backward with a repeated glm ; at each step, the predictor with the max pvalue would be

It might help if you provided the code you used. It's possible that you didn't use direction="backward" in stepAIC(). Or if you did, it was still running, so whatever else you try will still be slow. The statement "R provides only the pvalues for each level" is wrong: look at the anova() function. Bob On 29 June 2017 at 11:13, Beno?t PELE <benoit.pele at

Hola, Mira la versión de ggplot2 que estás usando por tenerla actualizada a la última... A mi me funciona introduciendo algún cambio salvo la última línea cuando cambias el tema... #----------- Lines <- " id Rangos Clase Asistentes 1 45d EPI 405 2 1a EPI 812 3 2a EPI 1639 4 3a EPI 2202 5 4a EPI 1994 6 5a EPI

?? Hola. Estoy intentando ordenar el gráfico tal que las columnas se apilen en el siguiente orden: debajo, los EPI; los SemiEPI en el medio y los NoEpi arriba de todo. Estos son mis datos: > tabla.barras Rangos Clase Asistentes 1 De 45d EPI 405 2 De 1a EPI 812 3 De 2a EPI 1639 4 De 3a EPI 2202 5 De 4a EPI 1994 6 De 5a EPI

Junto con saludar, Quer?a saber si los que han usado boxplot.stats, saben como ubicar los resultados que arroja como "out" (outliers) directamente en la base de datos, ya que, muestra los resultados en valor de cada dato at?pico. Por ejemplo, estoy usando una base de datos con 300 datos y 10 columnas, en este caso estoy buscando los datos at?picos de la columna del precio de x

I am using pyramid.plot() from the plotrix package. I have something like this ############################################ xy.pop<-dados$masfr xx.pop<-dados$femfr #agelabels<-dados$femlab xycol<-color.gradient(c(0,0,0.5,1),c(0,0,0.5,1),c(1,1,0.5,1),11) xxcol<-color.gradient(c(1,1,0.5,1),c(0.5,0.5,0.5,1),c(0.5,0.5,0.5,1),11) xylab<-dados$maslab xxlab<-dados$femlab

Hi, again i am stuck in my presentation, and i have never learn R before in my life but need this to be done, so please help me out for a favour: http://www.nabble.com/file/p20333155/kew.dat kew.dat run this in R and these comes up: Month Year Rain 1 Jan 1900 74.400000 2 Feb 1900 80.500000 3 Mar 1900 23.600000 4 Apr 1900 23.600000 5 May 1900 25.100000 6

On Mon, 24 May 2004 12:00:46 +0200 v.demartino2@virgilio.it wrote: > Working - among other things- in the field of (short & long term) electricity > forecast, * * * > we have to comply with the Tramo-seats closed-source procedure (http://www.bde.es/informes/be/docs/dt0014e.pdf) > to deal with seasonality of electricity monthly time-series, in line with > the methodology

Hello R users, I am new to R and am wondering if anyone can help me out with the following issue: I wrote a function to build ts models using different inputs, but when R displays the call for a model, I cannot tell which variables it is using because it shows the arguments instead of the real variables passed to the function. (e.g Call: lm(formula = dyn(dep ~ lag(dep, -1) + indep)) --->