similar to: How to calculate fold change and its confidence interval

Thanks for the proposition. As you see bellow, par("usr") is the same before and after the points() (the full code is bellow): .... > par("usr") [1] -0.2500000 1.2500000 -0.1666667 1.1666667 > # if you remove this points() function, axis will show nothing. > > points(1.5, 1.5, type="p") > p2 <- par(no.readonly=TRUE) > par("usr")

>>>>> Marc Girondot via R-help <r-help at r-project.org> >>>>> on Mon, 24 Jul 2017 09:35:06 +0200 writes: > Thanks for the proposition. As you see bellow, par("usr") is the same > before and after the points() (the full code is bellow): > .... >> par("usr") > [1] -0.2500000 1.2500000 -0.1666667

Thanks... I agree that the problem was explained in the documentation but I can't find a way to have axis() working even manipulating par("plt") or with graphics.reset = TRUE: - adding graphics.reset=TRUE does not allow axis() to be shown; - I see that par()$plt is involved but it is the not sufficient to explain why axis() works because if it is changed by hand, axes are not

Dear Rich, I think that it's generally a bad idea to give statistical (as opposed to simply technical) advice by email without knowing the context of the research. I think that you'd do well to seek help from a statistician, and not just do what I suggest below. Interpolating the data only makes sense if there's no random component to the response (mag in your data). Otherwise, it

Hi marc, Try: par("usr") before and after the call to points and see if it changes. Jim On Sat, Jul 22, 2017 at 12:05 AM, Marc Girondot via R-help <r-help at r-project.org> wrote: > It is known (several discussions on internet) that axis() cannot be used > after fields:::image.plot() (axis() shows nothing). > > However, if points(1.5, 1.5, type="p") is

It is known (several discussions on internet) that axis() cannot be used after fields:::image.plot() (axis() shows nothing). However, if points(1.5, 1.5, type="p") is inserted before the axis() finctions, it works. I have investigated what points(1.5, 1.5, type="p") is doing to allow axis to work and I don't find a solution. par() options are identical (p1 and p2 are

Hi Rich, If I understand your comment about "uniformly distributed along the log=x axis" then I think John's (second) set of commands needs a change to the definition of xx, as follows: xx <- exp(seq(from=log(min(x)),to=log(max(x)),length=50)) m <- lm(y ~ log(x)) yy <- predict(m, data.frame(x=xx)) points(xx, yy, col="red") HTH, Eric On Mon, Sep 25, 2017 at

Hi I believe this is a clipping bug in the 'graphics' package. A simpler version of the problem is this ... plot(1:10) mtext("margin-label", side=2, at=9, las=1, line=1, adj=0) par(mfg=c(1,1)) ## Only the text within the plot region is drawn mtext("margin-label", side=2, at=9, las=1, line=1, adj=0, col="red") I have committed a fix to the development

Dear Rich, Assuming that I understand what you want to do, try adding the following to your script (which, by the way, is more complicated that it needs to be): xx <- 10:50 m <- lm(y ~ x) yy <- predict(m, data.frame(x=xx)) lines(spline(xx, yy), col="blue") m <- lm(y ~ log(x)) yy <- predict(m, data.frame(x=xx)) points(xx, yy, col="magenta") The first set of

Hello everyone I have a couple of questions about the usage of the R function "read.table(.)". My point of departure is that I want to import a matrix (consisting of time and daily stock returns of many stocks) in R. Most of the data is numeric, however some values are missing (blanks) and in other cases I have the character "#DIV/0!" (from excel). My goal is to do some