[R] fft and significance

2012-10-25 Thread nuncio m
HI users, Is it possible to check the significance of peaks in a power spectrum in R Thanks and regards nuncio -- Nuncio.M Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 ph off 91 832 2525636 ph: cell 91 9890357423 [[altern

[R] Significance of spectral peaks

2013-12-23 Thread nuncio m
Dear useRs, I have a time series of length approcimately 55. Is it possible to find the significance of fft spectral peaks with R? thank you -- Nuncio.M Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 ph off 91 832 2525636 ph: cell 91 98903574

[R] timeseries

2010-07-05 Thread nuncio m
Dear useRs, I am trying to construct a time series using as.ts function, surprisingly when I plot the data the x axis do not show the time in years, however if I use ts(data), time in years are shown in the x axis. Why such difference in the results of both the commands Thanks nuncio -- Nuncio.

[R] acf

2010-07-06 Thread nuncio m
Hi list, I have the following code to compute the acf of a time series acfresid <- acf(residfit), where residfit is the series when I type acfresid at the prompt the follwoing is displayed Autocorrelations of series ‘residfit’, by lag 0. 0.0833 0.1667 0.2500 0. 0.4167 0.5000 0.5

[R] Wavelet

2010-07-07 Thread nuncio m
Hi useRs, Is it possible to get MORLET wavelet in R Thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version deleted]] __ R-hel

[R] acf significance levels

2010-07-14 Thread nuncio m
Dear useRs, How to save the correlations corresponding to the significance levels from ACF function Thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version dele

[R] tsdiag

2010-07-22 Thread nuncio m
HI list, I want to know whether tsdiag uses k-(p+q) as the lag in ljung box test. How is it possible to save those values nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version d

[R] auto.arima

2010-07-23 Thread nuncio m
HI list, I am using auto.arima from forecast package, I wonder whether its possible to save model orders to a seperate file Thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML ver

[R] sink function

2010-07-23 Thread nuncio m
I have the following code to write the output from auto.arima function. The issue is not in finding the model but to divert its out put fit to a file order_fit.txt. code runs but nothing is written to order_fit.txt where am I going wrong library(forecast) for (i in 1:2) { filen = paste("file",i,"

[R] arima

2010-10-13 Thread nuncio m
HI useRs, Is it required to remove mean before using ARIMA models thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version deleted]] __

[R] removing spatial auto correlation

2010-08-02 Thread nuncio m
Hi list, I am trying to fit arima model for a grid of 360x161x338 points, where 360x161 is the spatial dimension and 338 is the number of time steps I have, which is seasonal. For this purpose I used the auto.arima function in forecast package. After fitting residuals at each grid in spa

[R] band pass filter

2010-08-14 Thread nuncio m
Hello list, Is there any way to bandpass filter in R thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version deleted]] __ R

[R] removing columns

2011-05-03 Thread nuncio m
Hi list, I have a matrix with all elements of some columns are zeroes. Is it possible to remove these columns and create a new matrix nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML ve

Re: [R] documentation in R

2011-06-12 Thread nuncio m
Dear Arun, Prof. Rob Hyndman's forecast package has automated arima. You have to install 'forecast' library for that. in linux go to r-prompt and type install.packages(). Hope this helps nuncio On Mon, Jun 13, 2011 at 12:10 PM, siddharth arun wrote: > How we can call auto.arima in

[R] writing autocorrelation and partial auto correlation functions to a file

2010-05-20 Thread nuncio m
Dear All, I am very new to T. I need to fit a ARIMA model to my time series. So I found the auto correlation functions and partial auto correlation function in R. Now I want to save these valuse along with the significance levels to a file. How to do that?. I tried some function in R

[R] missing values in autocorelation

2010-05-31 Thread nuncio m
Hi all, I am trying to find the autocorrelation of some time series. I have say 100 files, some files have only missing values(-99.99, say). I dont want to exclude these files as they represent some points in a grid. But when the acf command is issued i get an error. Error in plot.window

[R] ARIMA order

2010-06-03 Thread nuncio m
Hi all, Is there any way in R to select the order of an ARIMA model automatically nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML version deleted]]

[R] prewhiten

2010-06-07 Thread nuncio m
HI all., I have some univariate time series that need to be prewhitened. HOw this can be performed in R. I am thinking of to fit an ARIMA model and substract this from the original series. Is this the correct way THanks in advance nuncio -- Nuncio.M Research Scientist National Center for Antarcti

[R] Extracting columns from a class

2011-03-17 Thread nuncio m
Hi list, I am not a frequent user of R. Recently I used R in principal component analysis and got the result as a class, which has information like standard deviation and principal components from 1 to 10. How is it possible to extract the column corresponding to first principal compon

[R] principal components

2011-03-31 Thread nuncio m
HI all, I am trying to compute the EOF of a matrix using prcomp but unable to get the expansion co-efficients. is it possible using prcomp or are there any other methods thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Go

[R] svd

2011-04-04 Thread nuncio m
Dear list, I searched the libraries but could not find means to compute the svd of a coupled field. Is it possible in R Thanks nuncio -- Nuncio.M Research Scientist National Center for Antarctic and Ocean research Head land Sada Vasco da Gamma Goa-403804 [[alternative HTML ve

Re: [R] svd

2011-04-04 Thread nuncio m
)%*%B. but this differs significantly from cov(A,B). Thanks nuncio 2011/4/4 Juan Carlos Borrás > m <- matrix(c(1:12), nrow=3, ncol=4) > svd(m) > > > > On Mon, Apr 4, 2011 at 11:51 AM, nuncio m wrote: > > Dear list, > >I searched the libraries but coul