2013-10-02 16 views
0

我正在製作一個很好的圖表,但我遇到了一個問題,我畫的線條的末端。我不知道爲什麼R直到maxxx(這裏爲10)都沒有繪製直線。您可以從下面什麼我談論(我不能把圖片設置爲我的名聲太低)碼圖中看到,該線的末端10之前停止:R plot - lines()在結束之前停止的小問題值

maxxx<-10 
plot(0,type="n",axes=FALSE,xlim=c(0,10),ylim=c(0,10),ylab="",xlab="") 
mtext("TROLOLOL",side=3,cex=3) 
axis(1,pos=0,at=c(0,10),labels=FALSE) 
mtext("R1 ",side=1,line=0,cex=3) 
axis(2,pos=0,at=c(0,10),labels=FALSE) 
mtext("R2",side=2,line=0,cex=3) 


Given_growth_rates<-c(0,0.3,0.5,1,2,5); 


K1<-7 
g1<-5.5 
m<-0.2 
R1_isoclines<-numeric(length(Given_growth_rates)) 
for (i in 1:length(Given_growth_rates)){ 
R1_isoclines[i]<-((Given_growth_rates[i]+m)*K1)/(g1-Given_growth_rates[i]-m) 
} 
R1_isoclines 

K2<-10 
g2<-7 
R2_isoclines<-numeric(length(Given_growth_rates)) 
for (i in 1:length(Given_growth_rates)){ 
R2_isoclines[i]<-((Given_growth_rates[i]+m)*K2)/(g2-Given_growth_rates[i]-m) 
} 
R2_isoclines 

for (i in 1:length(R1_isoclines)){ 
lines(rep(R1_isoclines[i],times=length(R2_isoclines[i]:maxxx)), R2_isoclines[i]:maxxx,    col=i+1, type="l") 
lines(R1_isoclines[i]:maxxx,rep(R2_isoclines[i],times=length(R1_isoclines[i]:maxxx)),col=i+1,type="l") 
} 

回答

1

嘗試更換四強通過您的代碼行:

for (i in 1:length(R1_isoclines)){ 
    lines(rep(R1_isoclines[i],times=2), c(R2_isoclines[i],maxxx), col=i+1, type="l") 
    lines(c(R1_isoclines[i], maxxx), rep(R2_isoclines[i],times=2),col=i+1,type="l") 
} 

enter image description here

也許已經做了你所期望的。

+0

是的,它的工作,謝謝:D但爲什麼我的代碼沒有工作? – Ouistiti

+0

由於「:」 - 運算符的功能,例如檢查0.5:10的結果,其最後的值將是9.5 – Sophia

+0

非常感謝!我會更仔細地使用這個操作符。祝你有美好的一天:D – Ouistiti

相關問題