我正在处理大量数据,并使用julia中的glm计算了逻辑回归。我所做的是检查找工作率的不同影响,我得到0或1。 现在我想绘制我的结果。我有吉尔曼/山丘书,他们做了类似的练习,但在R(也用glm)。 他们的绘图是:
fit.1 <- glm(switch ~ dist, family=binomial(link="logit"))
jitter.binary <- (a, jitt=.05){
ifelse(a==0, runif(length(a), 0, jitt), runif(length(a), 1-jitt, 1))
}
switch.jitter <- jitter.binary(X)
plot(dist, switch.jitter)
curve(invlogit (coef(fit.1)[1] + coef(fit.1)[2]*x), add=TRUE)
我的问题很基本我认为,我怎么能在朱莉娅做类似的事情,因为我不太喜欢R代码。 对于情节,我使用Plots.jl
我希望有人能够帮助我,如果您需要更多信息,请告诉我。
干杯
答案 0 :(得分:1)
using Plots
plot(x, predict(fit1), seriestype = :line) #where x is your predictor variable
更新以包含@pkofod的评论