similar to: noise in time series

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:

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

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

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

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

Hi list, I would like to normalize a matrix (two actually for comparison) using iterative proportional fitting. Using ipf() would be the easiest way to do this, however I can't get my head around the use of the function. More specifically, the margins settings... for a matrix: mat <- matrix(c(65,4,22,24,6,81,5,8,0,11,85,19,4,7,3,90),4,4) using fit <-

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

Is Pearson's Chi-Square test for contingency tables asymptotically unbiased for large tables (large degrees of freedom) regardless of the expected values in each cell? The rule of thumb is that Pearson's Chi-square should not be used when large numbers of cells have expected values < 5. However, I compared the results on 4x4 contingency tables for R's chisq.test using chi-square

Dear R experts and statisticians, I have some time series datasets, they are several years vegetation indices (about 50 data points per year) sampled from different station. These indices have similar dynamics with seasonal change. My questions are, 1) How can I compare the difference among the indices, and how can I say there is significant differnce between two time series. They

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,

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)

Dear list fellows, I want to study time series and use R to analyse time series of fishing data from several species (landings and cpue) investigating the correlation between them and with environmental factors (water temperature, wind, etc.). Searching at Amazon I found three books with examples in R: Time Series Analysis: With Applications in R by Jonathan D. Cryer and Jonathan D. Cryer

Hello, Anyone who knows a fast and accurate algorithm for generating draws from an underdispersed Poisson distribution. Or even better, if there is a package containing such an implementation. Thanks Rene

Hello, Guys: I'm from China, my English is poor and I'm new to R. The first message I sent to R help meets some problems, so I send again. Hope that I can get useful suggestions from you warm-hearted guys. Thanks. I builded a multiplicative seasonal ARMA model to a series named "cDownRange". And the order is (1,1)*(0,1)45 The regular AR=1; regular MA=1; seasonal AR=0; seasonal

I am a beginner in using R and I need help in the interpretation of AR result by R. I used 12 observations for my AR(2) model and it turned out the intercept showed 5.23 while first and second AR coefficients showed 0.40 and 0.46. It is because my raw data are in million so it seems the intercept is too small and it doesn't make sense. Did i make any mistake in my code? My code is as follows:

Dear all, I have the following dataset: each row corresponds to count of forest floor small mammal captured in a plot and vegetation characteristics measured at that plot > sotr plot cnt herbc herbht 1 1A1 0 37.08 53.54 2 1A3 1 36.27 26.67 3 1A5 0 32.50 30.62 4 1A7 0 56.54 45.63 5 1B2 0 41.66 38.13 6 1B4 0 32.08 37.79 7 1B6 0 33.71 30.62

If anyone wants to see a prime example of how inefficient it is to program in SAS, take a look at the SAS programs provided by the US Agency for Healthcare Research and Quality for risk adjusting and reporting for hospital outcomes at http://www.qualityindicators.ahrq.gov/software.htm . The PSSASP3.SAS program is a prime example. Look at how you do a vector product in the SAS macro

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,

Hi! I am calculating correlation between two variables: 1. X versus Y 2. X versus Y(with a 3 steps lag) I would like to test if the correlation increase/decrease from 1 to 2 is significant or not. Is there any function in R to do this? any hints? Thanks for help :) David Ria?o Center for Spatial Technologies and Remote Sensing (CSTARS) University of California 250-N, The Barn One Shields

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