为什么不在一张图中绘制多个函数？

Question

我在尝试根据逆高斯分布在一张图中绘制多个可靠性函数时遇到一些问题。我需要函数是线，而我得到的只是点，当尝试设置 type="l" 时，它恰好到处都是乱七八糟的画线。

这是代码

library("statmod")
x<-rinvgauss(90,0.000471176,0.0000191925)
y<-rinvgauss(90,0.000732085,0.000002982015)
z<-rinvgauss(180,0.000286672,0.00000116771)

den<-pinvgauss(x,0.000471176,0.0000191925)
dens<-pinvgauss(y,0.000732085,0.000002982015)
densi<-pinvgauss(z,0.000286672,0.00000116771)

rel<-1-den
reli<-1-dens
relia<-1-densi

plot(x,rel, xlim=c(0,0.002), col="red")
points(y,reli, col="blue")
points(z,relia, col="black")

如果能提供任何帮助，我将不胜感激！

Answer 1

问题是您的 x、y、z 值未排序...

library("statmod")
x <- sort(rinvgauss(90,0.000471176,0.0000191925))
y <- sort(rinvgauss(90,0.000732085,0.000002982015))
z <- sort(rinvgauss(180,0.000286672,0.00000116771))

den   <- pinvgauss(x,0.000471176,0.0000191925)
dens  <- pinvgauss(y,0.000732085,0.000002982015)
densi <- pinvgauss(z,0.000286672,0.00000116771)

rel   <- 1-den
reli  <- 1-dens
relia <- 1-densi

plot(x,rel, xlim=c(0,0.002), col="red", type="l")
lines(y,reli, col="blue")
lines(z,relia, col="black")

Answer 2

您的值未排序。这应该有效：

    x<-sort(rinvgauss(90,0.000471176,0.0000191925))
    y<-sort(rinvgauss(90,0.000732085,0.000002982015))
    z<-sort(rinvgauss(180,0.000286672,0.00000116771))

den<-sort(pinvgauss(x,0.000471176,0.0000191925))
dens<-sort(pinvgauss(y,0.000732085,0.000002982015))
densi<-sort(pinvgauss(z,0.000286672,0.00000116771))

rel<-1-den
reli<-1-dens
relia<-1-densi

plot(x,rel, xlim=c(0,0.002), col="red",type="l")
lines(y,reli, col="blue")
lines(z,relia, col="black")

为什么不在一张图中绘制多个函数？

Why isn't this plotting multiple functions in one graph?

r

plot

reliability