[R] Intervention/Impact analysis in time series

2003-07-23 Thread Kenneth Cabrera
Hi R users: Does any one knows about a R library for deal with intervention/impact analysis in time series (eg. Box-Tiao et. al. theory?). Thank you for your help -- Kenneth Roy Cabrera Torres Celular +57 (315) 405 9339 __ [EMAIL PROTECTED] mailing list

RE: [R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread M. Edward Borasky
I haven't gotten around to assembling the toolset required to build R on Windows, since most of what I do is smallish interactive problems. However, another possibility would be to load CygWin/XFree86 on your laptop (which I've done), then download Atlas 3.5.7 from SourceForge (which I've done), th

[R] pls regression - optimal number of LVs

2003-07-23 Thread Dowkiw, Arnaud
Dear R-helpers, I have performed a PLS regression with the mvr function from the pls.pcr package an I have 2 questions : 1- do you know if mvr automatically centers the data ? It seems to me that it does so... 2- why in the situation below does the output say that the optimal number of latent

Re: [R] trouble with maps

2003-07-23 Thread Ray Brownrigg
[EMAIL PROTECTED] wrote: > On Wed, 23 Jul 2003, Deborah Swayne wrote: > > > Has anyone else seen this behavior from the "maps" package? > > > > map('state', fill=TRUE) > > > > results in a lively mix of overlapping polygons inside a > > map of the US, but they have no obvious relationship to

[R] negative binomial

2003-07-23 Thread Bruno Paul Mmbando
I am trying to model counts data from three different sites where the number of zeros differs significantly from one sight to another and also the counts gets bigger in site with few number of zeros considerably, meaning that the k's are different. In my model I am including about 10 covariates

[R] unz( "x.zip", "y.csv" ) != pipe( "unzip -p x.zip y.csv" )

2003-07-23 Thread cberry
Not sure this is a bug in R. Maybe its a bug in my understanding of unz(). The character 'b2' (hexadecimal) is in position 535 of line 1 of 'naughty.csv'. This character appears as superscript '2' and came to me in an EXCEL file that I converted to text in a comma separated ( *.csv ) format.

RE: [R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread Liaw, Andy
> From: Jason Liao [mailto:[EMAIL PROTECTED] > > Thanks for Prof. Ripley and Andy for your technical > explantion. It seems that that the real CPU speed has not > advanced as fast as these Ghz or other performance indicator suggest. > > Yes, my program is totally CPU intensive. > > "> We do

Re: [R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread A.J. Rossini
Jason Liao <[EMAIL PROTECTED]> writes: > R can not really use dual CPU for one R session if I understand > correctly It certainly can, using message passing libraries or sockets. While it isn't technically "one session", it's awfully similar to that, for the user. best, -tony -- A.J. Rossini

Re: [R] error bars in color

2003-07-23 Thread Jerome Asselin
Yes. See ?plot, ?segments, ?lines, and in particular see the help on the "col" option in ?par. HTH, Jerome On July 23, 2003 03:13 pm, Heinrich Kestler wrote: > Hi, > is it possible to generate differently colored error bars in one plot? > > Thx in advance, > Heinrich > > [[alternative HTM

[R] RE: Shapefiles package upload

2003-07-23 Thread Benjamin . STABLER
Sorry, I meant to say ArcGIS does not permit periods in field names. >-Original Message- >From: STABLER Benjamin >Sent: Wednesday, July 23, 2003 3:29 PM >To: [EMAIL PROTECTED] >Subject: Shapefiles package upload > > >I uploaded version 0.3 of the shapefiles package to CRAN >earlier today

[R] Shapefiles package upload

2003-07-23 Thread Benjamin . STABLER
I uploaded version 0.3 of the shapefiles package to CRAN earlier today. Version 0.2 had a bug that omitted the decimal precision in numeric fields so some programs (such as ArcGIS) would not parse the fields correctly. I also added an argument to write.dbf to swap "." with "_" since ArcGIS does no

[R] error bars in color

2003-07-23 Thread Heinrich Kestler
Hi, is it possible to generate differently colored error bars in one plot? Thx in advance, Heinrich [[alternative HTML version deleted]] __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] position available

2003-07-23 Thread Joanne Chory
Bioinformatics Specialist(Full-Time) Biostatistician To work closely with Ph.D. level biologists in the area of quantitative genetics, including array analysis and software development. Strength in computer programming is required. Experience with R/bioconductor, perl, and C. Must be able to exe

Re: [R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread Jason Turner
Jason Liao wrote: I have been using a laptop computer of Pentium III 1.13 Ghz. I heard that AMD's Athlon has excellent floating point capacity. So I bought a Athlon 2200+ laptop yesterday. I expected that new Athlon 2200+ will be twice as fast as the P III 1.13 GB. I ran a R simulation program and

RE: [R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread Jason Liao
Thanks for Prof. Ripley and Andy for your technical explantion. It seems that that the real CPU speed has not advanced as fast as these Ghz or other performance indicator suggest. Yes, my program is totally CPU intensive. "> We do have a dual P4 Xeon 2.4 GHz with 8GB RAM, and jobs run more > tha

RE: [R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread Liaw, Andy
Overall performance depends on a few other things besides CPU clock speed (e.g., RAM speed and size, cache size, disk speed, etc.) Unless your code is spending great majority of the time in the CPU, you should not expect speed-up to be equal to ratio of clock speeds. (Also, as Prof. Ripley pointe

Re: [R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread Prof Brian Ripley
On Wed, 23 Jul 2003, Jason Liao wrote: > I have been using a laptop computer of Pentium III 1.13 Ghz. I heard > that AMD's Athlon has excellent floating point capacity. So I bought a > Athlon 2200+ laptop yesterday. I expected that new Athlon 2200+ will be > twice as fast as the P III 1.13 GB. I r

[R] Dismal R performance of Athlon moble CPU?

2003-07-23 Thread Jason Liao
I have been using a laptop computer of Pentium III 1.13 Ghz. I heard that AMD's Athlon has excellent floating point capacity. So I bought a Athlon 2200+ laptop yesterday. I expected that new Athlon 2200+ will be twice as fast as the P III 1.13 GB. I ran a R simulation program and the new computer i

RE: [R] paste and NAs

2003-07-23 Thread Benjamin . STABLER
I didn't realize that ifelse is vectorized. Thanks for the suggestion. >-Original Message- >From: Thomas Lumley [mailto:[EMAIL PROTECTED] >Sent: Wednesday, July 23, 2003 11:50 AM >To: STABLER Benjamin >Cc: [EMAIL PROTECTED]; GREGOR Brian J >Subject: Re: [R] paste and NAs > > >On Wed, 23 J

Re: [R] Condition indexes and variance inflation factors

2003-07-23 Thread John Fox
Dear Peter and Uwe, I don't have a copy of Belsley's 1991 book here, but I do have Belsley, Kuh, and Welsch, Regression Diagnostics (Wiley, 1980). If my memory is right, the approach is the same: Belsley's collinearity diagnostics are based on a singular-value decomposition of the scaled but un

Re: [R] S3 and S4 classes

2003-07-23 Thread John Chambers
Duncan Murdoch wrote: > > On Wed, 23 Jul 2003 14:53:56 +0200, Laurent Faisnel > <[EMAIL PROTECTED]> wrote : > >Could anyone point me out what's S3-like in > >the following sample and why it is not fully S4-compatible ? > > ># a function that objects of this class have > >perform <- function(.Obje

Re: [R] .ps files in R

2003-07-23 Thread Michael A. Miller
You can include in your latex with \includegraphics[angle=90]{foo.eps} __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] trouble with maps

2003-07-23 Thread Roger Bivand
On Wed, 23 Jul 2003, Deborah Swayne wrote: > Has anyone else seen this behavior from the "maps" package? > > map('state', fill=TRUE) > > results in a lively mix of overlapping polygons inside a > map of the US, but they have no obvious relationship to > state boundaries. (See attached jpeg.

Re: [R] .ps files in R

2003-07-23 Thread Ko-Kang Kevin Wang
On Wed, 23 Jul 2003, Johanna Hardin wrote: > in ghostscript. I've also tried many different rotating commands in LaTeX > (angle in \includegraphics, \rotate, \sideways,...) But, the picture seems > to be unaffected by any of these commands. I find it the rotation won't show up in the DVI file,

Re: [R] .ps files in R

2003-07-23 Thread Prof Brian Ripley
Did you read ?postscript ? You want an encapsulated PS (EPS) file to include in LaTeX, and that page tells you the options required. On Wed, 23 Jul 2003, Johanna Hardin wrote: > I have recently "printed" in R to a postscript file. I'm working on a SSH > without an X terminal. It was fairly aut

Re: [R] .ps files in R

2003-07-23 Thread Jerome Asselin
You can specify the postscript() option horizontal=F. See help on "horizontal" in ?postscript. HTH, Jerome On July 23, 2003 11:42 am, Johanna Hardin wrote: > I have recently "printed" in R to a postscript file. I'm working on a > SSH > > without an X terminal. It was fairly automatic: > > plot

Re: [R] paste and NAs

2003-07-23 Thread Thomas Lumley
On Wed, 23 Jul 2003 [EMAIL PROTECTED] wrote: > I understand how R treats NAs but in the situation below it would be nice to > have a na.skip argument to paste so it does not convert the NAs to "NA"s and > simply skips those elements of the vector for pasting. > > > x > [1] "2.13" "2.3" NA

[R] .ps files in R

2003-07-23 Thread Johanna Hardin
I have recently "printed" in R to a postscript file. I'm working on a SSH without an X terminal. It was fairly automatic: > plot(x,y) > dev.off() And then the default creates a file called Rplots.ps which I can ftp to my laptop and open in Ghostscript. I can see the file, and nothing looks odd

[R] paste and NAs

2003-07-23 Thread Benjamin . STABLER
I understand how R treats NAs but in the situation below it would be nice to have a na.skip argument to paste so it does not convert the NAs to "NA"s and simply skips those elements of the vector for pasting. > x [1] "2.13" "2.3" NA NA "2.83" NA > paste(x, "0", sep="") [1] "2.130"

Re: [R] Curious warning in R for OS X w/Xwindows

2003-07-23 Thread Bear F. Braumoeller
On Wednesday, July 23, 2003, at 05:03 AM, Peter Dalgaard BSA wrote: "Bear F. Braumoeller" <[EMAIL PROTECTED]> writes: As to the other questions, Sys.getenv("TERM") TERM "xterm" Sys.getenv("PAGER") PAGER "/usr/bin/less" options("pager") $pager [1] "/usr/local/lib/R/bin/pager" And l

[R] trouble with maps

2003-07-23 Thread Deborah Swayne
Has anyone else seen this behavior from the "maps" package? map('state', fill=TRUE) results in a lively mix of overlapping polygons inside a map of the US, but they have no obvious relationship to state boundaries. (See attached jpeg.) I reinstalled the maps and mapdata packages from ft

[R] odd behavior with the maps package

2003-07-23 Thread Deborah Swayne
Has anyone else seen this behavior from the "maps" package? map('state', fill=TRUE) results in a lively mix of overlapping polygons inside a map of the US, but they have no obvious relationship to state boundaries. (See attached jpeg.) I reinstalled the maps and mapdata packages from ftp:/

[R] Passing references to data objects into R functions

2003-07-23 Thread David Khabie-Zeitoune
Hi. I have the following question about reading from large data objects from within R functions; I have tried to simplify my problem as much as possible in what follows. Imagine I have various large data objects sitting in my global environment (call them "data1", "data2", ...). I want to write

Re: [R] Spectral Analysis

2003-07-23 Thread zhu wang
Take a look at http://onb.ent.psu.edu/onb1/ Look at the bottom of Software R/S-code. I believe spec.lomb is not a default package/function in Splus, that explains why you can not use spec.lomb unless you install the package "nlt". Hope this helps. __ [

Re: [R] Strange behaviour when running R from within Emacs on Winddows

2003-07-23 Thread Prof Brian Ripley
On Wed, 23 Jul 2003, Søren Højsgaard wrote: > Dear R-experts, > I run R in a shell under Emacs on Win2k using ESS. I get the following strange error > > > shell("copy c:\\file.txt c:\\newfile.txt") > warning: extra args ignored after 'copy' > Forkert syntaks for kommandoen. > Warning message: >

Re: [R] Condition indexes and variance inflation factors

2003-07-23 Thread Uwe Ligges
Peter Flom wrote: Has anyone programmed condition indexes in R? I know that there is a function for variance inflation factors available in the car package; however, Belsley (1991) Conditioning Diagnostics (Wiley) notes that there are several weaknesses of VIFs: e.g. 1) High VIFs are sufficient b

[R] Condition indexes and variance inflation factors

2003-07-23 Thread Peter Flom
Has anyone programmed condition indexes in R? I know that there is a function for variance inflation factors available in the car package; however, Belsley (1991) Conditioning Diagnostics (Wiley) notes that there are several weaknesses of VIFs: e.g. 1) High VIFs are sufficient but not necessary co

Re: [R] S3 and S4 classes

2003-07-23 Thread Duncan Murdoch
On Wed, 23 Jul 2003 14:53:56 +0200, Laurent Faisnel <[EMAIL PROTECTED]> wrote : >Could anyone point me out what's S3-like in >the following sample and why it is not fully S4-compatible ? ># a function that objects of this class have >perform <- function(.Object) UseMethod("perform", .Object); It

Re: [R] animal models and lme

2003-07-23 Thread Spencer Graves
If you can solve the problem for fixed rho = (va/ve) using gls, then you can call gls for many values of rho, plot the log(likelihood) contours vs. rho, construct confidence intervals, etc. You may even be able to write a function to return (-2)*log(likelihood) for a fixed rho and then use

Re: [R] Read trajectory file into R

2003-07-23 Thread Prof Brian Ripley
On Wed, 23 Jul 2003, Karim Elsawy wrote: > > Prof Brian Ripley wrote: > > > > On Wed, 23 Jul 2003, Karim Elsawy wrote: > > > > > I wonder if there is a way to read a molecular dynamic trajectory file ( > > > binary file) produced by CHARMM into R. Something like that in matlab. > > > Actually t

Re: [R] Read trajectory file into R

2003-07-23 Thread Thomas Lumley
On Wed, 23 Jul 2003, Karim Elsawy wrote: > Thanks a lot for your help, actually I do not know the exact file format > at the moment > all what I know is : > "The DCD files (the trajectory files) are single precision binary > FORTRAN files, so are transportable > between computer architectures. The

Re: Re[2]: [R] lattice: how to format axis labels?

2003-07-23 Thread Douglas Bates
High-level control of axes in xyplot is implemented by the scales argument to xyplot. You can include components 'at' and 'labels' in a list given as the scales argument. See ?xyplot. Wladimir Eremeev <[EMAIL PROTECTED]> writes: > jhcc> check out 'sprintf' for formating in a specific way. > >

Re: [R] animal models and lme

2003-07-23 Thread Douglas Bates
<[EMAIL PROTECTED]> writes: > Not convinced that responses so far have addressed the problem. The > model is > > y = mu + U + e > > where e is a vector of independendent errors with variance ve, and U > is a vector of random effects with covariance matrix va*A, where A is a > known matrix (which

[R] S3 and S4 classes

2003-07-23 Thread Laurent Faisnel
Hi helpers, I've been programming in R for a few months now but I still have doubts about my code - I would like it to be completely S4-compatible. The current code works fine but is probably 'unclean'. I read the interesting article in the last R News and it helped me understand the difference

[R] Strange behaviour when running R from within Emacs on Winddows

2003-07-23 Thread Søren Højsgaard
Dear R-experts, I run R in a shell under Emacs on Win2k using ESS. I get the following strange error > shell("copy c:\\file.txt c:\\newfile.txt") warning: extra args ignored after 'copy' Forkert syntaks for kommandoen. Warning message: cmd execution failed with error code 1 in: shell("copy c:\\fi

Re[2]: [R] lattice: how to format axis labels?

2003-07-23 Thread Wladimir Eremeev
Dear james, jhcc> check out 'sprintf' for formating in a specific way. This will not solve the problem. I will have to specify the argument like labels=(...). I would like to avoid it. I wonder if there a key or option to make automatically appearing labels be formatted in the mentioned way. I h

[R] multinomial logit discrete choice model

2003-07-23 Thread DJ
Hi, I'm struggling trying to specify a multinomial logit discrete choice model in R. Any help and/or code examples appreciated. I am specifically interested in specifying a model where no universal choice set exists and each choice set has a variable number of alternatives (one of which is chose

[R] category data (text categories)

2003-07-23 Thread Alan Davis
May I ask for a pointer to documentation of how to use category data (text categories), which I would like to graph as present/absent at stations along a transect. Thanks for any pointers, Alan Davis -- [EMAIL PROTECTED] 1-670-322-6580 Alan E. Davis, P

Re: [R] Read trajectory file into R

2003-07-23 Thread Karim Elsawy
Prof Brian Ripley wrote: > > On Wed, 23 Jul 2003, Karim Elsawy wrote: > > > I wonder if there is a way to read a molecular dynamic trajectory file ( > > binary file) produced by CHARMM into R. Something like that in matlab. > > Actually this will save tremendous effort in post processing. > > I

Re: [R] animal models and lme

2003-07-23 Thread iwhite
Not convinced that responses so far have addressed the problem. The model is y = mu + U + e where e is a vector of independendent errors with variance ve, and U is a vector of random effects with covariance matrix va*A, where A is a known matrix (which we can assume is a correlation matrix). If w

Re: [R] Boosting,bagging and bumping. Questions about R tools and predictions.

2003-07-23 Thread Gavin Simpson
Take a look at the randomForest package on CRAN: randomForest: Breiman's random forest for classification and regression Classification and regression based on a forest of trees using random inputs. Version: 3.9-6 Depends: R (>= 1.7.0) Author: Fortran original by Leo Breiman and Adele Cutler, R

Re: [R] Curious warning in R for OS X w/Xwindows

2003-07-23 Thread Peter Dalgaard BSA
"Bear F. Braumoeller" <[EMAIL PROTECTED]> writes: > I'm starting R with > > xterm -sb -rightbar -sl 1000 -bg black -fg blue -title R -e > /usr/local/bin/R > > -- but it also happens if I just start a vanilla terminal and type R. > > As to the other questions, > > > Sys.getenv("TERM") >

Re: [R] Read trajectory file into R

2003-07-23 Thread Prof Brian Ripley
On Wed, 23 Jul 2003, Karim Elsawy wrote: > I wonder if there is a way to read a molecular dynamic trajectory file ( > binary file) produced by CHARMM into R. Something like that in matlab. > Actually this will save tremendous effort in post processing. If you know the file format, yes. That's a

[R] Read trajectory file into R

2003-07-23 Thread Karim Elsawy
dear helpers, I wonder if there is a way to read a molecular dynamic trajectory file ( binary file) produced by CHARMM into R. Something like that in matlab. Actually this will save tremendous effort in post processing. best regards karim __ [EMAIL PROT

[R] lattice: how to format axis labels?

2003-07-23 Thread Wladimir Eremeev
Dear r-help, I draw graphics with xyplot() function. Labels on the y axis are appearing as follows: "1.5, 1, 0.5, 0" I'd like to have them to be "1.5, 1.0, 0.5, 0.0", i.e. with fixed number of digits after the dot (one in this case). Is there any way to do this without implicit specify

Re: [R] libblas.so.3

2003-07-23 Thread Prof Brian Ripley
On Tue, 22 Jul 2003, Francisco J Molina wrote: > > I am using linux rehat 9 > When I try to run R I get > > usr/lib/R/bin/R.bin: error while loading shared libraries: libblas.so.3: cannot open > shared object file: No such file or directory > > ( because I removed it ) > > I thought of compil

Re: [R] using getBioC()

2003-07-23 Thread Bill Barnard
There was a post in the Bioconductor mailing list regarding the same error message: https://stat.ethz.ch/pipermail/bioconductor/2003-February/000854.html One possible solution mentioned there was to remove the line in getBioC() that checks the capabilities. Since it seems likely that you do have

Re: [R] Processing a large number of files

2003-07-23 Thread Uwe Ligges
Douglas Bates wrote: I maintain the Devore5 package which contains the data sets from the 5th edition of Jay Devore's text "Probability and Statistics for Engineering and the Sciences". The 6th edition has now been published and it includes several new data sets in exercises and examples. In addi