0

ロジスティック回帰モデルの 20 個の切片と係数のセットがあります。それらすべてを1つのグラフにプロットしたいと思います（理想的には、ベースであり、ggplotではありません）。次のような単一の曲線をプロットします。

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

線形モデルでこれを行うには、次を使用できますmapply。

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

ループを使用したり、ループの使用を回避しmapplyたりする方法はありますか?curveplogis

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 に答える