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])
相关推荐
AC赳赳老秦8 小时前
科研数据叙事:DeepSeek将实验数据转化为故事化分析框架
开发语言·人工智能·数据分析·r语言·时序数据库·big data·deepseek
Kratzdisteln16 小时前
【MCM】((N_G,N_A)) 可行域; ((C_R,c_E)) 分区图
c语言·开发语言·r语言
张小凡vip17 小时前
数据挖掘(八)--让Jupyter notebook支持R语言环境
jupyter·数据挖掘·r语言
青灯照颦微2 天前
【R】三种方式安装R包
开发语言·r语言
weixin_462446232 天前
PaddleX 3.2 人脸识别实战:自定义人脸库 + CartoonFace 官方案例 Top-K 识别完整指南
开发语言·r语言
Tiger Z2 天前
《R for Data Science (2e)》免费中文翻译 (第19章) --- Joins(1)
r语言·编程·数据科学
Tiger Z2 天前
《R for Data Science (2e)》免费中文翻译 (第18章) --- Missing values
开发语言·r语言
带我去滑雪3 天前
R语言抑郁症状网络分析
r语言
AC赳赳老秦3 天前
R语言数据分析:DeepSeek辅助生成统计建模代码与可视化图表
开发语言·人工智能·jmeter·数据挖掘·数据分析·r语言·deepseek
czliutz7 天前
R语言gm玩音乐示例代码Rmarkdown
开发语言·r语言