[R] comparing ARIMA model to data

hi, i am trying to teach myself about ARIMA models. i have followed examples from a number of sources and have more or less got the hang of how it works. i would like to compare the output from the fitted model to the original data. is this possible? or even a meaningful thing to do? to be clear,

Re: [R] asterisk in subscript

magic! this does the trick: expression(paste(italic(n)[symbol("\052")])) thanks for the hint, ted! -- Andrew B. Collier Physicist Waves and Space Plasmas Group Hermanus Magnetic Observatory Honorary Senior Lecturer tel: +27 31 2601157 Space Physics Research Institu

Re: [R] asterisk in subscript

l: +27 31 2601157 Space Physics Research Institute fax: +27 31 2607795 University of KwaZulu-Natal, Durban, South Africagsm: +27 83 3813655 On Fri, 2011-01-28 at 09:48 +, Prof Brian Ripley wrote: > On Fri, 28 Jan 2011, Andrew Collier wrote: > > > hi,

[R] asterisk in subscript

hi, i am trying to label a plot axis with the equivalent of the latex $n_*$. i initially tried expression(paste(italic(n)["*"])) but this made the * absolutely tiny and centred about midway wrt the n. then expression(paste(italic(n)[textstyle("*")])) made the * about the right size but now it

[R] lines and points without margin

hi, i am sure that this is a trivial question but i have not been able to find an answer by searching the mailing lists. i want to plot points on a graph, joined by lines. the command that i am using is points(x, y, type = "b", pch = 21) this plots nice open circles at the data points and draws

Re: [R] scale caption on levelplot

hi peter and david, thanks for the excellent suggestions. here is something like what i am finally using (those fancy fonts were really tempting, but i chose something a little more mundane!): library(lattice) x <- sort(rnorm(100,50,10)) y <- sort(runif(100,0,20)) d <- expand.grid(x=x, y=y) d$z

[R] scale caption on levelplot

hi, i am trying to figure out how to put a caption on the colour scale of a levelplot. there does not seem to be an option for this in levelplot(). i tried using mtext() but as soon as you put the text far out enough on the right of the plot, it goes beyond the plot boundary. so i tried to extend

[R] save() with 64 bit and 32 bit R

hi, i have been using a 64 bit desktop machine to process a whole lot of data which i have then subsequently used save() to store. i am now wanting to use this data on my laptop machine, which is a 32 bit install. i suppose that i should not be surprised that the 64 bit data files do not open on m

[R] ts subscripting problem

hi, i am having trouble getting a particular time series to plot. this is what i have: > class(irradiance) [1] "ts" > irradiance[1:30] 197811 197812 197901 197902 197903 197904 197905 197906 1366.679 1366.729 1367.476 1367.739 1368.339 1367.883 1367.916 1367.055 197907 197908

[R] p values for polychor

hello, i have been using cor.test() for calculating the correlation coefficient and p values for some data. however, since the data consist of two dichotomous sequences (actually just binary data), i understand that simply using the pearson correlation is not sufficient. however, having done a

[R] cor.test() and binary sequences

peter, thanks for your help with my questions regarding cor.test(). i have another question though: does this function make any assumptions about the underlying distribution of the two sequences? does it assume that they have a gaussian distribution? i ask because the data that i am working wit

[R] (no subject)

hello, i am a bit of a statistical neophyte and currently trying to make some sense of confidence intervals for correlation coefficients. i am using the cor.test() function. the documentation is quite terse and i am having trouble tieing up the output from this function with stuff that i have r