作者ANIO1021 (anio)
看板Statistics
標題[程式] R在plot時候想限定y軸座標
時間Mon May 17 12:22:09 2010
[軟體程式類別]:
r
[程式問題]:
繪圖
[軟體熟悉度]:
低(1~3個月)
[問題敘述]:
我要plot出一個圖
控制圖的y軸在（0,0.02)區間是0.002
試了很多次plot(x,y,ylim=a+range(b))
a,b嘗試了很多組合 但仍舊沒辦法達到我的要求
煩請版上高手教我怎麼設定
[程式範例]:
E=matrix(0,100,25)
V1=matrix(0,100,25)
V2=matrix(0,100,1)
G=matrix(0,100,1)
N=25
k=1
w=1:100
p=0.01*w
for(i in 1:100)
{
for(j in 1:25)
{
E[i,j]=((j/N)^(1/k))*(factorial(N)/(factorial(j)*factorial(N-j)))*((1-0.01*i)^(k*j))*((1-(1-0.01*i)^k)^(N-j))
}
G=1-apply(E,1,sum)
}
for(a in 1:100)
{
for(b in 1:25)
{
V1[a,b]=((b/N)^(2/k))*(factorial(N)/(factorial(b)*factorial(N-b)))*((1-0.01*a)^(k*b))*((1-(1-0.01*a)^k)^(N-b))
}
V2=apply(V1,1,sum)-(1-G)^2
}
plot(p,V2,xlab=expression("true value of p"),ylab=expression("variance of
estimator p"),type="l",ylim=0.01+range(0.005));
我主要要把k值用1 2 3 5 7 10 15 25代入並且畫在同張圖上
但就是y軸部份我弄不到我想要得區間
請高手幫我
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 114.38.49.196
1F：→ clickhere:程式排版一下吧 05/17 13:48
2F：→ clickhere:你需要把迴圈中的結果存在陣列中 05/17 13:49
3F：→ clickhere:迴圈結束後,再用range()取得範圍,然後再一次畫出來. 05/17 13:50
不是很會排版＠＠
迴圈中的結果應該是存在V1跟E中
我用V2跟p做描點
C大 我沒弄懂你說的range要怎麼用？
※ 編輯: ANIO1021 來自: 114.38.49.196 (05/17 14:25)
4F：→ ANIO1021:最後只好用這招ylim=c(0,0.02)勉強算可以了 05/17 18:31