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 study notes[1]
r语言
Chef_Chen2 天前
从0开始学习R语言--Day57--SCAD模型
开发语言·学习·r语言
医工交叉实验工坊2 天前
R 语言绘制六种精美热图:转录组数据可视化实践(基于 pheatmap 包)
开发语言·信息可视化·r语言
AAIshangyanxiu2 天前
最新基于R语言结构方程模型分析与实践技术应用
开发语言·r语言·结构方程模型·生态统计学
biomooc3 天前
R拟合 | 一个分布能看到三个峰,怎么拟合出这三个正态分布的参数? | 高斯混合模型 与 EM算法
r语言
请你喝好果汁6413 天前
R中匹配函数
开发语言·r语言
Tiger Z3 天前
R 语言科研配色 --- 第 81 期 (附免费下载的配色绘图PPT)
r语言·科研·配色
恩喜玛生物6 天前
100个GEO基因表达芯片或转录组数据处理27 GSE83456
深度学习·r语言
Yolo566Q7 天前
基于现代R语言【Tidyverse、Tidymodel】的机器学习方法
开发语言·经验分享·机器学习·r语言
Tiger Z7 天前
《R for Data Science (2e)》免费中文翻译 (第0章) --- Introduction
r语言·数据科学·中文翻译