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])
相关推荐
Biomamba生信基地1 小时前
R语言基础| 基本统计分析
开发语言·r语言·生信
smx66666681 天前
在docker中运行R容器,并在Windows下的vscode中使用该R
windows·docker·r语言
清同趣科研1 天前
扩增子测序|R包microeco妙用之零模型计算群落确定性和随机性过程(NST,pNST,βNTI,RCbray生态过程计算)
开发语言·r语言
拓端研究室TRL1 天前
R-INLA实现绿地与狐狸寄生虫数据空间建模:含BYM、SPDE模型及PC先验应用可视化...
开发语言·数据库·r语言
smx66666682 天前
在 Windows 下的 Docker 中安装 R语言
windows·docker·r语言
米糕.2 天前
【R语言】ggplot2绘图常用操作
大数据·开发语言·数据分析·r语言
drebander2 天前
使用 Conda 管理 R 语言
开发语言·r语言·conda
十三画者4 天前
【R包】tidyplots----取代ggplot2的科研绘图利器
r语言·数据可视化
周周记笔记4 天前
R与RStudio简介及安装
开发语言·r语言