similar to: calculate phase/amplitude of fourier transform function in R

Dear R People: Suppose I want to estimate the parameters of the following AR model: (1 - phi_1 B - phi_2 B^2 - phi_9 B^9) x_t = a_t and I want to use the arima0 command from the ts library. How would I use the order subcommand, please? R Version 1.4.1 for Windows. Thanks! Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston -

I've looked around but I can't find the method in R for testing whether the resulting estimated coefficients of an AR model imply that the model is invertible. To quote from eric zivot's blue book : " the AR(p) is invertible provided the rots of the characteristic equation Phi(z) = 1 - phi_1*z - phi_2*z^2 = phi_3*z^3 - ..... Phi_p*z^p = 0 lie outside the complex circle".

Hello! Is it possible to use R time series to identificate a process which is subjected to known input? I.e. I have 2 sequences - one is measurements of black box's state and the second is the "force" by which this black box is driven (which is known too) and I want to fit thist two series with AR-process. The "ar" procedure from stats package expects that the force is

Hi, Is this method broken in R? I am using it to find roots of the following function: f(x) = 2.5*exp(-0.5*(2*0.045 - x)) + 2.5*exp(-0.045) + 2.5*exp(-1.5*x) - 100 It is giving an answer of -38.4762403 which is not even close (f(x) = 2.903809e+25 for x=-38.4762403). The answer should be around 0.01-0.1. This function should converge.. Even for a simple function like f(x) = exp(-x) * x, it gives

Hi, I tried adding the PHI nodes in BB_unique, and it works for the simple case described here, but in case the nodes were declared in some predecessors of ExitBB1 and used in ExitBB1_redirect and its successors, it won't work, unless I create entries for all of them in BB_unique. B1 (declares PHI_1) B3 | | B2

On Oct 11, 2010, at 9:05 AM, Xinfinity wrote: > > Hi, > > I am working on a pass aimed to unify multiple exits of a loop into a unique > basic block. The approach is straight forward: > I create a unique BasicBlock BB_unique that has as predecessors all the exit > blocks of the loop, it contains a phi instruction and a switch to redirect > the flow correctly.

Hi, I am working on a pass aimed to unify multiple exits of a loop into a unique basic block. The approach is straight forward: I create a unique BasicBlock BB_unique that has as predecessors all the exit blocks of the loop, it contains a phi instruction and a switch to redirect the flow correctly. Additionally, for each initial exit block I create an associated block that will jump to the

Tian It appears the attachment might not have worked so I'll embed Bill's message at the end. Peter Alspach > -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of Peter Alspach > Sent: Thursday, 17 August 2006 8:02 a.m. > To: T Mu; R-Help > Subject: Re: [R] confusing about contrasts concept

Good morning, (sorry for my english) i have some problems to put off by extimation ARIMA coefficients the ones not significatives. Exist a method to extimate only that significatives? i use the command: arima(). thanks to all Stefano [[alternative HTML version deleted]]

Dear R-helpers, I have time series data from 16 subjects and 2 treatment groups. The seasonal variation can be best described by two harmonics, I called the frequencies omega and omega2. I now want to test whether (1) the seasonal pattern differs between the treatments (some kind of overall test). If this is the case, (2) I want to conduct tests to find out whether the amplitude of the

Hi dear R community, I have up to 12 measures of a protein for each of 6 patients, taken every two or three days. The pattern of the protein looks periodic, but the height of the peaks is highly variable. It's something like this: patient <- data.frame( day = c(1, 3, 5, 8, 10, 12, 15, 17, 19, 22, 24, 26), protein = c(5, 3, 10, 7, 2, 8, 25, 12, 7, 20, 10, 5) ) plot(patient$day,

Dear list, I apologies first for my English, hope you will understand well my question. I am working on 1/2 hour piezometric data, time unit is second. They present daily oscillation when using the spectrum() function. What I am really interested in, is to find the amplitude corresponding to this oscillation. I work with a college using Matlab, and although we apply the same methodology, our

Hi dear R community, I have up to 12 measures of a protein for each of 6 patients, taken every two or three days. The pattern of the protein looks periodic, but the height of the peaks is highly variable. I'm testing for periodicity using a Monte Carlo simulation envelope approach applied to a cumulative periodogram. Now I want to predict the location of the peaks in time. Of course, the

Hi all, I'd like to fit a multivariate regression with the variance of the error term porportional to the predictors, like the WLS in the univariate case. y_1~x_1+x_2 y_2~x_1+x_2 var(y_1)=x_1*sigma_1^2 var(y_2)=x_2*sigma_2^2 cov(y_1,y_2)=sqrt(x_1*x_2)*sigma_12^2 How can I specify this in R? Is there a corresponding function to the univariate specification lm(y~x,weights=x)??

Dear All, Here is what I am trying to achieve: I would like to plot some data in 3D. Usually, one has a matrix of the kind y_1(x_1) , y_1(x_2).....y_1(x_i) y_2(x_1) , y_2(x_2).....y_2(x_i) ........................................... y_n(x_1) , y_n(x_2)......y_n(x_i) where e.g. y_2(x_1) is the value of y at time 2 at point x_1 (see that the grid in x is the same for the y values at all times).

Dear R users, I have a response variable in a csv file called "y" and a matrix of predictor variables in a csv file called "mat". I have used the function "nls" I have specified the nonlinear relation between these variable.The code I have witten is called Rprog which begins with the phrase: L.minor.m1<-nls(Y~a ....etc.. The program when I execute the program, I

Hi all, I'm not familiar with R programming and I'm trying to reproduce a result from a paper. Basically, I have a dataset which I would like to model in terms of successive increments, i.e. (y denote empirical values of y) y_1 = y1, y_2 = y1 + delta1, y_3 = y1 + delta1 + delta2. ... y_m = y1 + sum_2^m delta j where delta_j donote successive increments in the y-values, i.e. delta

Dear listers In R, if I have a sequence x(t), t=1,...N, fft(x) is actually giving us sum(x(t)exp(-i*omega*t)) at fourier frequency omega= 2*pi*i/N, i=0,1,...(N-1). The question is if I want to calculate sum(x(t)exp(-i*2*omega*t)), how can I do it? thanks a lot! ------------------------------- liyun (Lauren) Ma Dept of Statistics North Carolina State University

Hello, I have a series and I need to run a Fourier Transform for that series. I have done that using the function fft from stats package. However I am not sure whether the result I am getting is correct or not. Seems that the first value of the Fourier Transform list is the sum of all elements of the given series and all the following numbers represent a palindromic series for some how. I have

This posting has nothing to do with R (except maybe that I am using R very heavily in writing the paper to which the question pertains.) I simply wish to draw upon the impressive knowledge and wisdom of the R community. Since this question is way off topic, if anybody has the urge to reply, they should probably email me directly: rolf at math.unb.ca rather than via this list. My question