[程式] R在plot時候想限定y軸座標

看板Statistics作者 (anio)時間15年前 (2010/05/17 12:22), 編輯推噓0(004)
留言4則, 2人參與, 最新討論串1/1
[軟體程式類別]: 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

05/17 13:48, , 1F
程式排版一下吧
05/17 13:48, 1F

05/17 13:49, , 2F
你需要把迴圈中的結果存在陣列中
05/17 13:49, 2F

05/17 13:50, , 3F
迴圈結束後,再用range()取得範圍,然後再一次畫出來.
05/17 13:50, 3F
不是很會排版@@ 迴圈中的結果應該是存在V1跟E中 我用V2跟p做描點 C大 我沒弄懂你說的range要怎麼用? ※ 編輯: ANIO1021 來自: 114.38.49.196 (05/17 14:25)

05/17 18:31, , 4F
最後只好用這招ylim=c(0,0.02)勉強算可以了
05/17 18:31, 4F
文章代碼(AID): #1ByCHp_L (Statistics)