作者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