2017-04-05 44 views
0

我正在處理大量數據,只是使用glm在茱莉亞計算邏輯迴歸。我所做的是檢查就業率的不同影響,我有0或1. 現在我想繪製我的結果。我有gilman/hill書,他們在那裏做類似的練習,但在R(也是glm)。 他們繪製的是:朱莉婭:GLM陰謀logit函數

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的地塊

我希望有人能夠幫助我,如果您需要更多信息,請讓我知道。

乾杯

回答

1
using Plots 
plot(x, predict(fit1), seriestype = :line) #where x is your predictor variable 

更新以包括@pkofod的評論

+1

或只是 積(X,預測(FIT1)行=:線) – pkofod

+1

澄清:該行的默認值關鍵字是:路徑。路徑只會在第一對(x,y)到第二對,第三對,...最後一對之間劃一條線。如果你指定了line =:line,它將首先按x排序,這樣你就可以得到你的「平常」的線條圖,而不必先排序x和關聯的y。 – pkofod

+0

不錯的@pkofod! –