similar to: statistical help

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

Does anyone have experience importing image files into R? I would like to be able to import such files, use R to analyze and compare and perhaps export again a modified file which can be displayed by the usual image display programs. -------------- next part -------------- An HTML attachment was scrubbed... URL: https://stat.ethz.ch/pipermail/r-help/attachments/20010303/ad72932b/attachment.html

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

I was using the hist() function to create a frequency table of some network traffic data. The range in values is rather large, from 0 till just under 10e12. Calling hist(x, breaks=c(0,1000,1e6,1e9,1e12),plot=F,freq=T) causes hist() to return : $breaks [1] -1.0000e+05 1.0100e+05 1.1000e+06 1.0001e+09 1.0000e+12 Is this recalculation of the breaks by hist() intended? Maarten van Gelder.

According to the CRAN contributed packages page, the license for "KernSmooth" is *not* GPL. It is listed as "License: Unlimited distribution (from Dr Wand)." Are there any restrictions on this of any kind? Is there some written license document similar to the GPL license? The reason I'm asking is that the legal beagles in the company I work for are extremely paranoid about

Hi, I've seen that handle missing values makes difficult to manage data series within R. I've look at Shumway's program (ASTSA) and Brockwell & Davis pogram (ITSM), it seems that both handle NA's values without problem, or maybe, I'm getting an spurious results? My problem is with two variables: rain and river flow: the first has 7 NA values (is to say, seven months, non

Dear R users, I'm looking for a test of bimodality in order to make some decisions about how to procede with an analysis algorithm. I have not come across any such tests in my readings and discussions apart from the Rao which appears to be applicable to cyclic data. The data I'm interested in characterizing as uni- or bimodal are frequency x amplitude spectra of consonant speech sounds,

Hi Everyone, Does anyone know how to take the output of a 2D fft() and shift the DC component to the center? Is there a function that does that? if you know, please send me an email. thanks for any help! brad __________________________________________________ . http://im.yahoo.com -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

When I try to go to http://cran.us.r-project.org I get sent to a location called "domains 4 sale" or something of that ilk, indicating that some "URL scavenger" has flagged it as expired. Has something happened, or is it just wishful thinking on some scavenger's part that the URL is no longer pointing to the US CRAN mirror? Try this link:

I tried to run the same dataset with R-1.2.3 on Linux, compiled from the source tarball. After 7 hours, it was still grinding away and had not crashed. I had to kill it. "gtop" was showing it with a real size around 160 MB, same as the peak on the Windows 2000 environment. This is on the same hardware as the Windows run; my machine is dual-booted Windows 2000 and Red Hat Linux 7.1. If I

Dear all, I have a lot plot graphs and i would like to save them. One solution is to use the grab function with xv. But it is tedious.... So is there a R function that can save automatically the result of a plot in a gif or ps or bmp file Thanks Olivier -- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- Olivier MARTIN phone: (33) 04 76 61 53 55 Projet IS2

I have a Red Hat 7.0 system, plus glibc-2.2 and gcc-2.96.69. When I run the R demo 'demo("graphics") with KDE it runs fine, but when I try it with Gnome it croaks: demo(graphics) ---- ~~~~~~~~ Type <Return> to start : > opar <- par(ask = interactive() && (.Device %in% c("X11", "GTK", "windows", "Macintosh")))

Hello, I have a problem. I have to read a HUGE file which has to be line processed, so I would like to use scan like p <- scan(pfile,what=c(0,0,0,"",0,""),nlines=1) but it continues to read the first line, so I have to add a skip and increment the skip after each read. It takes forever to read a test file of 1 Mb, the real one is in 60's Mb Is there a way to read line

Hi, How the eigenvectors output by the eigen() function are ordered. The first column corresponds to the largest eigenvalue? or is the last column as in Octave? I'm performing a spatial-temporal analysis of some climatic variables so my matrices are MxN (locations*time)and I'm looking for the leading EOF's. As I have understand the eigenvectors columns represent those EOF's

I'm trying to do a modulus operator in R-1.2.2 on Windows -- specific version: R : Copyright 2001, The R Development Core Team Version 1.2.2 Patched (2001-03-22) Here's what the R language manual says: "R contains a number of operators. They are listed in the table below. [snip] % Modulus, binary" Here's what happens on the R console: > records <- 100000 >

Is it expected behavior for gc() to return shrinking values as it gets called multiple times? Here's what I've got: > gc() used (Mb) gc trigger (Mb) Ncells 221754 6.0 467875 12.5 Vcells 3760209 28.7 14880310 113.6 > gc() used (Mb) gc trigger (Mb) Ncells 221760 6.0 467875 12.5 Vcells 3016206 23.1 11904247 90.9 > gc() used (Mb) gc

Is there any way to control the format of the axis label numbers on a plot? More specifically, I have some plots that get axes with label numbers in exponential format, and I'd like to change that to non-exponential. Thanks!! -- M. Edward (Ed) Borasky, MS, MNLP, NST, FBG, PGS & PTA znmeb at borasky-research.net http://www.borasky-research.net