我有一组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
类似的内容来避免使用循环?
答案 0 :(得分: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)