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.