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])
相关推荐
Tiger Z2 天前
R 语言科研绘图第 20 期 --- 箱线图-配对
开发语言·程序人生·r语言·贴图
Coovally AI模型快速验证3 天前
目标检测新视野 | YOLO、SSD与Faster R-CNN三大目标检测模型深度对比分析
人工智能·yolo·目标检测·计算机视觉·目标跟踪·r语言·cnn
秦灏3 天前
十一、apply家族(4)
r语言
m0_748234084 天前
差异基因富集分析(R语言——GO&KEGG&GSEA)
开发语言·golang·r语言
xiao5kou4chang6kai44 天前
基于R语言的现代贝叶斯统计学方法(贝叶斯参数估计、贝叶斯回归、贝叶斯计算实践过程
回归·r语言·贝叶斯·线性回归·统计学
optiz5 天前
R语言绘图
开发语言·r语言
Biomamba生信基地5 天前
如何在没有root权限的情况下使用R语言
开发语言·r语言
码农老关【关东升】6 天前
告别 Excel,拥抱 R 语言:开启数据分析新时代
数据分析·r语言·excel
Biomamba生信基地6 天前
R语言基础| 方差分析
开发语言·r语言·生信
weixin_贾7 天前
R语言贝叶斯方法在生态环境领域中的高阶技术
r语言