我正在处理大量数据,刚刚在 julia 中使用 glm 计算了逻辑回归。我所做的是检查对求职率的不同影响,我有 0 或 1。现在我想绘制我的结果。我有 gilman/hill 书,他们在其中进行了类似的练习,但在 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)
我认为我的问题很基本,我怎么能在 julia 中做类似的事情,因为我不太喜欢 R 代码。对于情节,我使用 Plots.jl
我希望有人能够帮助我,如果您需要更多信息,请告诉我。
干杯