多項Logistic迴歸，不使用圖形顯示GGPLOT2

我已經從婦女的調查顯示獲得以下數據。我跑某些命令和I，接着[http://www.ats.ucla.edu/stat/r/dae/mlogit.htm] 我欣賞他們顯示的準則獲得的係數，P值和相對風險。多項Logistic迴歸，不使用圖形顯示GGPLOT2

我附帶以下命令

但是，我無法顯示圖表。

有人能幫助我嗎？謝謝。

這裏是數據

年齡：32 39 40 35 55 30 60 55 28 25 35 28 30 50 25 35 47 38 54 38 54 48 52 48 41 51 43 49 51 63 39 64 38 29 36 44 57 28 46 51 22 60 56 55 41 35 30 60 68 30 36 25 45 50 32 43 52 55 45 58 53 51 50 52 80 67 48 38 32 52 35 60 51 34 46 38 35 28 80 48 91 40 45 40 58 40 45 37 38 25 58 53 42 60 70 68 62 53 53 66 40 39 35 32 25 40 40 45 51 48 58 44 47 52 57 35 57 92 88 87 43 44 40 47 41 52 40 42 49 40 37 40 38 30 44 38 39 32 41 38 42 40 29 26 38 31 43 38 36 32 28 28 34 37 32 40 39 31 42 38 38 36

REL：克里斯克里斯克里斯克里斯島伊斯蘭島克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯Chris Chris Chris克里斯克里斯克里斯克里斯克里斯克里斯克里姆斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯Chris克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里姆斯克里斯克里斯島克里斯·克里斯 - 克里斯·克里斯島克里斯·克里斯克里斯·克里斯克里斯·克里斯克里斯·克里斯島克里斯·克里斯·克里斯 - 克里斯·克里斯 - 克里斯·克里斯 - 克里斯·克里斯 - 克里斯·克里斯 - 克里斯·克里斯 - 克里斯·克里斯 - 克里斯·克里斯·克里斯島島克里斯·克里斯·克里斯 - 克里斯·克里斯

寵物：三三二三三一T wo一三二二三三三二三三三三三三一三三三三三二二三三三二三三三一一一三個三個兩個三個三個三個三個兩個一個一個兩個兩個一個兩個一個一個兩個一個一個兩個一個三個一個兩個一個兩個兩個三個一個兩個兩個兩個兩個兩個兩個三個兩個兩個兩個兩個兩個一個三個一個三個三個三個三個一，三，二，三，三，三，三，三，兩，一，兩，三，三，兩，兩，三，兩，三，兩，三，兩，三，兩，三，三，三，兩，三，兩，三，三，兩，三，三，兩，三，三，三，三，三，三，三兩三三三三

Sec=read.table("Second.txt", header=TRUE) 
Sec 
attach(Sec)  
library(nnet) 
library(ggplot2) 
library(reshape2) 
m1 <- multinom(pet ~ age+rel, data = Sec) 
summary(m1) 
Relrisk=exp(coef(m1)) 
Relrisk 
z <- summary(m1)$coefficients/summary(m1)$standard.errors 
z 
p <- (1 - pnorm(abs(z), 0, 1))*2 
p 
frame=data.frame(age=rep(seq(22,92,by=0.434),3),rel=rep(c("Chris","Islam","Trad"),each=54)) 
frame 
attach(frame) 
pred=cbind(frame,predict(m1,newdata=frame,type="probs", se= TRUE)) 
pred 

Preprobmelt=melt(pred,id.vars=c("rel","age"),value.name="Probability") 
Preprobmelt 

Graph=ggplot(Preprobmelt,aes(x = age,y=Probability, colour=rel,group=age)) 

geom_line() 

facet_grid(variable~ .,scales = "free") 

Graph

No lines displayed

來源

2016-11-22 Nasamu Bawa

這裏是[下一次]幾個技巧（http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example）你想發佈一口井形成了問題。 –

我沒有看到正確構建ggplot對象的必要加號。它應該是：

Graph=ggplot(Preprobmelt,aes(x = age,y=Probability, colour=rel,group=age)) + 
    geom_line() + 
    facet_grid(variable~ .,scales = "free")

注意仔細如何三件事情正在用加號（+）標誌加在一起。

你或許應該也明確在腳本打印ggplot對象：

print(Graph)

來源

2016-11-22 17:40:25 Spacedman

剛剛嘗試這一點（你不需要連續可變「年齡」作爲ggplot AES組（））：

Graph=ggplot(Preprobmelt,aes(x = age,y=Probability, colour=rel)) + 
geom_line(lwd=1.2) + 
facet_grid(variable~ .,scales = "free") 

Graph

來源

2016-11-22 18:06:19

多項Logistic迴歸，不使用圖形顯示GGPLOT2

回答

相關問題