【问题标题】:Superimposing a logistic line of fit叠加逻辑拟合线
【发布时间】:2014-03-21 21:18:57
【问题描述】:

如何在绘图上叠加逻辑曲线。

Temp<-c(27.2,28.3,29.9) 
Temp 
male<-c(0,8,8) 
male 
female<-c(10,4,2)
female
table=read.table("E:\\Book1.txt",header=T)
attach(table) 
table 
Y=cbind(male,female) 
Y 
mylogit <- glm(Y ~ Temp, family = "binomial",table)
summary(mylogit) 

我需要为男性和女性叠加一条逻辑曲线

curve(predict(mylogit,data.frame(male=x),type="resp"),add=TRUE) 
title(main="Males Temperature with Fitted GLM Logistic Regression Line") 

任何帮助?

【问题讨论】:

    标签: r curve-fitting logistic-regression


    【解决方案1】:

    你的代码有点混乱,所以这可能不是你想要的。

    逻辑回归模拟事件发生的概率。因此,在您的情况下,您正在根据数据对males 的概率进行建模,该数据在三个不同的Temp 处给出males 的比例,例如0/10、8/12 (.66) 和 8/10 (0.8)。因此,要将模型与数据进行比较,您必须绘制预测响应与男性比例。

    Temp    <- c(27.2,28.3,29.9) 
    male    <- c(0,8,8) 
    female  <- c(10,4,2)
    Y       <- cbind(male,female)  
    mylogit <- glm(Y ~ Temp, family = "binomial")
    plot(Temp,predict(mylogit,type="resp"),
         type="b",col="blue",lty=2, 
         ylim=c(0,1),ylab="Fraction of Males",
         main="Males Temperature with Fitted GLM Logistic Regression Line")
    points(Temp,male/(male+female),pch=16, col="red")
    

    【讨论】:

    • 您的帖子很有帮助。但是,我相信您可以叠加曲线而不是最佳拟合线。我遇到了在以下链接上绘制曲线的代码。 cran.r-project.org/web/packages/HSAUR/vignettes/…。我也试过下面的代码,(curve(predict(Temp,male=x,type="resp"), add=TRUE)) 但它给了我一个错误。
    • 我承认代码令人困惑,尤其是在阅读表格时。我只是在尝试不同的解决方案。
    • 我设法使用以下代码叠加曲线。 ###问题1 a Temp
    猜你喜欢
    • 1970-01-01
    • 2019-10-13
    • 2018-10-02
    • 1970-01-01
    • 2021-08-07
    • 2019-07-14
    • 2021-01-13
    • 1970-01-01
    • 2015-04-26
    相关资源
    最近更新 更多