多项 Logistic 回归，使用 ggplot2 无图形显示

Question

我从一项针对女性的调查中获得了以下数据。当我遵循 [http://www.ats.ucla.edu/stat/r/dae/mlogit.htm] 显示的指南时，我运行了一些命令并获得了系数、p 值和相对风险 我很欣赏他们。

我带来了以下命令

但是，我无法显示图表。

有人可以帮帮我吗？ 谢谢。

这是数据

年龄：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 551 456 35 52 35 60 54 46 32 52 40 60 70 50 53 42 60 70 70 54 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 Trad Islam Trad Chris Chris Chris Chris Chris Chris Chris Chris Chris Chris Chris Chris Islam Chris Chris Chris Chris Chris Islam Chris Trad Chris Chris Islam Chris Chris Chris Chris Chris Chris Chris克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·伊斯拉姆·克里斯·克里斯·克里斯·克里斯·克里斯·伊斯拉姆·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·克里斯·特拉德·克里斯·C克里斯克里斯克里斯 克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯

宠物：三三二三三一二一三一二三三三二三三三三三三三一三三三三三三二二三三三二三三三三一一三三三一 三 三 一 一 三 三 一 一 三 三 三 三 二 三 三 三 三 二 一 一 二 二 一 二 一 一 二 一 一 一 二 一 三 一 二 一 二 二 三 一 二 二 二 二 二 二 三 二 二二 二 二 一 一 三 一 三 三 三 三 一 三 三 二 三 三 三 三 三 一 三 二 一 二 三 三 三 二 二 一 三 二 三 二 三 二 三 二 三 三 三 三 二 三 二 三 三 二 二三 三 二 三 三 三 三 三 三 三 三 三 二 三 三 三 三 T三

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

Answer 1

我没有看到正确构建 ggplot 对象所需的加号。应该是：

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

仔细注意三个东西是如何与加号 (+) 相加的。

您可能还应该在脚本中显式打印ggplot 对象：

 print(Graph)

Answer 2

试试这个（你不需要连续变量'age'作为ggplot aes（）中的组）：

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

Graph