Re: [R] Odd behaviour of step (and stepAIC)?

2004-03-20 Thread Berwin A Turlach
JO == Jari Oksanen [EMAIL PROTECTED] writes: JO There seems to be some obscure features in step() when you JO have interaction terms: their interpretation is order JO sensitive. [...] Just because R internally decides to order JO terms differently than in the scope (this may

Re: [R] substitute question

2004-03-20 Thread Peter Dalgaard
Tony Plate [EMAIL PROTECTED] writes: Earlier in this discussion, Peter Dalgard stated you can only do substitutions on language objects and then used the function is.language() in an example, which I took at that time to imply that substitute() would go inside objects for which is.language()

Re: [R] labels on axis(4) and adaptation of legend to size of the plot

2004-03-20 Thread Andrew Robinson
Jan, try ?mtext Cheers Andrew On Saturday 20 March 2004 06:35, Jan Verbesselt wrote: Hi weekend R helpers, *Is it possible to get a ylabel on the right hand side axis == axis(4)? Thankx, Jan opar - par(mfrow = c(2,1)) plot(ts.Origi[,1],ylab='NDII', main=name) # Add the

Re: [R] labels on axis(4) and adaptation of legend to size of the plot

2004-03-20 Thread Patrick Burns
Two tricks to know: You need to have enough room in that margin, so you probably need to set the mar graphics parameter. For example: par(mar=c(5,4,4,4) +.1) You then want to use mtext, like: mtext(side=4, line=2, the other ylab) Patrick Burns Burns Statistics [EMAIL PROTECTED] +44 (0)20 8525

[R] minimum values

2004-03-20 Thread Rogério Rosa da Silva
Hi, In my data set I have created a matrix of distances. How can I get a list of minimum value in each row? thanks a lot, Rogério __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

RE: [R] minimum values

2004-03-20 Thread John Fox
Dear Rogerio, If the distances is the distance matrix, apply(distances, 1, min) should give you what you want. I hope this helps, John -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Rogério Rosa da Silva Sent: Saturday, March 20, 2004 7:10 AM

Re: [R] minimum values

2004-03-20 Thread Marc Schwartz
On Sat, 2004-03-20 at 06:09, Rogrio Rosa da Silva wrote: Hi, In my data set I have created a matrix of distances. How can I get a list of minimum value in each row? thanks a lot, Rogrio If your matrix is 'mat': apply(mat, 1, min) Example: mat - matrix(rnorm(30), ncol = 3) mat

[R] newbie question on clustering performance

2004-03-20 Thread Xin Gong
May I ask you some questions on clustering analysis? How do I get jacard score from cluster object? or after I call kmanes or pam function how do I know which number of clusters sounds the best. Another question on heatmap, if you do not mind, how do I retrive the name of rows when some area

[R] timing a function

2004-03-20 Thread Fred J.
Hello is there a way to time how long it takes to run a code in R. somthing like tic toc in matlab as such? help.search(timing) put out nothing. and while I got you, debugging the code, is there a step through and the rest of the debugging tools working with ESS. thanks

Re: [R] timing a function

2004-03-20 Thread Marc Schwartz
On Sat, 2004-03-20 at 12:11, Fred J. wrote: Hello is there a way to time how long it takes to run a code in R. somthing like tic toc in matlab as such? help.search(timing) put out nothing. See ?system.time and while I got you, debugging the code, is there a step through and the rest of the

[R] contrast lme and glmmPQL and getting additional results...

2004-03-20 Thread Paul Johnson
I have a longitudinal data analysis project. There are 10 observations on each of 15 units, and I'm estimating this with randomly varying intercepts along with an AR1 correction for the error terms within units. There is no correlation across units. Blundering around in R for a long time, I

[R] Operating on windows of data

2004-03-20 Thread Sean Davis
I have a data set that is comprised of, for simplicity, a vector of numbers that I want to march across in overlapping windows of say 10 values each, computing a couple of values for each window. Is there a vectorized way to do this, or do I truly need to resort to looping--I think so? Any

RE: [R] Operating on windows of data

2004-03-20 Thread Gabor Grothendieck
See ?embed and possibly ?filter Date: Sat, 20 Mar 2004 15:16:38 -0500 From: Sean Davis [EMAIL PROTECTED] To: R-Help [EMAIL PROTECTED] Subject: [R] Operating on windows of data I have a data set that is comprised of, for simplicity, a vector of numbers that I want to

Re: [R] Operating on windows of data

2004-03-20 Thread Marc Schwartz
On Sat, 2004-03-20 at 14:16, Sean Davis wrote: I have a data set that is comprised of, for simplicity, a vector of numbers that I want to march across in overlapping windows of say 10 values each, computing a couple of values for each window. Is there a vectorized way to do this, or do I

Re: [R] contrast lme and glmmPQL and getting additional results...

2004-03-20 Thread kjetil
On 20 Mar 2004 at 11:06, Paul Johnson wrote: From what you say, it seems like you have a linear normal (mixed) model. This is what lme is made for, and there is no reason to use glmmPQL (which calls lme iteratively). Kjetil Halvorsen I have a longitudinal data analysis project. There are 10

[R] setwd() permenent

2004-03-20 Thread Fred J.
Hello when I close and reopen R it gives back a different wd than what I used in setwd. how can I get it to permenently use dirname in setwd(dirname) and not the other name it keeps defaulting to? I am using W2K and ESS thanks __ [EMAIL PROTECTED]