使用mapply绘制多条逻辑曲线

Question

对于logistic回归模型，我有一组20个截距和系数。我想把它们都画在一个图上(理想情况下，使用base而不是ggplot)。我绘制一条曲线，如下所示：

x = datapoints
curve(plogis(intercept + coef*x), min(x), max(x), add=TRUE) 

要对线性模型执行此操作，我可以使用mapply

a = intercepts
b = coefs
mapply(abline, a,b)

有没有一种方法可以使用mapply或类似的curve和plogis来避免使用循环？

Answer 1

您没有提供可重现的示例，但这是一个概念证明：

x = runif(20)/2
plot(x, bty="l", ylim=c(0,1))

a = runif(20)/100
b = runif(20)/4
my.func <-function(intercept,coefs) curve(plogis(intercept+coefs*x), add=TRUE)
mapply(my.func,a,b)

​

​