Re: [Rd] ifelse() woes ... can we agree on a ifelse2() ?

2016-08-07 Thread Uwe Ligges
On 06.08.2016 17:30, Duncan Murdoch wrote: On 06/08/2016 10:18 AM, Martin Maechler wrote: Dear R-devel readers, ( = people interested in the improvement and development of R). This is not the first time that this topic is raised. and I am in now state to promise that anything will result from

[Rd] table(exclude = NULL) always includes NA

2016-08-07 Thread Suharto Anggono Suharto Anggono via R-devel
This is an example from https://stat.ethz.ch/pipermail/r-help/2007-May/132573.html . With R 2.7.2: > a <- c(1, 1, 2, 2, NA, 3); b <- c(2, 1, 1, 1, 1, 1) > table(a, b, exclude = NULL) b a 1 2 11 1 22 0 31 0 1 0 With R 3.3.1: > a <- c(1, 1, 2, 2, NA, 3); b <- c(2,

[Rd] problem with abine(lm(...)) for plot(y~x, log='xy')

2016-08-07 Thread Spencer Graves
Hello: In the following plot, the fitted line plots 100 percent above the points: tstDat <- data.frame(x=10^(1:3), y=10^(1:3+.1*rnorm(3))) tstFit <- lm(log(y)~log(x), tstDat) plot(y~x, tstDat, log='xy') abline(tstFit) I can get the correct line with the following: tstPredDa

Re: [Rd] problem with abine(lm(...)) for plot(y~x, log='xy')

2016-08-07 Thread peter dalgaard
Try log10()... -pd > On 07 Aug 2016, at 21:03 , Spencer Graves wrote: > > Hello: > > > In the following plot, the fitted line plots 100 percent above the > points: > > > tstDat <- data.frame(x=10^(1:3), y=10^(1:3+.1*rnorm(3))) > tstFit <- lm(log(y)~log(x), tstDat) > plot(y~x, tstDat,