[R] sas vs. R

2004-11-21 Thread ivo_welch-rstat8783
SAS * better manuals. * tech support for most universities contracted into the price, thus for researchers. * batch orientation. if you have to handle data sets that are as large as your memory, SAS generally does it better. It seems to be an n-pass design. Years ago, when memory was

[R] non-linear solve?

2004-11-01 Thread ivo_welch-rstat8783
hi: could someone please point me to a function that allows me to solve general non-linear functions? irr.in - function(r, c1, c2, c3 ) { return(c1+c2/(1+r) + c3/(1+r)^2); } solve.nonlinear( irr.in, -100, 60, 70 ); 0.189 If someone has written an irr function, this would be helpful,

[R] pdf device --- axes labels text disappeared?

2004-10-15 Thread ivo_welch-rstat8783
Dear R Wizards: Running R 1.9.1. on amd64. Promise- c(0,20,40); Expect- c(0, 20, 0.2*20+.8*40 ); # this omits printing numbers on the axes labels. pdf(file = bug.pdf ); plot(Promise, Expect, type=b, ylim=c(0,60)); dev.off(); # this works postscript(file = bug.eps ); plot(Promise, Expect, type=b,

Re: [R] pdf device --- axes labels text disappeared?

2004-10-15 Thread ivo_welch-rstat8783
Hi: Thank you everyone. You were all correct. This turns out to be a viewer bug---and even more likely a viewer installation bug, not an xpdf bug, much less an R bug. Apologies for having wasted everyone's time. Regards, /ivo --- ivo welch __

[R] adding observations to lm for fast recursive residuals?

2004-09-15 Thread ivo_welch-rstat8783
dear R community: i have been looking but failed to find the following: is there a function in R that updates a plain OLS lm() model with one additional observation, so that I can write a function that computes recursive residuals *quickly*? PS: (I looked at package strucchange, but if I am

[R] loadhistory() in .Rprofile ?

2004-08-21 Thread ivo_welch-rstat8783
dear wizards: my .Rprofile has just one command for testing, loadhistory(~/.Rhistory) but this gives me an error on R startup: Error: couldn't find function loadhistory Invoking loadhistory() as the first interactive command works fine; incidentally, I believe loadhistory() in the .Rprofile

[R] R Cookbook

2004-08-14 Thread ivo_welch-rstat8783
is anyone writing an R cookbook (ala the perl cookbook)? this would be more for programming and graphics task than a statistics textbook. This seems more like a manufacturing defect than a random occurrance. if not, if I can fit it into my schedule, I may start one slowly on my website

[R] expression + paste + arguments + ...

2004-07-29 Thread ivo_welch-rstat8783
dear R wizards: I would like to write a function that roughly places the equivalent of the following latex text into the current plot: \newcommand{ \placesigma }[4]{ \put(\#1,\#2){ \sigma_{A , #3} = #4 } I cannot figure out how to do this. I know I have to use a function that uses

[R] omit complete cases

2004-07-07 Thread ivo_welch-rstat8783
thanks everyone. all solutions were better than what I had, and simple. R is an interesting experience. Extremely powerful and awe-inspiring for its elegance; things work like I would never have believed how elegantly they work. The IQ in the subsetting alone is superbly clever. And

[R] options --- na.rm , save-on-exit

2004-07-03 Thread ivo_welch-rstat8783
dear R wizzards: two options() related questions. [a] options(na.action): many of my variables have missing variables, so I would like to set a default to ignore, especially in univariate functions like mean() and median() . alas, without na.rm=T in the functions themselves, I always get an