绘图时par(）的作用置疑

在做图时，无意中发现了这样一个问题：

1. op<-par(bg="light blue")
   
2. plot(c(-2, 2), c(-2, 2))
   
3. rect(-2, -2, 1, 1, col="cornsilk")
   
4. par(op)

复制代码

发现，最后一行（par(op)无论加与不加，结果都是一样的，都出现下图：

怪哉怪哉！

这个不怪,不管加不加最后一行都会出现你所给的四幅图,但是你没理解par(op)
的作用,他是把当前图形设备恢复到op <- par(mfrow=c(2,2), mar=c(4,4,3,2))
之前的设置.

为了方便你理解,你可以运行一下.

1. library(Zelig)d <- data.frame(  
   
2.   y = crabs$sp,  
   
3.   x1 = crabs$FL,  
   
4.   x2 = crabs$RW)r <- zelig( y ~ x1 + x2, model="probit", data=d )op <- par(mfrow=c(2,2), mar=c(4,4,3,2))plot(r) plot(1:10)
   

复制代码

就能明白,先运行plot(r)看一下图形
然后再运行plot(1:10)

kaifengedu 发表于 2013-10-4 10:39
这个不怪,不管加不加最后一行都会出现你所给的四幅图,但是你没理解par(op)
的作用,他是把当前图形设备恢复 ...

我好像明白了，第一个op<-par()，是设定图形参数，而第二个par(op)表面意思是再设一次，实际功能是恢复到原来状态，是这样么？老师？
尝试了一下，好像我理解的还是错误，见下面运行：
par(bg="light gray")
plot(c(-2,  2), c(-2,  2))  
rect(-2,  -2,  1,  1,  col="cornsilk")
par(bg="light gray")
运行之后，发现还是没有恢复。我的疑惑是：非得把par()赋值给一个名称，如op(我尝试了，用o也可以），才可以吗？死记当然可以，就是有些不理解。
非常感谢您！

耕耘使者 发表于 2013-10-4 20:45
我好像明白了，第一个op

你的没恢复是什么意思?
你看你的两条命令都是
par(bg="light gray")
"我好像明白了，第一个op<-par()，是设定图形参数，而第二个par(op)表面意思是再设一次，实际功能是恢复到原来状态，是这样么？老师？"
这个理解是对的.
"我的疑惑是：非得把par()赋值给一个名称，如op(我尝试了，用o也可以），才可以吗？死记当然可以，就是有些不理解。"
其实把par() 赋值给一个名称.个人感觉也不太需要怎么记啊.
就相当于开始i=0,你中间赋值了i=1,想用i=0时间再恢复一下即可.

天哪，楼主竟然有这么多论坛币

kaifengedu 发表于 2013-10-4 21:14
你的没恢复是什么意思?
你看你的两条命令都是
par(bg="light gray")

副教授，我理解的还是不对，看下面代码，恢复不了了！！
op<-par(bg="light blue")
plot(c(-2,  2), c(-2,  2))  
rect(-2,  -2,  1,  1,  col="cornsilk")
par(op)
plot(1:10)

tmdxyz 发表于 2013-10-5 05:54
天哪，楼主竟然有这么多论坛币

这是原始资本积累。从入坛太久了，而开始时，赚取论坛币比较容易。
当然，容易也不相对的，我上传了大量的数据年鉴。

耕耘使者 发表于 2013-10-5 08:26
副教授，我理解的还是不对，看下面代码，恢复不了了！！
op

我给你个截图你看看吧,我不知道你哪里没理解,呵呵.
运行了

1. op<-par(bg="light blue")
   
2. plot(c(-2,  2), c(-2,  2))  
   
3. rect(-2,  -2,  1,  1,  col="cornsilk")

复制代码

生成了图像


运行

1. par(op)
   
2. plot(1:10)

复制代码

生成了



BTW:以后如果方便请把代码放在代码框里方便别人复制粘贴,呵呵.

kaifengedu 发表于 2013-10-5 09:02
我给你个截图你看看吧,我不知道你哪里没理解,呵呵.
运行了

又尝试了一下，没有问题了。上次软件出了bug。R就是这点不好，有时莫名其妙的出bug，关闭重新做又好了。
多谢！
顺便请教，怎样把代码放在代码框中？

耕耘使者 发表于 2013-10-5 12:39
又尝试了一下，没有问题了。上次软件出了bug。R就是这点不好，有时莫名其妙的出bug，关闭重新做又好了。
...

你可以用回复里笑脸旁边的<>,一点就可以了.