Displaying 20 results from an estimated 400 matches similar to: "arima on defined lags"
2009 Mar 04
0
R under Citrix and access to Lotus notes
Dear All,
1. Does anyone have experience of running R on a server inside a Citrix
shell - I'd like to get R onto the server and would be greatful for any
tips or direction on the matter.
2. This may seem like a silly question so forgive my ignornace.
Most of the data I currently work with is held on a number of Lotus Notes
(LN) Databases (well it's called a DB here but it's really a
2008 Nov 28
1
confidence interval for glm
Hi all,
simple Q:
how do I extract the upper and lower CI for predicted probabilities
directly for a glm - I'm sure there's a one line to do it but I can't find
it.
the predicted values I get with the predict (.. "response")
Thanks
Gerard
**********************************************************************************
The information transmitted is intended only for
2008 Nov 20
1
binomial glm???
Hi everyone,
newbee query!
I've installed R 2.8.0 and tried to run this simple glm -
x is no of cars in a given year, y is the number voted in an election
that year while n is the population 18+:
votes <- data.frame(x = c(0.62,0.77,0.71,0.74,0.77,0.86,1.13,1.44),
+ y=c(502,542,711,653,771,806,934,1123), n=
2008 Dec 22
2
queue simulation
Hi all,
I have a multiple queing situation I'd like to simulate to get some idea of
the distributions - waiting times and allocations etc.
Does R has a package available for this - many years ago there used to be a
language called "simscript" for discrete event simulation and I was
wondering if R has an equivalent (or hopefully with graphics, something
better!).
Apologies if there
2009 Jan 13
1
deviance in polr method
Dear all,
I've replicated the cheese tasting example on p175 of GLM's by McCullagh
and Nelder. This is a 4 treatment (rows) by 9 ordinal response (cols)
table.
Here's my simple code:
#### cheese
library(MASS)
options(contrasts = c("contr.treatment", "contr.poly"))
y = c(0,0, 1, 7, 8,8,19, 8,1, 6,9,12,11, 7,6, 1, 0,0, 1,1, 6, 8,23,7,
2009 Jan 14
1
loglm fitting
Dear all,
sorry to bother you all with this but I've been trying to use the loglm in
MASS package (v2.8.0) and cannot get any sensible output.
I'm wondering am I doing something very foolish or missing something
obvious.
For example, I tried the documentation help(loglm) example - here's the
code
# Case 1: frequencies specified as an array.
sapply(minn38,
2008 Dec 09
2
for loop query
Hi all,
apologies if this is obvious - but I can't see it and would appreciate some
quick help!
the matrix mhouse is 26x3 and I'm computing odds ratios. The simple code
below "should" compute the odds vector for every pair (325) i.e. 26C2 in
cols 1 and 2.
On the first i=1 outer loop the inner j loop runs from 2 to 26 ok and then
I get the error (Error: subscript out of bounds)
2003 Apr 30
2
Bug in arima?
I'm using the fixed argument in arima. Shouldn't ar4, ar5, and ar6
display as zero in the output?
Call:
arima(x = window(log(hhprice), start = c(1990, 1), end = c(2003, 3)),
order = c(7,
1, 0), xreg = window(ts.union(exa1 = lag(exa, -1), exa12 = lag(exa,
-12), exb1 = lag(exb, -1), exc1 = lag(exc, -1), exc12 = lag(exc,
-12)), start = c(1990, 1), end = c(2003, 3)),
2011 Sep 09
2
Different results with arima in R 2.12.2 and R 2.11.1
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
2011 Sep 12
1
Difference in function arima estimation between 2.11.1 and R 2.12.2
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
2003 Nov 24
0
link between arima and arma fit
Hi dear sirs,
I am wondering why the fit of the time serie x with an arima and the fit of
diff(x) with an arma (same coeff p & d) differ one from another
here are the output of R:
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
> modelarma<-arma(diff(x),c(7,5))
> modelarma
Call:
arma(x = diff(x), order = c(7, 5))
Coefficient(s):
ar1 ar2 ar3 ar4 ar5 ar6 ar7 ma1 ma2
0.06078
2011 Feb 16
0
Arima contents
Hello,
I'm running a number of arima models using the "arima" function. Often,
when lag length gets too high, these model don't converge and an error
message appears as this:
> reg <- arima(y,order=c(7,0,7),xreg=isr)
Warning message:
In arima(y, order = c(7, 0, 7), xreg = isr) :
possible convergence problem: optim gave code=1
In this case, when you print the results
2009 Feb 20
0
residuals from a fractional arima model and other questions
Dear list and Martin,
I'm testing different approaches to fit an electricity demand time series and come upon the fracdiff package (v 1.3-1) for fitting fractional ARIMA models. The following questions are motivated by this package.
1. Despite having a help page, the residuals and fitted functions don't seem to have implementation, or did i miss something obvious? Alternatively, having a
2006 Nov 22
2
problems with garchFit
Hi all,
I post it on both r-help and r-finance since I don't know where is most
appropriate for this topic. Sorry if it bothers you.
I did garch fitting on S&P500 monthly returns with garchFit from fSeries. I
got same coefficients from all cond.dist except normal. I thought that is
probabaly usual for the data. But when I play with it, I got another
question.
I plot skew normal with
2008 Nov 09
3
Arms Race
hey can anybody help me? i have to simulate the richardson Arms race
model on R.. for my simulation class...
2004 Jan 14
2
Fixed parameters in an AR (or arima) model
Hello
I want to fit an AR model were two of the coefficients are fixed to zero
(the second and third ar-coefficients).
I used the "arima" function with the "fixed" argument but the ar3
coefficient is not set to zero:
==============================================
> arima(Y, order=c(4,0,0), xreg=1:23, fixed=c(NA,0,0,NA,NA,NA))
Call:
arima(x = Y, order = c(4, 0, 0), xreg =
2009 Jan 27
2
optim() and ARIMA
dhabby wrote:
Last week I run in to a lot a problems triyng to fit an ARIMA model to a
time series. The problem is that the internal process of the arima
function
call function "optim" to estimate the model parameters, so far so good...
but my data presents a problem with the default method "BFGS" of the
optim
function, the output error looks like this:
2009 May 07
2
Linear least squares fit with errors in both x and y values.
HI,
I'd like to perform a weighted linear least squares fit with R on data
with varying errors on both vectors. I can do this with one axis using
lm, but have no idea where to go from here. I've tried googling, but no
idea. Any suggestions?
Thanks,
James
2009 Feb 03
1
SAS language to R :interview
Dear List,
Please find a frank interview with Phil Rack, creator of Bridge to R (
from both SAS and WPS interfaces).
For those unaware of WPS- it is basically a SAS language compiler
(read SAS code,writes SAS code,Reads and writes SAS datasets) ,priced
at 660 $ a licence ( or estimated 10 times cheaper than Base SAS.
The UK based WPC held, WPS doesnt have advanced statistical facilities
like
2009 Jun 03
2
how can I ordinal regression??
What function and package I use to conduct ordinal regression??
My data is composed 2colums and 180rows.
The first colum indicate level of mass and second colum is intensity.
So, I want to calculate how much intensity are related mass.
[[alternative HTML version deleted]]