R语言详解predict函数

R语言中predict函数在建立模型,研究关系时常用。但是不同type得到的结果常常被混为一谈,接下来,探讨predict得到的不同结果。

复制代码
#数据
set.seed(123)
n<-1000
age<-rnorm(n,mean=50,sd=10)
gender<-rbinom(n,1,0.5)
disease<-rbinom(n,1,plogis(-3+0.05*age+0.5*gender))
data<-data.frame(age,gender,disease)

1、type="response"得到的时预测概率

复制代码
glm_model<-glm(disease~age,family=binomial(link="logit"),data=data)
a1=predict(glm_model,type="response")
print(a1[1:10])

2、type="link"得到log-odds (发生概率与不发生概率比值取自然对数)

复制代码
glm_model<-glm(disease~age,family=binomial(link="logit"),data=data)
a1=predict(glm_model,type="link")
print(a1[1:10])

3、type="terms"得到一个矩阵,给出线性预测量表上模型公式中每个项的拟合值。

复制代码
a3=predict(glm_model,type="terms")
print(a3[1:10])
相关推荐
探序基因1 天前
R语言读取单细胞转录组基因表达矩阵loom文件
开发语言·r语言
房开民3 天前
paddle 文本检测识别模型转为onnx
开发语言·r语言·paddle
全栈开发圈3 天前
干货分享|R语言聚类分析2
人工智能·机器学习·r语言
全栈开发圈4 天前
干货分享|R语言聚类分析1
开发语言·r语言
AI科技星7 天前
空间光速螺旋动力学:统一质量、引力、电磁与时空本源的公理化理论与全现象验证
c语言·开发语言·opencv·算法·r语言
zhangfeng113311 天前
提示 R for Windows front-end 怎么被防火墙 阻止了 Rscript.exe` 和 `R.exe`区别
windows·r语言·php
全栈开发圈11 天前
新书速览|R语言医学数据分析与可视化
开发语言·数据分析·r语言
木与长清12 天前
人鼠同源基因离线转换
数据库·矩阵·数据分析·r语言
HP-Patience12 天前
【Rmarkdown】快速入门
r语言
HP-Patience12 天前
【Data Mining】01抽样技术
人工智能·数据挖掘·r语言