我试图仅在一个图中绘制图形 . 我使用一个函数有4个不同的图 . 这是我的代码:
hazard.plot.w2p(beta = beta.spreda, eta = eta.spreda, time = exa1.dat$time, line.colour = "blue")
hazard.plot.w2p(beta = 1.076429, eta = 26.21113, time = exa1.dat$time, line.colour = "blue")
hazard.plot.w2p(beta = 5, eta = 32.97954, time = exa1.dat$time, line.colour = "blue")
hazard.plot.w2p(beta = 2, eta = 32.9795, time = exa1.dat$time, line.colour = "blue")
这是我用来获取输出的函数:
hazard.plot.w2p <- function(beta, eta, time, line.colour, nincr = 500) {
max.time <- max(time, na.rm = F)
t <- seq(0, max.time, length.out = nincr)
r <- numeric(length(t))
for (i in 1:length(t)) {
r[i] <- failure.rate.w2p(beta, eta, t[i])
}
plot(t, r, type = 'l', bty = 'l',
col = line.colour, lwd = 2,
main = "", xlab = "Time",
ylab = "Failure rate",
las = 1, adj = 0.5,
cex.axis = 0.85, cex.lab = 1.2)
}
我想仅在一个图中绘制所有4个图 .
这是一个示例数据集:
fail time
a 4.55
a 4.65
a 5.21
b 3.21
a 1.21
a 5.65
a 7.12
1 回答
您可以在每个绘图后使用
par(new=TRUE)
. 或者,您可以将它放在函数的末尾 .