逻辑回归卡方拟合优度检验 r语言

老师写的代码
r语言做逻辑回归时的皮尔森拟合优度检验

AZTAIDS <- expand.grid(AZT=factor(c("Yes","No"),levels=c("No","Yes")),Race=factor(c("White","Black"),levels=c("Black","White")))
AZTAIDS <- data.frame(AZTAIDS,Yes=c(14,32,11,12), No=c(93,81,52,43))
AZTAIDS <- AZTAIDS[,c(2,1,3,4)]
AZTAIDS#逻辑回归建模
myglm6 <- glm(cbind(Yes,No) ~ AZT + Race , family=binomial, data=AZTAIDS)
summary(myglm6)fitYes <- myglm6$fitted.values*(AZTAIDS$Yes+AZTAIDS$No)
fitNo <- (AZTAIDS$Yes+AZTAIDS$No)-fitYesfit <- c(fitYes,fitNo)
obs <- c(AZTAIDS$Yes,AZTAIDS$No)Gsquare <- 2*sum(obs*log(obs/fit));Gsquare
Xsquare <- sum((obs-fit)^2/fit);XsquareGXsquare <- function(model,data){fitYes <- model$fitted.values*(data$Yes+data$No)fitNo <- (data$Yes+data$No)-fitYesfit <- c(fitYes,fitNo)obs <- c(data$Yes,data$No)obs <- ifelse(obs==0,obs+10^(-8),obs)Gsquare <- 2*sum(obs*log(obs/fit));GsquareXsquare <- sum((obs-fit)^2/fit);Xsquaredata.frame(Gsquare,Xsquare)
}
GXsquare(model = myglm6,data = AZTAIDS)

逻辑回归卡方拟合优度检验 r语言相关推荐

UA MATH566 统计理论推导卡方拟合优度检验
UA MATH566 统计理论推导卡方拟合优度检验卡方拟合优度检验主要是检验categorical data的,假设一共有ddd种category,每一种理论比例为pip_ipi,满足 ∑i=1 ...
python卡方拟合优度检验_如何使用Python中的科学库执行卡方拟合优度检验？
假设我有一些我凭经验得到的数据: from scipy import stats size = 10000 x = 10 * stats.expon.rvs(size=size) + 0.2 * np ...
卡方 python_Python卡方拟合优度检验得到最优分布
给定一组数据值,我试图得到描述数据的最佳理论分布.经过几天的研究,我想出了下面的python代码.在import numpy as np import csv import pandas as pd ...
R语言逻辑回归、方差分析、伪R平方分析
目录怎么做测试假设条件并非所有比例或计数都适用于逻辑回归分析过度分散伪R平方测试p值 Logistic回归示例模型拟合系数和指数系数方差分析伪R平方模型的整体p值标准化残差图 ...
R语言glm.nb函数构建负二项分布回归模型、epiDisplay包的poisgof函数对拟合的负二项分布回归模型进行拟合优度检验、即模型拟合的效果、验证模型是否有过度离散overdispersion
R语言glm.nb函数构建负二项分布回归模型(negative binomial).使用epiDisplay包的poisgof函数对拟合的负二项分布回归模型进行拟合优度检验.即模型拟合的效果.验证模型 ...
r语言中残差与回归值的残差图_R语言逻辑回归、方差分析、伪R平方分析
原标题:R语言逻辑回归.方差分析.伪R平方分析原文链接:http://tecdat.cn/?p=9589 目录怎么做测试假设条件并非所有比例或计数都适用于逻辑回归分析过度分散伪R平方测试 ...
多元有序logistic回归_医学统计与R语言：多分类logistic回归HosmerLemeshow拟合优度检验...
微信公众号:医学统计与R语言如果你觉得对你有帮助,欢迎转发输入1:multinominal logistic regression install.packages("nnet" ...
r语言实现岭回归_预测分析：R语言实现
预测分析:R语言实现作者:(希)鲁伊·米格尔·福特(Rui Miguel Forte) 著出版日期:2016年10月文件大小:53.55M 支持设备: ￥40.00仅供试读适用客户端: 言商书 ...
毕业设计实用模型（四）——回归模型的实现（R语言）
目录 0引言 1.课本介绍 1.1理论的书 1.2 R语言的书 2.构造数据 3.相关性分析 4.多元回归模型的建立 4.1建立模型 5.2模型分析 5.3方差分析表 5.变量选择 5.1 逐步回归 ...
格兰杰因果关系检验r语言_R语言系列第四期：R语言单样本双样本差异性检验
之前详细介绍了利用R语言进行统计描述,详情点击:R语言系列第三期:①R语言单组汇总及图形展示.R语言系列第三期:②R语言多组汇总及图形展示.R语言系列第三期:③R语言表格及其图形展示从这个部分我们就 ...

逻辑回归卡方拟合优度检验 r语言

逻辑回归卡方拟合优度检验 r语言相关推荐

最新文章

热门文章