Re: [R] Scanning tab-separated numbers

2004-03-02 Thread Prof Brian Ripley
On Mon, 1 Mar 2004, Robert W. Baer, Ph.D. wrote: The problem would seem to be that the Windows GUI does not accept the TAB key as input. There is no whitespace between your values to be recongnized. Take your sequesnce and save it as a tab separated file in Excel or through Notepad, and you

[R] filled contour

2004-03-02 Thread Abdou Ali
Hello, how can I fill map contour with R without setting the legend. I cannot see an option to skip the legend for filled.contour. Thank for your response. __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE

Re: [R] passing a string from .C()

2004-03-02 Thread Prof Brian Ripley
On Mon, 1 Mar 2004, Vadim Ogranovich wrote: Could someone please point to an example of passing strings from .C() calls back to R? There are not many, as such things are better done using .Call. I want to be able to do something like this: str - .C(return_foo_string,

[R] Problem with Integrate

2004-03-02 Thread Anon.
The background: I'm trying to fit a Poisson-lognormal distrbutuion to some data. This is a way of modelling species abundances: N ~ Pois(lam) log(lam) ~ N(mu, sigma2) The number of individuals are Poisson distributed with an abundance drawn from a log-normal distrbution. To fit this to data, I

[R] Import range of cells from Excel

2004-03-02 Thread Jean-Noel
Dear all, I would like to import a range of cells (e.g. F10:K234) from an Excel worksheet to R. I have looked for documentation on RODBC and RDCOMClient but I was not able to find enough information to solve my problem and all the examples I have seen were dealing with an entire worksheet, not a

Re: [R] non-negative least-squares

2004-03-02 Thread Simon Wood
I think pcls will do what you want. You'll need to explicitly set up the constraint matrix, since pcls is designed to deal with general linear inequality constraints (i.e. it will do more than just non-negativity). Also the penalty matrix/matrices for pcls would be t(L)%*%L, I think. btw, there's

Re: [R] boxplot notches

2004-03-02 Thread Christoph Scherber
Dear colleagues, I think it would be a good idea to include a short note in the R boxplot() help file, stating exactly how the confidence levels are calculated (the notches are +/- 1.58 IQR/sqrt(n)) - at least as a guidance for users not advanced enough to directly interpret the code. Would

Re: [R] Problem with Integrate

2004-03-02 Thread Spencer Graves
Have you done a search of www.r-project.org - search - R site search for hermite quadrature? I just got 11 hits on this, the second of which referred to gauss.quad {statmod}. This said, See also ... gauss.quad.prob {statmod}. I haven't tried it, but it looks like what you want.

[R] R programming

2004-03-02 Thread Andersson, Henrik
I have been using R for a few months to plot my data, and fit statistical models and so on. It is stated that R is not only a package for statistcs and graphics but also a programming language. Currently I am working with Fortran 90 to do numerical simulations (1D reactive -transport), and

Re: [R] R programming

2004-03-02 Thread Sean O'Riordain
Hi Henrick, if you checked the posting guide (refer the bottom of all emails) you'd have spotted a link # R Language Definition (aka `R Language Manual') describes the R language, data objects, etc. http://cran.r-project.org/doc/manuals/R-lang.pdf cheers, Sean ps. R is indeed a full and

[R] dynamic linking

2004-03-02 Thread Samuel Kemp
Hi, I managed to fixed my problem with linking a .a file (an archive of .o files) with the following. g++ -I/usr/local/lib/R/include -I/usr/local/include -mieee-fp -fPIC -g -O2 -c fileGT.cc -o fileGT.o (i.e. the beginning bit of R CMD SHLIB) then I linked (underlined) the files

Re: [R] boxplot notches

2004-03-02 Thread Prof Brian Ripley
On Mon, 1 Mar 2004, David James wrote: Prof Brian Ripley wrote: On Mon, 1 Mar 2004, Martin Maechler wrote: TL == Thomas Lumley [EMAIL PROTECTED] on Mon, 1 Mar 2004 09:54:48 -0800 (PST) writes: TL On Mon, 1 Mar 2004, Christoph Scherber wrote: Dear list members,

