similar to: How to reference R in a scientific document ?

Hi All, I have a time series called 'upwelling', like this: >upwell Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov 1984 494.7 303.7 220.8 288.4 188.1 125.5 215.7 1985 56.0 127.4 165.8 189.4 261.6 223.7 186.3 150.2 107.8 120.3 91.1 from 1984 until 2000. How do I could extract, i.e, Jan or Dec columns from the data? Or how do I

Hi, How do I can make a lagged scatterplot of two variables: Yt (nao) versus Xt-h (mei) if they have the following structure: >series mei nao Jan 1950 -1.036 0.55 Feb 1950 -1.133 3.31 Mar 1950 -1.259 0.81 Apr 1950 -1.027 1.60 May 1950 -1.399 -1.73 Jun 1950 -1.366 1.26 Jul 1950 -1.300 -0.87 . . . I've tried with lag.plot but I don't understanf how to use it Thanks in

Hi All, How do I can create a new vector, i.e. 'seasonal' or 'trend' from the resultant seasonal or trend component of the Time.Series object produced by STL, and how I could superpose in the same graphic, i.e. original data and trend or seasonality? Thanks in advance! Antonio Antonio Rodr?guez Verdugo CICEM Agua del Pino Huelva Oceanography and Coastal Resources, PhD Program,

Hi, I have a data frame with 13 time series (same length). Two of them have several NA's values, so when I've tried to run: >acf(flow, na.action=na.omit) I get >Error in na.omit.ts(as.ts(x)) : time series contains internal NAs I've read through help files of acf, na.action, na.omit, and I don't understand how to manage NA's values. If I list the variable (eg: rain,

Hi, Is there anybody who could explain a little bit the use of s.window and t.window in STL. If I decompose a series using the default values, I get a very irregular trends, but if I play with the options, I get less rare trends. If anybody can explain me the concepts behind these options I will appreciate a lot. I'm not a mathematician ;-) Thanks in advance Cheers, Antonio Antonio

Hi, I have recently installed R-1.2.2 for windows (16MB RAM, P-166) and I getting the following message after processing my data (6 variables and 1200 observations): >Error: cannot allocate vector of size 4 Kb >In addition: Warning message: >Reached total allocation of 15Mb: see help(memory.size) Then, the program close. With the last version, 1.1.1 (I think) I didn't have this

Hi, Are there some texts about time series that teatches it with R? And site 'bout time series "for dummies" , can anyone indicate them for me? thanks a lot Rafael Bertola Undergraduate in Statistics at UNICAMP - Brasil -------------------------------------------------- --------------------------------------------------

Dear all, I'm writing a lgrind definition for R. Lgrind is a program to produce pretty program listings for inclusion in LaTeX documents. Does anybody here have any experience with lgrindefs? My definition now looks like this: R:\ :pb=^\d?\p\d<-\dfunction\(\a\):\ :bb={:be=}:cb=\d#:ce=$sb=":se=\e":lb=':le=\e':id=.:\ :zb=@:ze=@:tb=%%:te=%%:mb=%\$:me=\$%:vb=%\|:ve=\|%:\

Dear all, I'm writing a lgrind definition for R. Lgrind is a program to produce pretty program listings for inclusion in LaTeX documents. Does anybody here have any experience with lgrindefs? My definition now looks like this: R:\ :pb=^\d?\p\d<-\dfunction\(\a\):\ :bb={:be=}:cb=\d#:ce=$sb=":se=\e":lb=':le=\e':id=.:\ :zb=@:ze=@:tb=%%:te=%%:mb=%\$:me=\$%:vb=%\|:ve=\|%:\

Hi I'm looking to a function in R equivalent to S-PLUS' "interp". Is there any implementation in R? If not, is there any similar interpolation function such that at the data points the interpolated values concides with the data? Thanks in advance P.J. Paulo Justiniano Ribeiro Jr Dept Maths & Stats - Fylde College Lancaster University Lancaster LA1 4YF - U.K.

Dear R members, Is there a way I can include more than one variable in a histogram? Many thanks in advance, Liqing, Liqing Zhang Department of Ecology & Evolutionary Biology 321 Steinhaus Hall U.C. Irvine Irvine, CA 92697-2525 Phone:(949)-824-7703 Fax: (949)-824-2181 Email: lzhang at uci.edu ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

How would you proceed with the following problem: I have an data frame containing a table relating raw scores on an iq test to 'decile-like' scores. data<-data.frame(age =c( 8, 9, 10, 11, 12, 13, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70), pers90=c(33, 39, 43, 47, 50, 54, 54, 52, 50, 49, 47, 46, 43, 42, 39, 37), perc75=c(27, 31, 37,

Dear all, There has been some time since I asked about interpolation in higher (>2) dimensions, and I must admit I failed to write a function to do this myself the last time, but eventually ended up doing it in MATLAB. I tried to translate the MATLAB code, but MATLAB code is so much more opaque than R (S) code, so I failed that too, mainly because I could only get one MATLAB session, I would

Dear all, I am interested by the various messages sent to the r-help mailing list these last few days about documentation on R. Yet, there are a lot of issues actually open on this topic, and it is not clear to me what is the "urgent need for documentation" on R. R comes with the 'Introduction to R' by Venables et al. I recommend the colleagues and students in my lab who are

Dear all, I have been trying to understand what runif() is telling me. I am generating lots of numbers (billions and billions (wow, I''ve dreamed about saying that for many years... :-) ), for a distribution that has the following quantile function: 1 / (2 * sqrt(1 - p)) (that is, the distribution has a lower cutoff) As you can imagine, this has rather heavy upper tail. I was

Sorry, I pasted the wrong file in earlier... this is the correct one: pValCalculator(b, n=20, m=20) { ind <- 1:min(c(n,m)) prob <- (1-pnorm(b,sd=std*sqrt(ind))) prob1 <- sum((n-ind+1)*(m-ind+1)*prob) prob1 } inputData <- scan("/users/lvssso/projects/LAMA/output/pValLamaScore.tmp", list(block1 = "",block2 = "",width1 = 0,width2 = 0,alignment = 0,score

I am trying to install R on a unix workstation (Sun-SPARC) but I am not sure which of the binaries to use. The only binaries available are for linux, nextstep, osf, and windows. Which one of them will work for a SUN sparc station?? Thanks. PS. I tried to compile the source code but it looks like I have a probelm with fortran. So, if I could install one of the binaries it would be a lot

Hi Rers, W/o trying to make a long story short (I hope you read that correctly), I have been working on a nice little R function ... Somewhere along the line I got an error that refers to a max function in my function (I think this is a vector vs. scalar issue): Error in max(..., na.rm = na.rm) : invalid "mode" of argument If that's not bad enough, even when I q() R and say

Hi all, how can I remove objects from the workspace which starts with a certain pattern, e.g. lm (lm1, lm2, lm3 etc). Wildcards won?t work. Thanks, Sven -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the

Dear all, I think I must be going crazy.... If I do > mean(rpois(1000000, 14)) I get: [1] 13.50420 and again: > mean(rpois(1000000, 14)) [1] 13.49896 > mean(rpois(1000000, 14)) [1] 13.50161 > mean(rpois(1000000, 15)) [1] 14.49250 > mean(rpois(1000000, 15)) [1] 14.49897 > mean(rpois(1000000, 14.5)) [1] 13.99689 > mean(rpois(1000000, 14.5)) [1] 13.9963 This was on a platform