2016-04-15 32 views
0

我有四個向量,我正在嘗試繪製:此圖的問題。軸未正確

z1 <- rlnorm(100, meanlog=log(6), sdlog=0.3) 
z1 <-ceiling(z1) 
s4 <- seq(0,max(z1),0.1) 
d4 <- dlnorm(s4, meanlog=log(6), sdlog=0.3) 
hist(z1, prob=T, main="", col=rgb(1,0,0,alpha=0.1), ylim=c(0,max(d4)), xlim=c(0, max(z1))) 
par(new=TRUE) 

z2 <- rlnorm(100, meanlog=log(4.2), sdlog=log(2.1)) 
z2 <- ceiling((z2)+15) 
s5 <- seq(0,max(z2),0.1) 
d5 <- dlnorm(s5, meanlog=log(4.2), sdlog=log(2.1)) 
hist(z2, prob=T, main="", col=rgb(1,0,0,alpha=0.15), ylim=c(0,max(d5)), xlim=c(0, max(z2))) 
par(new=TRUE) 

m1 <- rlnorm(100, meanlog=2, sdlog=0.5) 
m1 <- ceiling((m1)+5) 
s1 <- seq(0,max(m1),0.1) 
d1 <- dlnorm(s1, meanlog=log(2), sdlog=0.5) 
hist(m1, prob=T, main="", col=rgb(0,1,0,alpha=0.1) , ylim=c(0,max(d1)), xlim=c(0, max(m1))) 
par(new=TRUE) 

m2 <- rnorm(100, mean=16, sd=5) 
m2 <- ceiling(m2) 
s2 <- seq(0,max(m2),0.1) 
d2 <- dnorm(s2, mean=16, sd=5) 
hist(m2, prob=T, main="", col=rgb(0,1,0,alpha=0.15) , ylim=c(0,max(d2)), xlim=c(0, max(m2))) 
par(new=TRUE) 

我有密謀這個graph問題。 軸不正確。

謝謝你的幫助。

+0

你想達到什麼目的?你的代碼的第一部分會拋出一個錯誤,因爲'm1'沒有被定義。你可能想要清理你的代碼。 – epo3

+0

我會比較這些組,但該圖的比例是錯誤的。 x軸現在是:5,10,10,20 .....並且不能讀取y軸。 – Sisi

回答

0

使用的add=TRUE參數,可以將所有直方圖繪製在同一軸上。

hist(z1, prob=T, main="", col=rgb(1,0,0,alpha=0.1), ylim=c(0,max(d4)), xlim=c(0, max(c(z1,z2, m1, m2)))) 
hist(z2, prob=T, main="", col=rgb(1,0,0,alpha=0.15), add=T) 
hist(m1, prob=T, main="", col=rgb(0,1,0,alpha=0.1), add=T) 
hist(m2, prob=T, main="", col=rgb(0,1,0,alpha=0.15), add=T)