Warning: file_get_contents(/data/phpspider/zhask/data//catemap/4/r/73.json): failed to open stream: No such file or directory in /data/phpspider/zhask/libs/function.php on line 167

Warning: Invalid argument supplied for foreach() in /data/phpspider/zhask/libs/tag.function.php on line 1116

Notice: Undefined index: in /data/phpspider/zhask/libs/function.php on line 180

Warning: array_chunk() expects parameter 1 to be array, null given in /data/phpspider/zhask/libs/function.php on line 181
R 如何在Ggplot 2中的一个图上绘制多个逻辑回归曲线_R_Ggplot2_Regression_Logistic Regression - Fatal编程技术网

R 如何在Ggplot 2中的一个图上绘制多个逻辑回归曲线

R 如何在Ggplot 2中的一个图上绘制多个逻辑回归曲线,r,ggplot2,regression,logistic-regression,R,Ggplot2,Regression,Logistic Regression,我的数据框示例如下: 我想做3个图,根据预测概率为每个二进制文件(sent、corp、smb)绘制一条最合适的线——如果可能的话,我还想为预测概率添加点。到目前为止,我已经能够创建3个单独的绘图,但我想把所有三个放在一个绘图上。以下是我目前的情况: 以下是公司地块的代码: corp.line你无法绘制逻辑回归得到的“S”形曲线,因为你没有一个连续变量来绘制。相反,您只能绘制预测值,并围绕这些预测值绘制CI 在数据框中创建包含ent、corp和smb的列 newdata3<-read.tab

我的数据框示例如下:

我想做3个图,根据预测概率为每个二进制文件(sent、corp、smb)绘制一条最合适的线——如果可能的话,我还想为预测概率添加点。到目前为止,我已经能够创建3个单独的绘图,但我想把所有三个放在一个绘图上。以下是我目前的情况:

以下是公司地块的代码:


corp.line你无法绘制逻辑回归得到的“S”形曲线,因为你没有一个连续变量来绘制。相反,您只能绘制预测值,并围绕这些预测值绘制CI

在数据框中创建包含ent、corp和smb的列

newdata3<-read.table("clipboard", header=T)
newdata4<-unique(newdata3)[-4,] #different lower limits for smb... removing the second smb LL


newdata4$NewVar<-rep("",length(newdata[,1]))
newdata4$NewVar[which(newdata3$ent==1)]<-"ent"
newdata4$NewVar[which(newdata3$corp==1)]<-"corp"
newdata4$NewVar[which(newdata3$smb==1)]<-"smb"

windows(5,5)
ggplot(newdata4, aes(NewVar, PredictedProb, colour=NewVar)) + geom_point() +
    geom_errorbar(aes(ymin=LL, ymax=UL), width=.1, size=1)

newdata3可能有助于解决this@ckluss谢谢你的链接!虽然它没有解决我的问题,但它在如何组合情节方面确实很有帮助。嘿,谢谢你的帮助。在我考虑了更多之后,我得出了同样的结论:没有一个连续的预测因子。另外,如果您有时间,可以解释一下哪些命令和独特的命令正在执行什么操作吗?再一次,谢谢
unique()
,应用于整个数据帧,删除所有重复的行。如果你不这样做,你会一遍又一遍地画同样的东西,这可能是一个问题,也可能不是一个问题,你的图形如何结束。运行
unique(newdata3)
命令,查看
R
返回的内容。使用
which()
函数,我查看了您的每个列(ent、corp和smb),发现哪些行等于1,并将该列名指定给我的
NewVar
中的特定行。我解开了你的指示变量。运行
which(newdata3$ent==1)
查看返回的内容。还可以查看每个函数的
帮助()。
newdata3<-read.table("clipboard", header=T)
newdata4<-unique(newdata3)[-4,] #different lower limits for smb... removing the second smb LL


newdata4$NewVar<-rep("",length(newdata[,1]))
newdata4$NewVar[which(newdata3$ent==1)]<-"ent"
newdata4$NewVar[which(newdata3$corp==1)]<-"corp"
newdata4$NewVar[which(newdata3$smb==1)]<-"smb"

windows(5,5)
ggplot(newdata4, aes(NewVar, PredictedProb, colour=NewVar)) + geom_point() +
    geom_errorbar(aes(ymin=LL, ymax=UL), width=.1, size=1)