1
我试图绘制对同一地块不同的方程一些累积分布函数,这是我写的代码:情节没有显示CDF
library (triangle)
library(lattice)
library(latticeExtra)
n = 1000
W1 = rtriangle(n,3128,3250)
W2 = rtriangle(n,3065,3149)
SO = rtriangle(n,0.2,0.3)
MCtab <- data.frame(W1,W2,SO)
set.seed(1)
for (n in 1:n) {
NPV30 <- (1*W1 + 2*W2 + 3*SO)}
set.seed(1)
for (n in 1:n) {
NPV50 <- ecdf((4*W1 + 5*W2 + 6*SO))}
set.seed(1)
for (n in 1:n) {
NPV70 <- ecdf((7*W1 + 8*W2 + 9*SO))}
MCtab2 <- data.frame(NPV30,NPV50,NPV70)
plot(NPV30, verticals=TRUE, main= 'Polymer NPV CDF', do.points=FALSE, col='red')
plot(NPV50, verticals=TRUE, do.points=FALSE, add=TRUE, col='brown')
plot(NPV70, verticals=TRUE, do.points=FALSE, add=TRUE, col='orange')
只有一个曲线如图不能似乎找出原因,或者如果有人有更好的方法会很好。谢谢
一种可能性是,后来的值太大初始待观察情节维度。您可以添加包含所有数据点的最小值和最大值的ylim参数。 – lmo
@lmo谢谢,这似乎是最有可能的原因 – John
对我的答案的任何反馈? –