Re: [R] boxplot notches

2004-03-02 Thread kjetil
On 1 Mar 2004 at 9:54, Thomas Lumley wrote: On Mon, 1 Mar 2004, Christoph Scherber wrote: Dear list members, Can anyone tell me how the notches in boxplot(Y~X,notch=T) are calculated? What do these notches represent exactly? I´d suppose they are Conficence Intervals for the median,

Re: [R] filled contour

2004-03-02 Thread Petr Pikal
Hi On 2 Mar 2004 at 9:11, Abdou Ali wrote: Hello, how can I fill map contour with R without setting the legend. I cannot does image() with added contour() do what you want? Cheers Petr see an option to skip the legend for filled.contour. Thank for your response.

Re: [R] Import range of cells from Excel

2004-03-02 Thread Petr Pikal
Hi What about to copy a range, paste it to notepad, save it as a *.txt file and read with read.delim(). Cheers Petr On 2 Mar 2004 at 10:10, Jean-Noel wrote: Dear all, I would like to import a range of cells (e.g. F10:K234) from an Excel worksheet to R. I have looked for documentation on

re: [R] boxplot notches

2004-03-02 Thread Michael Friendly
I think John Tukey's idea was that this formula (or just the fact of using median and quartiles) is still often approximately correct for quite a few kinds of moderate contaminations... It may be approximately correct for the width of a CI (and when I checked it was only appproximately

Re: [R] superimposing two scatterplots

