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])
相关推荐
喜欢就别17 小时前
平均风险:Rényi DP 和零集中 DP
开发语言·r语言
Bioinfo Guy4 天前
高分Panel复现系列|甲基化分型可视化:怎么画这样的ComplexHeatmap分块热图
r语言·生物信息学·生信可视化
临床数据科学和人工智能兴趣组8 天前
R读数据总出错?一份文件读取全指南
r语言·r语言-4.2.1
GZ同学1 个月前
单双变量Ripley’s K函数 R 语言实现
开发语言·r语言
popcorn_min1 个月前
共享单车需求预测:时间特征工程 + 随机森林,R² 达到 0.931
随机森林·r语言·kotlin
m0_497048931 个月前
.NET10+Avalonia跨平台截屏工具解析
r语言
dongf20191 个月前
R语言朴素贝叶斯算法---iris数据集
开发语言·算法·数据分析·r语言
All_Will_Be_Fine噻1 个月前
重建R环境
开发语言·r语言
梦想的初衷~1 个月前
基于R语言APSIM模型高级应用及批量模拟
r语言·apsim·作物生长模型·农业模拟·气象数据处理·土壤水分平衡
dongf20191 个月前
R 语言 逻辑斯蒂回归
开发语言·数据分析·回归·r语言