[Rd] discrepancy between lm and MASS:rlm

2011-03-14 Thread Vadim Ogranovich
Dear R-devel, There seems to be a discrepancy in the order in which lm and rlm evaluate their arguments. This causes rlm to sometimes produce an error where lm is just fine. Here is a little script that illustrate the issue: library(MASS) ## create data n - 100 dat -

Re: [Rd] discrepancy between lm and MASS:rlm

2011-03-14 Thread Vadim Ogranovich
Dunlap [mailto:wdun...@tibco.com] Sent: Monday, March 14, 2011 1:39 PM To: Vadim Ogranovich; r-devel@r-project.org Subject: RE: [Rd] discrepancy between lm and MASS:rlm -Original Message- From: r-devel-boun...@r-project.org [mailto:r-devel-boun...@r-project.org] On Behalf Of Vadim

[Rd] calling browser on error

2010-10-15 Thread Vadim Ogranovich
Dear R-developers, I am trying to figure out a way to call browser() when an error occur, and naturally I want the browser() to be called in the environment of the error. I tried something simple in vain: f - function() { x - 1; stop('ok') } tryCatch(f(), error=browser()) Called from:

Re: [Rd] calling browser on error

2010-10-15 Thread Vadim Ogranovich
48839 language R version.string R version 2.9.1 (2009-06-26) -Original Message- From: Joshua Ulrich [mailto:josh.m.ulr...@gmail.com] Sent: Friday, October 15, 2010 1:31 PM To: Vadim Ogranovich Cc: r-devel@r-project.org Subject: Re: [Rd] calling browser on error I believe options

Re: [Rd] [R] converting result of substitute to 'ordidnary' expression

2010-06-26 Thread Vadim Ogranovich
...@tajo.ucsd.edu] Sent: Friday, June 25, 2010 9:16 PM To: Vadim Ogranovich Cc: 'r-h...@r-project.org' Subject: Re: [R] converting result of substitute to 'ordidnary' expression On Fri, 25 Jun 2010, Vadim Ogranovich wrote: Dear R users, As substitute() help page points out: Substituting

[Rd] Save all objects in all environment stack

2009-09-14 Thread Vadim Ogranovich
Dear R Developers, I am writing a function that would save all objects in all environments 'visible' from the point where the function is called (very much like save.image, but kind of saving the entire stack). I figured that the call lapply (seq(0, sys.nframe()), function(i) ls(sys.frame(i),

Re: [Rd] bug in approx crashes R

2009-07-28 Thread Vadim Ogranovich
Thank you Martin! The new 'rule'-s is exactly what I wanted. It's intuitive and consistent with the scalar use. Thanks, Vadim -Original Message- From: Martin Maechler [mailto:maech...@stat.math.ethz.ch] Sent: Tuesday, July 28, 2009 5:00 AM To: Vadim Ogranovich Cc: 'William Dunlap'; r

Re: [Rd] bug in approx crashes R

2009-07-27 Thread Vadim Ogranovich
message when yleft, yright, and f are set to non-scalars is better than silently returning NA. Thanks, Vadim -Original Message- From: William Dunlap [mailto:wdun...@tibco.com] Sent: Monday, July 27, 2009 12:14 PM To: Vadim Ogranovich; r-devel@r-project.org Subject: RE: [Rd] bug in approx

[Rd] bug in approx crashes R

2009-07-21 Thread Vadim Ogranovich
Dear R-devel, The following line crashes R approx(1, 1, 0, method='const', rule=2, f=0, yleft=NULL, ties='ordered')$y Process R:2 exited abnormally with code 5 at Tue Jul 21 14:18:09 2009 version _ platform i386-pc-mingw32 arch i386 os mingw32

Re: [Rd] lsfit w/ rank-deficient x

2009-03-17 Thread Vadim Ogranovich
Actually, the correct permutation is given by the inverse of qr$pivot: foo$coefficients[foo$qr$pivot] - foo$coefficients Here foo is an object returned by lsfit, see below. -Original Message- From: Vadim Ogranovich Sent: Friday, March 13, 2009 5:25 PM To: 'r-devel@r-project.org

[Rd] lsfit w/ rank-deficient x

2009-03-13 Thread Vadim Ogranovich
Dear R-devel, It seems that lsfit incorrectly reports coefficients when the input matrix 'x' is rank-deficient, see the example below: ## here values of 'b' and 'c' are incorrectly swapped x - cbind(a=rnorm(100), b=0, c=rnorm(100)); y - rnorm(100); lsfit(x, y)$coef Intercept a

[Rd] weekly update 20081208

2008-12-08 Thread Vadim Ogranovich
In Progress: * EURO DOLLAR data for trading Z* and vice versa * scripts for computing asset correlations w/ application to ED * improving R code for collecting stats accross dates and groups of data. Completed: * New signals for US Equity (based on PK alpha). Brandon will simulate and deploy if

[Rd] Rscript thank you

2007-12-19 Thread Vadim Ogranovich
Hi, I'd like to say thank you to the developers of Rscript. It is a nice tool that allows integration of R into multi-tool work flow and, in my experience so far, it does exactly what is expected of a utility program. Thank you very much! Vadim [[alternative HTML version