R语言线性模型glm()logistic回归模型

R语言广义线性模型glm()函数
glm(formula, family=family.generator, data,control = list(…))
formula数据关系,如y~x1+x2+x3
family:每一种响应分布(指数分布族)允许各种关联函数将均值和线性预测器关联起来。
常用的family:

binomal(link=’logit’) —-响应变量服从二项分布,连接函数为logit,即logistic回归

binomal(link=’probit’) —-响应变量服从二项分布,连接函数为probit

poisson(link=’identity’) —-响应变量服从泊松分布,即泊松回归

control:控制算法误差和最大迭代次数

glm.control(epsilon = 1e-8, maxit = 25, trace = FALSE)

 -----maxit:算法最大迭代次数,改变最大迭代次数:control=list(maxit=100)

源自
实战:

source("setwd.r")
##下面read.csv会出错,因为表头中文乱码,删掉表头,重新定义
Data=read.csv("bankloan.csv")[2:701,]
##查看前几行tail是后几行
head(Data)
colnames(Data)<-c("x1","x2","x3","x4","x5","x6","x7","x8","y")
tail(Data)
#logistic回归模型
glm<-glm(y~x1+x2+x3+x4+x5+x6+x7+x8,family = binomial(link = "logit"),data = Data)
glm
summary(glm)
###逐步寻优法 forward前向选择法backward后向选择法
logit.step<-step(glm,direction = c("both"))
summary(logit.step)
  • 10
    点赞
  • 111
    收藏
    觉得还不错? 一键收藏
  • 7
    评论
使用R语言glm函数可以拟合logistic回归模型glm函数是广义线性模型的函数,可以拟合多种类型的回归模型,包括logistic回归。 下面是一个示例代码,展示如何使用glm函数拟合logistic回归模型: ```R # 加载必要的包 library(ggplot2) # 创建一个示例数据集 data <- data.frame( x1 = rnorm(100), x2 = rnorm(100), y = sample(c(0, 1), 100, replace = TRUE) ) # 使用glm函数拟合logistic回归模型 model <- glm(y ~ x1 + x2, data = data, family = binomial) # 查看模型摘要 summary(model) # 进行预测 new_data <- data.frame( x1 = rnorm(10), x2 = rnorm(10) ) predictions <- predict(model, newdata = new_data, type = "response") # 可视化预测结果 ggplot(data, aes(x = x1, y = x2, color = as.factor(y))) + geom_point() + geom_contour(aes(z = predictions), alpha = 0.5) + theme_minimal() ``` 在上述代码中,首先加载了ggplot2包,并创建了一个示例数据集。然后使用glm函数拟合了logistic回归模型,其中y作为因变量,x1和x2作为自变量。在family参数中指定了binomial,表示拟合二分类模型。接着,可以使用summary函数查看模型的摘要信息。 最后,我们创建了一个新的数据集new_data,并使用predict函数进行预测。最后,使用ggplot2包可视化了原始数据和预测结果。 需要注意的是,以上只是一个简单的示例,实际应用中可能需要根据具体情况进行相应的数据预处理、模型优化等操作。
评论 7
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值