R 多项式逻辑回归，使用ggplot2无图形显示

我从一项妇女调查中获得了以下数据。我运行了一些命令并获得了系数、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 64 38 29 36 44 57 28 46 51 22 60 56 55 41 35 30 60 60 60 36 25 50 32 43 52 55 58 53 51 50 52 80 67 38 32 52 35 60 51 34 38 35 28 80 48 91 40 58 40 4045 37 38 25 58 53 42 60 70 68 62 53 53 66 40 39 35 32 25 40 45 51 48 44 47 52 57 57 92 88 87 44 40 47 41 42 49 40 37 40 40 38 38 38 38 32 28 34 32 40 39 38 38 38 36 36

rel:Chris Chris Chris Trad Islam Trad Chris Chris Chris Chris Chris Chris Chris Chris Chris Islam Chris Chris Chris Chris Chris Islam Chris Chris Chris Chris伊斯兰教Chris Trad Chris Islam Chris Chris Chris克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯克里斯

宠物：三三二三三三一二一三一二三三三三三三三三三三一三三三三三三二二三三三三三三三三三三三三三三三三三一一一一三三三三三一一一三ree一一三三三三三三三三三二一一二一二一二一一二一二三一二二三一二二二二二二二三二二二二二二二二一三三三三三三一三二三三三三三三三三三三三三三一三二一二三三三二二一三二三二三二三二三三三二三二三二三三二三三二三二三三二三三三三三三二三三三二三三三三三二三三三三三三三

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

Sec=read.table（“Second.txt”，header=TRUE）
秒
附加（第二节）
图书馆（nnet）
图书馆（GG2）
图书馆（E2）
m1我看不到正确构造ggplot对象所需的加号。它应该是：

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

请仔细注意三件事是如何加上加号（
+
）的

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

 print(Graph)

我看不到正确构造ggplot对象所需的加号。它应该是：

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

请仔细注意三件事是如何加上加号（
+
）的

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

 print(Graph)

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

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

以下是一些提示，用于发布格式良好的问题。以下是一些提示，用于发布格式良好的问题。