ggplotで項目特性曲線を描く

IRT（項目反応理論）の項目特性曲線をRのggplot2で描画するスクリプトです。2パラメータ・ロジスティックモデルを想定しています。

library(ggplot2)

a=1;b=0
p=function(theta) 1/(1+exp(-a*(theta-b)))

g=ggplot(data=data.frame(X=c(-4,4)), aes(x=X))
 +stat_function(fun=p,size=1.5,color="orange")
 +ggtitle("ICC")+labs(x=expression(theta), y="prob")
 +theme(text=element_text(size=15))
plot(g)

下のようなグラフが描画されます。