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])
相关推荐
烟锁池塘柳03 小时前
【R语言】R 语言中打印含有双引号的字符串时会出现 “\” 的原因解析
r语言
全栈开发圈4 天前
干货分享|如何从0到1掌握R语言数据分析
开发语言·数据分析·r语言
小杜的生信筆記7 天前
基于R语言,“上百种机器学习模型”学习教程 | Mime包
开发语言·学习·机器学习·r语言·sci
在打豆豆的小潘学长7 天前
【R语言】多样本单细胞分析_SCTransform+Harmony方案(2)
开发语言·r语言
TS的美梦7 天前
ROGUE: 【张院士团队R包】一种基于熵的用于评估单细胞群体纯度的度量标准
开发语言·r语言
weixin_493202639 天前
R语言代码加密(1)
r语言
Tiger Z9 天前
《R for Data Science (2e)》免费中文翻译 (第3章) --- Data transformation(2)
r语言·数据科学·中文翻译
星石传说11 天前
使用R将nc文件转换为asc文件或者tif文件
r语言·生信
Mister Leon12 天前
机器学习Adaboost算法----SAMME算法和SAMME.R算法
算法·机器学习·r语言
Tiger Z13 天前
R 语言科研绘图第 67 期 --- 箱线图-显著性
r语言·论文·科研·绘图·研究生