plotCI : 두 변수의 플롯을 오버레이하는 방법

Question

나는 신뢰 구간으로 포식자와 먹이의 집단을 시간의 경과에 따라 플롯하려고합니다. 나는이 두 그래프를 어떻게 구별 할 수 있습니까? 동일한 그래프에 그래프를 그리는 방법은 무엇입니까?

#take mean, number, and create se of prey(d) 
d.means=tapply(mydata$prey,mydata$week, mean) 
d.n=tapply(mydata$prey,mydata$week, length) 
d.se=tapply(mydata$prey,mydata$week, sd)/sqrt(d.n) 
#plot with se using plotrix 
plotCI(as.numeric(row.names(d.means)),d.means,d.se,ylim=c(0,400),pch=19,gap=0,xlab="Week",ylab="d, w population") 




#take mean, number, and create se of predator(w) 
w.means=tapply(mydata$pred,mydata$week, mean) 
w.n=tapply(mydata$pred,mydata$week, length) 
w.se=tapply(mydata$pred,mydata$week, sd)/sqrt(w.n) 
#plot with se using plotrix 
plotCI(as.numeric(row.names(w.means)),w.means,w.se,ylim=c(0,400),pch=19,gap=0,xlab="Week",ylab="d, w population") 

Answer 1

첫 번째 플롯 후, 다음 플롯을 플롯하기 전에 아래의 코드를 사용

par(new=T) 

당신이 두 플롯을 수용하기 xlim 및 ylim을 설정할 수 있는지 확인합니다. 그리고 axes = F 및 ann = F 옵션을 사용해야합니다.

이러한 그래픽 기능은 전자 책 "R Fundamentals & 그래픽"에서 자세히 설명합니다. 데스크 레퍼런스로 사용할 수도 있습니다.

Answer 2

#take mean, number, and create se of prey(d) 
d.means=tapply(mydata$prey,mydata$week, mean) 
d.n=tapply(mydata$prey,mydata$week, length) 
d.se=tapply(mydata$prey,mydata$week, sd)/sqrt(d.n) 


#take mean, number, and create se of predator(w) 
w.means=tapply(mydata$pred,mydata$week, mean) 
w.n=tapply(mydata$pred,mydata$week, length) 
w.se=tapply(mydata$pred,mydata$week, sd)/sqrt(w.n) 

여기에 당신은 당신이 필요로하는 모든 변수를 만든하지만 당신은 그들이 육식 동물 또는 먹이 어떤가를 나타내는 변수로 키가 큰 데이터 세트에있을 그들을 필요 ggplot 사용하여 플롯 할 수 있습니다. 나는 또한 시간 변수를 추가했다. 나는 네가 일주일이라고 생각한다. 나는 당신의 데이터가 그래서 여기처럼 정확히 모르는

x=data.frame(means=c(w.means,d.means), 
      n=c(w.n,d.n), 
      se=c(w.se,d.se), 
      role=c(rep("pred",length(w.n)),rep("prey",length(d.n))), 
      time=c(1:length(w.n),1:length(d.n)) 
      ) 

난 그냥 형식을 설명하기 위해 최대 요리 가짜입니다.

 means n   se role time 
1 0.9874234 10 0.16200575 pred 1 
2 1.4120207 12 0.08895026 pred 2 
3 2.7352516 8 0.07991036 pred 3 
4 1.1301248 11 0.05481813 prey 1 
5 2.4810040 13 0.28682585 prey 2 
6 3.1546947 9 0.22126054 prey 3 

일단 데이터가 ggplot을 사용하여이 멋진 형식으로되어 있으면 매우 쉽습니다. 이 제공

ggplot(x, aes(x=time, y=means, colour=role)) + 
    geom_errorbar(aes(ymin=means-se, ymax=means+se), width=.1) + 
    geom_line() 

: