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])
相关推荐
zhangfeng11339 小时前
BiocManager下载失败 R语言 解决办法
开发语言·r语言
Tiger Z3 天前
《R for Data Science (2e)》免费中文翻译 (第7章) --- Data import(1)
r语言·数据科学·中文翻译
魔力之心3 天前
R notes[2]
开发语言·r语言
Tiger Z5 天前
R 语言科研绘图第 71 期 --- 散点图-边际
r语言·论文·科研·绘图·研究生
小丑尽欢5 天前
R语言根据经纬度获得对应样本的省份
开发语言·r语言
czhc11400756637 天前
Linux 830 shell:expect,ss -ant ,while IFS=read -r line,
linux·运维·r语言
Morpheon7 天前
Intro to R Programming - Lesson 4 (Graphs)
开发语言·r语言
医工交叉实验工坊10 天前
R 语言 ComplexUpset 包实战:替代 Venn 图的高级集合可视化方案
算法·r语言
青春不败 177-3266-052010 天前
R语言贝叶斯方法在生态环境领域中的高阶技术应用
r语言·贝叶斯·生态学·科研绘图·结构方程·环境科学·混合效应
love530love11 天前
怎么更新 cargo.exe ?(Rust 工具链)
人工智能·windows·python·rust·r语言