2004-03-02 Thread Robert W. Baer, Ph.D.
Look at: ?plot ?points - Original Message - From: Raphael Schoenle [EMAIL PROTECTED] To: [EMAIL PROTECTED] Sent: Monday, March 01, 2004 4:35 PM Subject: [R] superimposing two scatterplots Hi, How can I plot two scatterplots on the same, one panel? I have two times series (price data

Re: [R] boxplot notches

2004-03-02 Thread Christoph Scherber
In McGill et al. (1978) there´s a description of the calculation as follows (p. 16): The widths [are] computed from the midspread or interquartile range (R) of the data (...), and the number of observations (N) for each group. The Gaussian-based asymptotic approximation (Kendall and Stuart

Re: [R] boxplot notches

2004-03-02 Thread P. B. Pynsent
A Google search showed that all this was discussed in April 1988 with an extensive reply to the question from M Maechler. I, as a non-statistician, blindly believed what was written in the boxplot() help file, I am sure many would be grateful to this help being modified. I still do not

Re: [R] Problem with Integrate

2004-03-02 Thread Thomas Lumley
On Tue, 2 Mar 2004, Anon. wrote: The background: I'm trying to fit a Poisson-lognormal distrbutuion to some data. This is a way of modelling species abundances: N ~ Pois(lam) log(lam) ~ N(mu, sigma2) The number of individuals are Poisson distributed with an abundance drawn from a

Re: [R] dev.print and X11(canvas = black)

2004-03-02 Thread Jeffrey J. Hallman
Prof Brian Ripley [EMAIL PROTECTED] wrote: On Mon, 1 Mar 2004 [EMAIL PROTECTED] wrote: In Splus, I often use graphics windows with a black background and white foreground. The S print.graph() function sends the current plot to my printer but with a white background and black

Re: [R] Problem with Integrate

2004-03-02 Thread Anon.
Thomas Lumley wrote: On Tue, 2 Mar 2004, Anon. wrote: The background: I'm trying to fit a Poisson-lognormal distrbutuion to some data. This is a way of modelling species abundances: N ~ Pois(lam) log(lam) ~ N(mu, sigma2) The number of individuals are Poisson distributed with an abundance drawn

Re: [R] boxplot notches

2004-03-02 Thread Thomas Lumley
On Tue, 2 Mar 2004, P. B. Pynsent wrote: A Google search showed that all this was discussed in April 1988 with an extensive reply to the question from M Maechler. I, as a non-statistician, blindly believed what was written in the boxplot() help file, I am sure many would be grateful to this

Re: [R] Problem with Integrate

2004-03-02 Thread Thomas Lumley
On Tue, 2 Mar 2004, Anon. wrote: Thomas Lumley wrote: The help page for integrate() says When integrating over infinite intervals do so explicitly, rather than just using a large number as the endpoint. This increases the chance of a correct answer - any function whose

Re: [R] boxplot notches

2004-03-02 Thread P. B. Pynsent
A Google search showed that all this was discussed in April 1988 with an extensive reply to the question from M Maechler. Sorry above should be 1998 On 2 Mar 2004, at 10:18, Christoph Scherber wrote: Dear colleagues, I think it would be a good idea to include a short note in the R boxplot()

re: [R] boxplot notches

2004-03-02 Thread Prof Brian Ripley
On Tue, 2 Mar 2004, Michael Friendly wrote: I think John Tukey's idea was that this formula (or just the fact of using median and quartiles) is still often approximately correct for quite a few kinds of moderate contaminations... It may be approximately correct for the width

Re: [R] dev.print and X11(canvas =

2004-03-02 Thread Gabor Grothendieck
Try this: dev.control(displaylist=enable) plot(1:4) identify(1:4) # perform some interactions which label points myplot - recordPlot() plot(1:10) # overwrite first plot myplot # at this point the original plot is back with your annotations --- Date: Tue, 02 Mar 2004 10:17:39 -0500 From:

Re: [R] boxplot notches

2004-03-02 Thread Martin Maechler
P == P B Pynsent [EMAIL PROTECTED] on Tue, 2 Mar 2004 15:16:23 + writes: P A Google search showed that all this was discussed in April 1998 with P an extensive reply to the question from M Maechler. Yes, indeed: http://finzi.psych.upenn.edu/R/Rhelpold/archive/0839.html ,

[R] Jackknife after bootstrap influence values in boot package?

2004-03-02 Thread Samuelson, Frank*
Is there a routine in the boot package to get the jackknife-after- bootstrap influence values? That is, the influence values of a jackknife of the bootstrap estimates? I can see how one would go about it from the jack.after.boot code, but that routine only makes pretty pictures. It wouldn't be

[R] how to delete a matrix column

2004-03-02 Thread Aimin Yan
Hello, I am new to R, How to delete a matrix column. Thanks, __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Re: [R] dev.print and X11(canvas = black)

2004-03-02 Thread Prof Brian Ripley
On Tue, 2 Mar 2004, Jeffrey J. Hallman wrote: Jim Lemon made essentially the same suggestion: run the code that created the plot twice, once with X11 as the active device, and again with postscript as the active device. The trouble with this is that is doesn't handle plots that the user has

RE: [R] how to delete a matrix column

2004-03-02 Thread Harold Doran
You can use dataset$columnname=NULL HCD -- Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington Street, Suite 220 Alexandria, Virginia 22314 703.647.1628 -Original Message- From: Aimin Yan [mailto:[EMAIL PROTECTED] Sent: Tuesday, March

[R] an extension to 'array'

2004-03-02 Thread Heywood, Giles
1. I want to extend the 'array' class, and prefer to use S4 in the belief that this is the best structure for new projects (as the documentation says). I actually wish to do something similar to the excellent Oarray by Jonathan Rougier, but as this class is S3, I can't see how to extend it by S4.

[R] row.names are dropped when extracting one column ?

2004-03-02 Thread ryszard . czerminski
Apparently row names are dropped when I extract single column from a data frame. Why this behaviour ? y - as.matrix(df[,1:2]); length(row.names(y)) [1] 324 y - as.matrix(df[,1:1]); length(row.names(y)) [1] 0 Best regards, Ryszard __ [EMAIL

Re: [R] how to delete a matrix column

2004-03-02 Thread Rolf Turner
Aimin Yan wrote: I am new to R, How to delete a matrix column. e.g. y - x[,-5] deletes column 5 I was going to point out to the inquirer where this feature is documented --- but I can't bleedin' find it!!! Can anyone tell me? I tried ?[ but got no joy that I could discern.

Re: [R] row.names are dropped when extracting one column ?

2004-03-02 Thread Prof Brian Ripley
On Tue, 2 Mar 2004 [EMAIL PROTECTED] wrote: Apparently row names are dropped when I extract single column from a data frame. Why this behaviour ? y - as.matrix(df[,1:2]); length(row.names(y)) [1] 324 y - as.matrix(df[,1:1]); length(row.names(y)) [1] 0 Why are you converting a subsetted

Re: [R] row.names are dropped when extracting one column ?

2004-03-02 Thread Marc Schwartz
On Tue, 2004-03-02 at 11:59, [EMAIL PROTECTED] wrote: Apparently row names are dropped when I extract single column from a data frame. Why this behaviour ? y - as.matrix(df[,1:2]); length(row.names(y)) [1] 324 y - as.matrix(df[,1:1]); length(row.names(y)) [1] 0 Best regards,

[R] gls anova wald test calculations

2004-03-02 Thread Alan Swanson
I have a question about the Wald test F-statistics that are calculated when the anova() command is used on a singular gls or lme object. As I recall from my linear models class, the Wald test examines H0: C'B = d0 vs Ha: C'B != d0. Does anybody know how this C matrix is constructed in R? Is

Re: [R] how to delete a matrix column

2004-03-02 Thread Prof Brian Ripley
On Tue, 2 Mar 2004, Rolf Turner wrote: Aimin Yan wrote: I am new to R, How to delete a matrix column. e.g. y - x[,-5] deletes column 5 I was going to point out to the inquirer where this feature is documented --- but I can't bleedin' find it!!! Can anyone tell me? I tried ?[ but

[R] possible bug in aov?

2004-03-02 Thread Arditi, Aries
Hi, I'm interested in doing a repeated measures anova using aov. The procedure is nicely described in section 6.7.1, pp. 24-27 of Baron and Li's Notes on the use of R for psychology experiments and questionnaires, and I've reproduced their example exactly. My own problem is almost identical to

[R] error() and C++ destructors

2004-03-02 Thread Vadim Ogranovich
Hi, I am writing C++ functions that are to be called via .Call() interface. I'd been using error() (from R.h) to return to R if there is an error, but then I realized that this might be not safe as supposedly error() doesn't throw an exception and therefore some destructors do not get called and

Re: [R] boxplot notches

2004-03-02 Thread Frank E Harrell Jr
On Tue, 2 Mar 2004 07:24:41 -0800 (PST) Thomas Lumley [EMAIL PROTECTED] wrote: On Tue, 2 Mar 2004, P. B. Pynsent wrote: A Google search showed that all this was discussed in April 1988 with an extensive reply to the question from M Maechler. I, as a non-statistician, blindly believed

Re: [R] dev.print and X11(canvas = black)

2004-03-02 Thread Paul Murrell
Hi Prof Brian Ripley wrote: On Tue, 2 Mar 2004, Jeffrey J. Hallman wrote: Jim Lemon made essentially the same suggestion: run the code that created the plot twice, once with X11 as the active device, and again with postscript as the active device. The trouble with this is that is doesn't

[R] time series with diff. y axes

2004-03-02 Thread Erin Hodgess
Dear R People: Someone had asked a question about plotting 2 time series with different y-axes. There is a nice example in S+, but here is a toy example (modified from the S+ example) which someone may find useful: y1 - rnorm(48) y2 - rnorm(48) y1.ts - ts(y1,start=c(1996,1),freq=12) y2.ts -

RE: [R] using object reference

2004-03-02 Thread Rajarshi Guha
On Tue, 2004-03-02 at 17:16, Andy Bunn wrote: Don't you just need to use return? x - 1 attr(x,'a') - 'some text' f - function(z) { attr(z,'a') - 'some new text' return(z) } y - f(x) y Well what I'm trying to is to pop up a Tk window and allow the user to add some text which gets

[R] Problem untarring R-1.8.1.tgz

2004-03-02 Thread Alistair Campbell
Hi, I have been trying to untar the source file for R and having no success. I have used decompression utilities that handle .tar and .tgz file but I keep getting an error that there is a Header CRC Error. So, I am wondering if this is just me and I am not using the right decompression

[R] Re: Untarring R-1.8.1.tgz

2004-03-02 Thread Alistair Campbell
Hi, Disregard my post. I used another utility and it worked perfectly. The others are obviously not up to the task! Cheers Alistair Campbell __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] Problem untarring R-1.8.1.tgz

2004-03-02 Thread Peter Dalgaard
Alistair Campbell [EMAIL PROTECTED] writes: Hi, I have been trying to untar the source file for R and having no success. I have used decompression utilities that handle .tar and .tgz file but I keep getting an error that there is a Header CRC Error. So, I am wondering if this is just

[R] R2HTML adding caption to data-frame

2004-03-02 Thread Matthew.Kelly
Hi, I was wondering if there was an easy way to add a caption to a data frame when it is being sent to a html page using R2HTML. currently the following command will produce a reasonably formatted table in a html file HTML(mydataframe) It would be good if you could send caption with the table

RE: [R] using object reference

2004-03-02 Thread Gabor Grothendieck
The double arrow in the example you cite is essential; you can't leave it out since that's what tells R to search through its parents' environments: attr(z,'a') - 'some new text' --- Date: Tue, 02 Mar 2004 16:52:20 -0500 From: Rajarshi Guha [EMAIL PROTECTED] To: R [EMAIL PROTECTED]

[R] Find out the day of week for a chron object?

2004-03-02 Thread Ajay Shah
I know that this is correct: library(chron) x = dates(01-03-04, format=d-m-y, out.format=day mon year) print(x) It gives me the string 01 Mar 2004 which is correct. I also know that I can say: print(day.of.week(3,1,2004)) in which case he says 1, for today is monday. My

[R] Stuck in trying to convert repetitive code into a function

2004-03-02 Thread Ajay Shah
Folks, I have the following repetitive code which works correctly: A = read.table(file=junior.csv, sep=,, col.names=c(date, l)); A$date = chron(as.character(A$date), format=m/d/y); r.junior = levels2weeklyret(lastfriday, A$date, A$l); plot(A$date, A$l, type=l, col=red, main=Junior

Re: [R] Nonparametric test of randomness (Run Test)

2004-03-02 Thread Ajay Shah
Suppose you have a returns vector r, which contains values like +3, -2, -7, +3, etc. Then say: require(tseries); test = runs.test(factor(sign(r))) print(test) If you want to access the prob value of the test from within a program, this is found in test$p.value -- Ajay Shah

Re: [R] Find out the day of week for a chron object?

2004-03-02 Thread Dirk Eddelbuettel
On Mon, Mar 01, 2004 at 10:18:33PM +0530, Ajay Shah wrote: I know that this is correct: library(chron) x = dates(01-03-04, format=d-m-y, out.format=day mon year) print(x) It gives me the string 01 Mar 2004 which is correct. I also know that I can say:

[R] Confusion about coxph and Helmert contrasts

2004-03-02 Thread Christian Lederer
Hi, perhaps this is a stupid question, but i need some help about Helmert contrasts in the Cox model. I have a survival data frame with an unordered factor `group' with levels 0 ... 5. Calculating the Cox model with Helmert contrasts, i expected that the first coefficient would be the same as if

[R] why is rcmd shlib including files not in the list?

2004-03-02 Thread Jeff D. Hamann
I've been trying to create a shlib file (using mingw instead of the cygwin -- after reconfiguring my machine, yet again) and I'm getting the following behaviour that I think is incorrect: When using rcmd shlib and explicitly listing the files to include in the build such as: