similar to: ARIMA - h-step ahead errors

Hi, When I run this code a part of my x-axis and y-axis changes color. Can somebody tell me what is wrong? Also, is there a way to control the color of the average lines? monthplot(AirPassengers+500, ylim=c(min(AirPassengers), max(AirPassengers+500)), ylab="") par(new=T) monthplot(AirPassengers, col="blue", ylim=c(min(AirPassengers), max(AirPassengers+500)),

Dear colleagues, Can anyone of you tell me how to write a "nabla" symbol in an R graph? Thanks in advance, Nuno ______________________________________________ Centro de Oceanografia - IO-FCUL, Portugal Center for Quantitative Fisheries Ecology - ODU, USA [[alternative HTML version deleted]]

Can someone explain me this spec.pgram effect? Code: period.6<-c(0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10 ,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10,0,0,0,0,0,10) period.5<-c(0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10 ,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0,0,0,0,0,0,10,0,0,0,0,10,0) par(mfrow=c(2,1))

Dear colleagues, I am new to R and statistics so please keep that in mind. I have doubts on the df calculation of Ljung-Box test (Box.test). The function seems to use always the df=lag=m and not df=m-p-q like suggested in Ljung and Box (1978) paper (that is referenced). Do you agree with this? If so, is there an R package function that computes Ljung-Box test with the degrees of

Dear colleagues, I am relatively new to R and time series and so I am experiencing difficulties in interpreting the output of "arima" in MA models (but not in AR models). I cannot make sense of the 1st innovations returned by "arima". In an AR(1) model I expect data[t]=phi1*data[t-1]+a[t] and in a MA(1) model data[t]=a[t]+theta1*a[t-1]. My interpretation from R-help is

Dear R colleagues, This seems pretty straight forward but I have been banging my head on this for some time and can't seem to find a solution suppose I have something like a1<-1; a2<-2; a3<-3; a4<-4; b1<-3; b2<-4 I would like to quickly assign to objects with a certain pattern, e.g., those in ls(pattern="a") a specific value, e.g., "99", without

I have been trying to read in a large data set using read.table, but I've only been able to grab the first 50,871 rows of the total 122,269 rows. > f <- read.table("http://www.cs.odu.edu/~fmccown/R/Tchange_rates_crawled.dat", header=TRUE, nrows=123000, comment.char="", sep="\t") > length(f$change_rate) [1] 50871 From searching the email archives,

Hi, I am far from experienced in both R and time series hence the question. The code for spec.pgram() seems to involve a circularity of the kernel (see below) yielding new power estimates to all frequencies computed by FFT. " if (!is.null(kernel)) { for (i in 1:ncol(x)) for (j in 1:ncol(x)) pgram[, i, j] <- kernapply(pgram[, i, j], kernel, circular = TRUE)

Hi, James S. Clark - Statistical computation for environmental sciences in R - mentions a function centralm (pg25) that I believe should be present in the base package but I can't find it. It is supposed to calculate means, variances, skewness, kurtosis of arrays. Does it exist in some other package? Thanks, Nuno Prista [[alternative HTML version deleted]]

I am using R 2.3.0, Is there a way to get prediction intervals using predict.gls? Any ideas would be appreciated. thanks, Spencer On 11/2/06, Frank McCown <fmccown@cs.odu.edu> wrote: > > I was wondering if anyone knows who should be contacted to add to the R > user-contributed documentation at > > http://www.r-project.org/other-docs.html > > There doesn't

Hello, I am analyzing some climate time series data using the Mann Kendall package and was wondering if there was a way to calculate the trend using Sen's nonparametric estimator slope in R? Thank you in advance, Barry _________________________ Barry Baker, Ph.D. Global Climate Change Initiative The Nature Conservancy 2424 Spruce St., Suite 100 Boulder, CO 80302 Tel: (303)-541-0322 Fax:

My stat textbook tells me that using Shapiro-Wilk test for each variable one by one is not equal to a test for multivariate normality as a whole. Does R have a function of testing for multivariate normality? Thanks. Hongsheng (Hank) Liao, Ph.D. Lab Manager Center for Quantitative Fisheries Ecology 800 West 46th Street Old Dominion University Norfolk, Virginia 23508 Phone:757.683.4571

I have a power law plot in log-log scale using plot(x,y,log="xy") but I can't get the regression line to plot correctly. abline(lm(log(y)~log(x))) plots a straight line that looks promising but the intercept with the y-axis seems to be way off. summary(lm(log(y)~log(x))) gives the values but i cant really make any sense out of it. fitted(..) just confirms the values that seem to be

Hi All I don't know whether i am asking on a correct list or not i have a query that how to know that what system call are called when packet enters the system or when we send a packet from a system? Thanks

Attached is a histogram chart with a break on y-axis which my friend made for me using phot shop. Do anyone know how to make such a break using R? Thanks in advance. (See attached file: 2007 age distribution.doc) Hongsheng (Hank) Liao, Ph.D. Lab Manager Center for Quantitative Fisheries Ecology 800 West 46th Street Old Dominion University Norfolk, Virginia 23508 Phone:757.683.4571

Where can I download Bhat package? I need to use plkhci function in this package, which estimates confidence interval using likelihood ratio. Thanks. Hongsheng (Hank) Liao, PhD. Lab Manager Center for Quantitative Fisheries Ecology Old Dominion University 757-683-4571 [[alternative HTML version deleted]]

I was wondering if there was any way to completely remove the borders around each bar when using the barplot function. I understand how to change the color of the border (border="white"), but I need the border to be invisible so I can plot many skinny bars right next to each other with only the bar's color visible. Thanks, Frank -- Frank McCown Old Dominion University

I setup samba on our Solaris 2.5 machine in ten minutes and everything has worked wondefully. Our linux box has been another story so I hope someone can help. We're now running redhat 5.1 and samba 1.9.18p10-51.1 on a pentium 200. The daemons appear to be running nicely. In going through the diagnosing document, everything is great until step 4. I do nmblookup -B courses__SAMBA__ (courses is

I am trying to fit a logistic model to my 10 year data (1999-2008) by year. Codes like below: Year <- c(1999: 2008) for(y in 1:length(year)) { file.input <- paste("C:\\", year[y], "\\data.csv", sep="") table <- read.csv(file=fileinput, header=TRUE, as.is=TRUE, na.strings=c("")) initial <- getInitial(percent ~ SSlogis(age, Asym,

This seems like a very basic question so I have searched the FAQ site, several manuals, and the R-Seek site as well as Googling for an answer but can't seem to come up with one. I am trying to import an Excel file that resides on a University network. The path name of the file is to a network drive and includes blank spaces (e.g. K:/science/next directory science2/mydata/). The setwd() works