【R语言】多项式逻辑回归logistic

require(foreign)
require(nnet)
require(ggplot2)
require(reshape2)# 失效下载本地数据
ml <- read.dta("https://stats.idre.ucla.edu/stat/data/hsbdemo.dta")
# 初步描述 社会经济地位
with(ml, table(ses, prog))
# 写作和选择倾向
with(ml, do.call(rbind, tapply(write, prog, function(x) c(M = mean(x), SD = sd(x)))))
# 设置参考组/对照组
ml$prog2 <- relevel(ml$prog, ref = "academic")
# final negative log-likelihood 179.981726
test <- multinom(prog2 ~ ses + write, data = ml)
summary(test)
# 系数除以标准误等于t值
z <- summary(test)$coefficients/summary(test)$standard.errors;z
# 2-tailed z test
p <- (1 - pnorm(abs(z), 0, 1)) * 2;p
# extract the coefficients from the model and exponentiate
# 从模型中提取系数并求指数相当于OR值
exp(coef(test))#  【predicted probabilities】
# fitted（拟合）是在*给定样本上*做预测，
#  而predict（预测）是在*新的样本上*做预测。# 使用预测的概率来帮助您理解模型{1}
head(pp <- fitted(test))
# 创建一个变量变化的小数据集，同时保持另一个常数
dses <- data.frame(ses = c("low", "middle", "high"), write = mean(ml$write))
predict(test, newdata = dses, "probs")# 使用预测的概率来帮助您理解模型{2}
# 使用预测概率理解模型的另一种方法是
# 查看每个Ses级别内连续预测变量写入的不同值的平均预测概率
summary(ml$write)
dwrite <- data.frame(ses = rep(c("low", "middle", "high"), each = 41), write = rep(c(30:70),3))
pp.write <- cbind(dwrite, predict(test, newdata = dwrite,type = "probs", se = TRUE))## calculate the mean probabilities within each level of ses
by(pp.write[, 3:5], pp.write$ses, colMeans)# 宽数据转长
lpp <- melt(pp.write, id.vars = c("ses", "write"), value.name = "probability")
head(lpp)
# 作图看概率趋势
ggplot(lpp, aes(x = write, y = probability, colour = ses)) + geom_line() + facet_grid(variable ~., scales = "free")

【R语言】多项式逻辑回归logistic相关推荐

R语言使用逻辑回归分类算法
R语言使用逻辑回归分类算法逻辑回归属于概率统计的分类算法模型的算法,是根据一个或者多个特征进行类别标号预测.在R语言中可以通过调用logit函数执行逻辑回归分类算法并预测输出概率.通过调用glm函数 ...
R语言用逻辑回归建立用户付费模型
原文链接:http://tecdat.cn/?p=967 对于某企业新用户,会利用大数据来分析该用户的信息来确定是否为付费用户,弄清楚用户属性,从而针对性的进行营销,提高运营人员的办事效率(点击文末& ...
r语言的逻辑回归分类
iris 是r语言内置的数据集 head(iris) # 与python的不同iris.head() Sepal.Length Sepal.Width Petal.Length Petal.Width ...
R语言惩罚逻辑回归、线性判别分析LDA、广义加性模型GAM、多元自适应回归样条MARS、KNN、二次判别分析QDA、决策树、随机森林、支持向量机SVM分类优质劣质葡萄酒十折交叉验证和ROC可视化
最近我们被客户要求撰写关于葡萄酒的研究报告,包括一些图形和统计输出. 介绍数据包含有关葡萄牙"Vinho Verde"葡萄酒的信息.该数据集有1599个观测值和12个变量,分别是 ...
R语言有序逻辑回归-因变量是等级资料
ordinal logistic regression适用于因变量为等级资料.使用课本例16-4的数据. 随机选取84例患者做临床试验,探讨性别和治疗方法对该病的影响.变量赋值为:性别(X1,男=0, ...
R语言LR逻辑回归实例
二分类实例去掉setosa类 index <- which(iris$Species == 'setosa') iris <- iris[- index,]training <-i ...
R语言主成分回归（PCR）、多元线性回归特征降维分析光谱数据和汽车油耗、性能数据...
原文链接:http://tecdat.cn/?p=24152 什么是PCR?(PCR = PCA + MLR)(点击文末"阅读原文"获取完整代码数据). • PCR是处理许多 x ...
R语言逻辑回归Logistic回归分析预测股票涨跌
最近我们被客户要求撰写关于逻辑回归的研究报告,包括一些图形和统计输出. 视频:R语言逻辑回归(Logistic回归)模型分类预测病人冠心病风险逻辑回归Logistic模型原理和R语言分类预测冠心病风 ...
R语言逻辑回归logistic regression对付费用户进行预测
对于某企业新用户,会利用大数据来分析该用户的信息来确定是否为付费用户,弄清楚用户属性,从而针对性的进行营销,提高运营人员的办事效率. 最近我们被客户要求撰写关于付费预测的研究报告,包括一些图形和统计输 ...

【R语言】多项式逻辑回归logistic

【R语言】多项式逻辑回归logistic相关推荐

最新文章

热门文章