[R] variable values in plotmath expressions

I tried: t - sample text x - 333 plot( 1, 1 , main=expression( main[x1] ) , xlab=expression( paste( t, xlab[x2] ) ) , sub=parse( text = paste( sub[, x, ] ) ) ) The displayed labels for ``main'' and ``sub'' are as expected. But ``xlab'' shows only t

Re: [R] variable values in plotmath expressions

?substitute, as in plot(1, 1, main = expression(main[x1]), xlab = substitute(paste(t, xlab[x2]), list(t=t)), sub = parse(text = paste(sub[, x, ])) ) [Why do you think ( needs a space but = does not? There is a canonical way to display R code: see `Writing R Extensions' --

RE: [R] Problems with heatmap.2

Hi I'm afraid not! heatmap.2(as.matrix(d),Rowv=as.dendrogram(hc.gene),Colv=1:4,scale=row, trace=none,col=greenred.colors(79)) Produces exactly the same results as if Colv=FALSE - my columns get re-ordered (by a dendrogram no less). Clearly this isn't meant to be happening (by the way,

RE: [R] Problems with heatmap.2

Ploughing on and solving my own problem, I find that setting: dendrogram=row In heatmap.2 swicthes off column ordering. Except that column ordering should have been switched off anyway, by Colv=FALSE, shouldn't it? Oh and dendrogram is spelled dendogram in the help :-) Now all I need to know

Re: [R] newsgroup on R

At 08:51 02/09/2004, you wrote: HI, by newsgroup, I mean a kind of user forum like MATLABs http://newsreader.mathworks.com/WebX?14@@/comp.soft-sys.matlab where anyone can post a question, without receiving unnecessary emails daily Hi Bin Jiang! I visit the url:

Re: [R] new user and non parametric test

Il ven, 2004-09-03 alle 11:32, Christoph Lange ha scritto: etc. - it's true. I didn't find an easy way to find certain test in R either. But e.g. apropos(wilcox) thanks a lot, i'll seach also for the book here in the biblio, the simplest way to load an .xls file ? export to txt ? thanks

Re: [R] new user and non parametric test

(Reply to Paolo Ariano) Paolo, [...] thanks a lot, i'll seach also for the book here in the biblio, the simplest way to load an .xls file ? export to txt ? to .cvs and read into R with read.table(...) -cl ps. BTW: couln't reply to your address, perhaps you fix the problem with

[R] seq

Hi everyone, I've tried the below on R 1.9.1 and the 2004-08-30 builds of R 1.9.1 Patched and R 2.0.0 on Windows 2000, and the results are consistent. seq(0.5, 0, by = -0.1) [1] 0.5 0.4 0.3 0.2 0.1 0.0 seq(0.7, 0, by = -0.1) [1] 7.00e-01 6.00e-01 5.00e-01 4.00e-01

Re: [R] seq

On Fri, 3 Sep 2004, Henric Nilsson wrote: Hi everyone, I've tried the below on R 1.9.1 and the 2004-08-30 builds of R 1.9.1 Patched and R 2.0.0 on Windows 2000, and the results are consistent. seq(0.5, 0, by = -0.1) [1] 0.5 0.4 0.3 0.2 0.1 0.0 seq(0.7, 0, by = -0.1) [1]

Re: [R] seq

On Fri, 3 Sep 2004, Roger Bivand wrote: On Fri, 3 Sep 2004, Henric Nilsson wrote: Hi everyone, I've tried the below on R 1.9.1 and the 2004-08-30 builds of R 1.9.1 Patched and R 2.0.0 on Windows 2000, and the results are consistent. seq(0.5, 0, by = -0.1) [1] 0.5 0.4 0.3 0.2

[R] Standard correlation

Hi Is there a function for computing the standard correlation coefficient (not pearson) in R? Thanks Mick __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide!

Re: [R] confidence intervals

Robert Waters wrote: Dear R users; Im working with lme and Id like to have an idea of how can I get CI for the predictions made with the model. Im not a stats guy but, if Im not wrong, the CIs should be different if Im predicting a new data point or a new group. Ive been searching through the web

[R] strptime problems

Hi, I'm experiencing a problem with strptime. (R 1.9.1 on a Win2000 machine) I have a large list containing 6 columns and 161800 rows. One column contains dates that I want to convert in order to compare the different dates. Some dates work just fine while others become NA. I don't see any

Re: [R] seq

Thanks to all who replied. I'm not surprised that R uses floating point arithmetics, nor am I surprised that an object and it's printed version may differ due to e.g. rounding taking place in the print() function. I am, however, surprised of the fact that seq() prints it's internal

RE: [R] Standard correlation

OK I better clarify what I mean as it appears it may not be a standard test. The pearson correlation coefficient, in laymans terms, uses the shape of a curve around that curve's average to compare two curves. The standard correlation coefficient measures the shape of a curve around zero, and

Re: [R] strptime problems

I suspect rather a locale problem which Date and chron will also have. You are looking for an abbreviated month name (%b). What language is the object ccc in? What language is your computer in? See the example in ?strptime ## read in date info in format 'ddmmm' ## This will

[R] Changing the value of an object's slot in a method

Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: 7bit Received-SPF: none (hypatia: domain of [EMAIL PROTECTED] does not designate permitted sender hosts) X-Virus-Scanned: by amavisd-new at stat.math.ethz.ch X-Spam-Checker-Version: SpamAssassin 2.63 (2004-01-11) on

RE: [R] Standard correlation

Hi! ?scale ?var ?cov /E *** REPLY SEPARATOR *** On 9/3/2004 at 2:20 PM michael watson (IAH-C) wrote: OK I better clarify what I mean as it appears it may not be a standard test. The pearson correlation coefficient, in laymans terms, uses the shape of a curve around that

Re: [R] seq

Henric Nilsson wrote: I am, however, surprised of the fact that seq() prints it's internal representation of 0 as 0.0 when `from = 0.5', but as -1.110223e-16 when `from = 0.7'. These second one isnt a zero! seq(from=0.5,by=-.1,len=6)==0 [1] FALSE FALSE FALSE FALSE FALSE TRUE - in this case

[R] 64 bit benefit for dual G5 powermac

Hi, Can dual G5 Powermac have benefit of 64 bit in terms of memory and maximum matrix size? I have read several positive postings on 64 bit machine. If yes, what will be configure looked like other than { ./configure --with-blas='-framework vecLib' --with-lapack --with-aqua } ? Thank you very much

[R] Different Index behaviors of Array and Matrix

Dear all, I found a difference between the indexing of an array and that of a matrix when there are NA's in the index array. The screen copy is as follows. A - array(NA, dim=6) A [1] NA NA NA NA NA NA idx - c(1,NA,NA,4,5,6) B - c(10,20,30,40,50,60) A[idx] - B A [1] 10 NA NA 40 50 60

[R] arima et graphique

bonjour, je rencontre quelques soucis au niveau de l'utilisation des fonctions arima.forecast/predict; en effet elles sont dites inconnues alors que j'ai bien installé et chargé le package ts. Aussi, j'aimerai savoir comment visualiser le graphique des prévisions avec arima et celui des

[R] arima et graphique

bonjour, je rencontre quelques soucis au niveau de l'utilisation des fonctions arima.forecast/predict; en effet elles sont dites inconnues alors que j'ai bien installé et chargé le package ts. Aussi, j'aimerai savoir comment visualiser le graphique des prévisions avec arima et celui des

[R] Problem with seq.dates in chron

I get faulty output from seq.dates() if I specify a length that is too long. For example, I ask for 129 months in the following call to the function, but it returns 131: R.version.string [1] R version 1.9.1, 2004-06-21 startdatetime - chron(dates=01/01/1995, times=00:00:00)

Re: [R] strptime problems

Prof Brian Ripley ripley at stats.ox.ac.uk writes: : : I suspect rather a locale problem which Date and chron will also have. I don't think chron uses locales so it would likely also be a problem with Date but not with chron. In this example, we change the locale to Danish and then display

Re: [R] allocating memory in C, not in R

Sigal Blay wrote: Thanks Sundar - but how do I return several vector objects of different modes wrapped in a list object using .Call? Sigal Have you looked at the documentation as I suggested? It's there in Section 4.8 of Writing R Extensions. Take a look at the examples for lapply or look at

Re: [R] arima et graphique

The language of the list is English. R no longer has a package `ts', so please read the posting guide and FAQ and get a current version of R. ?UKDriverDeaths and MASS/script/ch14.R will show you examples. On Fri, 3 Sep 2004, sokhna ndao wrote: bonjour, je rencontre quelques soucis au niveau

Re: [R] Problem with seq.dates in chron

I think we just discussed this 2 days ago: https://www.stat.math.ethz.ch/pipermail/r-help/2004-September/055323.html Waichler, Scott R Scott.Waichler at pnl.gov writes: : I get faulty output from seq.dates() if I specify a length that is too : long. : For example, I ask for 129 months in

Re: [R] windowing strategies

jmoreira at fe.up.pt writes: : Does anyone has implemented a function for evaluating models using windowing : strategies, such as growing window or sliding window ones? : The aim is to evaluate regression models on a time series data. I do not use : cross-validation once data sorted in a radom

RE: [R] seq

Note however the result of 0.1*(5:0) == 0.3) [1] FALSE FALSE FALSE FALSE FALSE FALSE as opposed to ==0.5, ==0.4, ==0.2, ==0.1 or ==0! I hope this helps, Ted. E-Mail: (Ted Harding) [EMAIL PROTECTED]

[R] how to debug a sudden exit in non-interactive mode

Hi, I have a piece of R code that calls mgcv::gam. The code runs fine in the interactive mode, but terminates R w/o a single message when run non-interactively. Though I think I should be able to locate the problem by brute force I'd appreciate an advise how to do it more intelligently using R

Re: [R] how to debug a sudden exit in non-interactive mode

If you look in ?debugger, especially right at the end, you will see how to attempt a post-mortem dump of the frams to a file. If that gets run (and this might be a fatal error from C level, but it is not reporting so), you can reload the dump and use debugger on it. On Fri, 3 Sep 2004, Vadim

[R] contrast/test statements in polynomial mixed model

__ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

[R] image() with color key?

Is there an easy way to plot a color key next to a color image (with image() in graphics package)? The color key should also include a numerical scale, so that the colors can be cross-referenced with image intensity values. I see that levelplot has a facility for color keys, but with image it

Re: [R] image() with color key?

On Friday 03 September 2004 14:42, Michael Jerosch-Herold wrote: Is there an easy way to plot a color key next to a color image (with image() in graphics package)? The color key should also include a numerical scale, so that the colors can be cross-referenced with image intensity values. I see

RE: [R] image() with color key?

In some applications, e.g., visualizing the values of the wells in a 96 well (8 x 12) assay plate, the matrix of results is truly discrete and you want to visualize it that way. image() does this, but not filled.contour() AFAIK. If that's correct, the only way I know how to do it is to overlay two

[R] debugging an S4 method

Does anyone know how to use the equivalent of debug() on an S4 method? I would like R to enter the browser not for the generic function, but for the method of the class that I specify. Thanks, David _ David Bickel http://davidbickel.com Research Scientist Pioneer

Re: [R] image() with color key?

Hi, Michael Did you try 'image.plot' in the 'fields' package? Christophe - Original Message - From: Michael Jerosch-Herold [EMAIL PROTECTED] To: [EMAIL PROTECTED] Sent: Friday, September 03, 2004 9:42 PM Subject: [R] image() with color key? Is there an easy way to plot a color key

[R] ML vs. REML with gls()

Hello listmembers, I've been thinking of using gls in the nlme package to test for serial correlation in my data set. I've simulated a sample data set and have found a large discrepancy in the results I get when using the default method REML vs. ML. The data set involves a response that is

Re: [R] debugging an S4 method

On Fri, 2004-09-03 at 21:54, Bickel, David wrote: Does anyone know how to use the equivalent of debug() on an S4 method? I would like R to enter the browser not for the generic function, but for the method of the class that I specify. There's non (that I know) ... you can insert the

[R] R question

Hi, Would you help me solve the following question? Thanks. Question: If I try to set the probability=0.05 and find the approximate x. (The answer should be somewhere between 2.1782 and 2.1783.) I write about this R program as follows but I don¡¦t know how to get the value of x which is