Re: [R] Bug in stepAIC?

2006-10-12 Thread Martin C. Martin
new to R and you don't really know where the problem is. On Wed, 11 Oct 2006, Martin C. Martin wrote: Hi, First of all, thanks for the great work on R in general, and MASS in particular. It's been a life saver for me many times. However, I think I've discovered a bug. It seems

[R] Bug in stepAIC?

2006-10-11 Thread Martin C. Martin
Hi, First of all, thanks for the great work on R in general, and MASS in particular. It's been a life saver for me many times. However, I think I've discovered a bug. It seems that, when I use weights during an initial least-squares regression fit, and later try to add terms using

Re: [R] read.table() and scientific notation

2006-10-10 Thread Martin C. Martin
I think the colClasses argument to read.table() is what you need. Either that, or explicitly cast columns in the data.frame that's returned by read.table(). That's how you get data types that aren't directly supported by read.table(), like various date formats. - Martin January Weiner wrote:

[R] update.default evaluating in wrong environment?

2006-10-10 Thread Martin C. Martin
Hi all, update.default, which is the method used to update lm objects (among others), extracts the call element from it's first argument, updates it, then evaluates it in the parent.frame(). Shouldn't it be evaluated in environment(formula(object)), if that's non-NULL? I ask because I call lm

Re: [R] update.default evaluating in wrong environment?

2006-10-10 Thread Martin C. Martin
be better default behavior than the formula's environment. - Martin On 10/10/06, Martin C. Martin [EMAIL PROTECTED] wrote: Hi all, update.default, which is the method used to update lm objects (among others), extracts the call element from it's first argument, updates it, then evaluates

Re: [R] update.default evaluating in wrong environment?

2006-10-10 Thread Martin C. Martin
Great; thanks for the detailed explanation! - Martin Thomas Lumley wrote: On Tue, 10 Oct 2006, Martin C. Martin wrote: Thanks, or even just: e - environment(formula(f.lm)) But this was more of a bug report. Is update.default wrong? Should it be changed? I don't see how evaluating

[R] Stopping ctrl-\ from qutting R

2006-10-06 Thread Martin C. Martin
Hi, In the Linux (FC3) version of R, ctrl-\ quits R. This wouldn't be so bad, but on my keyboard, it's right next to ctrl-p and I tend to hit it by accident. Is there any way to turn that off? Best, Martin __ R-help@stat.math.ethz.ch mailing list

Re: [R] Stopping ctrl-\ from qutting R

2006-10-06 Thread Martin C. Martin
previous command. Very helpful when you want to repeat something you recently did, perhaps after modifying it a bit. - Martin Alberto Monteiro wrote: Martin C. Martin wrote: In the Linux (FC3) version of R, ctrl-\ quits R. This wouldn't be so bad, but on my keyboard, it's right next

Re: [R] Stopping ctrl-\ from qutting R

2006-10-06 Thread Martin C. Martin
Thanks! I put it in my .bashrc, now it will never bother me again. Marc Schwartz (via MN) wrote: BTW, you might want to consider updating your FC distro, as FC3 is now EOL Unfortunately, this is my work machine and they're still using FC3. Although I've heard various grumblings about it,

[R] Computing sums of the columns of an array

2005-08-05 Thread Martin C. Martin
Hi, I have a 5x731 array A, and I want to compute the sums of the columns. Currently I do: apply(A, 2, sum) But it turns out, this is slow: 70% of my CPU time is spent here, even though there are many complicated steps in my computation. Is there a faster way? Thanks, Martin

Re: [R] Computing sums of the columns of an array

2005-08-05 Thread Martin C. Martin
Thanks everyone. Gavin Simpson wrote: But neither is that slow on my system. What is A? It's in the middle of a loop. I'm doing some maximum likelihood estimation, and having to evaluate my model over and over again for different parameter values. Thanks, Martin

Re: [R] Statistical significance of a classifier

2005-08-05 Thread Martin C. Martin
Liaw, Andy wrote: From: Martin C. Martin Hi, I have a bunch of data points x from two classes A B, and I'm creating a classifier. So I have a function f(x) which estimates the probability that x is in class A. (I have an equal number of examples of each, so p(class) = 0.5.) One way

[R] Adding sum to derivatives table

2005-08-04 Thread Martin C. Martin
Hi, Trying this: deriv(expression(sum(x)), x) Gives the error message: Function 'sum' is not in the derivatives table I'd like to add it, is this difficult? If not, where is the derivatives table? However, give how basic sum is, I suspect it would have been added if it were

Re: [R] Adding sum to derivatives table

2005-08-04 Thread Martin C. Martin
Rolf Turner wrote: deriv(expression(sum(x)), x) does not make any sense. Good point. But this does: deriv(expression(sum(log(a*x))), a) where a is a scalar. - Martin [[alternative HTML version deleted]] __

Re: [R] Help: how to stop the process when there is a mishandling

2005-08-04 Thread Martin C. Martin
The R for Windows FAQ has many useful tips, including this one: http://cran.r-project.org/bin/windows/base/rw-FAQ.html It's well worth at least reading the contents, so you know what questions it can answer. - Martin Duncan Murdoch wrote: On 8/4/2005 11:39 AM, Shengzhe Wu wrote: and how

[R] Converting a number of minutes to a difftime

2005-05-12 Thread Martin C. Martin
Hi, I have a variable m that contains the number of minutes that something lasted, e.g. m - 139 I'd like to convert it to a difftime, so I can add it to the POSIXct start time. But as.difftime seems to want a character string, with at most two characters for the minutes. All other conversion

[R] Ploting a function of two arguments

2005-03-10 Thread Martin C. Martin
Hi, I've written a function: myfun - function(x, y) { // blah blah } and I want to graph it. My plan is to use persp(), and my question is: how do I create the array of values? One possibility is: x - seq(0, 10, by=.1) y - seq(0, 10, by=.1) inputs - somehow create an array of x,y pairs

Re: [R] trouble loading R function code

2005-03-10 Thread Martin C. Martin
try: source(fhv.R) David Reinke wrote: I created an R function offline with a text editor and saved it as fhv.R When I type in load(file=fhv.R) I get the message Error: bad restore file magic number (file may be corrupted)-- no data loaded I've checked the file, and it opens up fine in Notepad

Re: [R] sample function

2005-03-10 Thread Martin C. Martin
hist is lumping things together. Try: sum(temp == 0) compare to the height of the left most bar. Is this a bug in hist? - Martin mirage sell wrote: Hi everyone, I need help. I want to have a uniform kind distribution. When I used sample function I got almost twice many zeros compared to other