Clive Nicholas
2014-Jun-01 20:20 UTC
[R] Plotting issues with -par(mfrow=c(r,c))- after -boot-
R 3.1.0 / RStudio 0.98.507 / OpenSUSE Linux 13.1 I'm having difficulties getting -par(mfrow=c(r,c))- to work as it should after running bootstrapped regression models using -boot-. An example (tested): test=data.frame(A=rnorm(500, mean=2.72, sd=5.36), B=sample(c(12,20,24,28,32),size=500,prob=c(0.333,0.026,0.026,0.436,0.179),replace=TRUE), C=sample(c(0,1),size=500,prob=c(0.5,0.5),replace=TRUE),D=sample(c(0,1),size=500,prob=c(0.564,0.436),replace=TRUE)) library(boot) bs=function(formula, data, indices) { test=data[indices,] fit=lm(formula, data=test) return(coef(fit)) } results=boot(data=test, statistic=bs, R=10000, formula=A~B+C+D+C*D) results When I run par(mfrow=c(2,2)) and then plot(results, index=2) plot(results, index=3) plot(results, index=4) plot(results, index=5) the plots display in full size only, as they are called, and not as one combined image of four plots. This problem also happened to me last year under earlier versions of R in Kubuntu Linux. This command worked fine when I first started using R in Linux, I don't think I'm doing much wrong in this instance and I would welcome any explanation as to what exactly is going on here, as well as a solution. -- Clive Nicholas "My colleagues in the social sciences talk a great deal about methodology. I prefer to call it style." -- Freeman J. Dyson [[alternative HTML version deleted]]
Jeff Newmiller
2014-Jun-01 21:07 UTC
[R] Plotting issues with -par(mfrow=c(r,c))- after -boot-
Read ?plot.boot, in particular the Side Effects section. Not a very friendly function... you will have to modify it if you wish to proceed. --------------------------------------------------------------------------- Jeff Newmiller The ..... ..... Go Live... DCN:<jdnewmil at dcn.davis.ca.us> Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research Engineer (Solar/Batteries O.O#. #.O#. with /Software/Embedded Controllers) .OO#. .OO#. rocks...1k --------------------------------------------------------------------------- Sent from my phone. Please excuse my brevity. On June 1, 2014 1:20:45 PM PDT, Clive Nicholas <clivelists at googlemail.com> wrote:>R 3.1.0 / RStudio 0.98.507 / OpenSUSE Linux 13.1 > >I'm having difficulties getting -par(mfrow=c(r,c))- to work as it >should >after running bootstrapped regression models using -boot-. An example >(tested): > >test=data.frame(A=rnorm(500, mean=2.72, sd=5.36), >B=sample(c(12,20,24,28,32),size=500,prob=c(0.333,0.026,0.026,0.436,0.179),replace=TRUE), >C=sample(c(0,1),size=500,prob=c(0.5,0.5),replace=TRUE),D=sample(c(0,1),size=500,prob=c(0.564,0.436),replace=TRUE)) >library(boot) >bs=function(formula, data, indices) { > test=data[indices,] > fit=lm(formula, data=test) > return(coef(fit)) >} >results=boot(data=test, statistic=bs, R=10000, formula=A~B+C+D+C*D) >results > >When I run > >par(mfrow=c(2,2)) > >and then > >plot(results, index=2) >plot(results, index=3) >plot(results, index=4) >plot(results, index=5) > >the plots display in full size only, as they are called, and not as one >combined image of four plots. This problem also happened to me last >year >under earlier versions of R in Kubuntu Linux. > >This command worked fine when I first started using R in Linux, I don't >think I'm doing much wrong in this instance and I would welcome any >explanation as to what exactly is going on here, as well as a solution.