0

我有一组 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来避免使用循环?

4

1 回答 1

3

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

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)

在此处输入图像描述

于 2013-05-21T00:44:50.283 回答