ロジスティック回帰分析

ロジスティック回帰分析（2値）

#データの読み込み
library(readxl)
dat1 <- read_excel("C:/Users/2ken/Desktop/Umeda3.xlsx")
View(dat1)
dat1=na.omit(dat1)

#基礎集計
library(psych)
describe(dat1)

#データ型の確認と変更
str(dat1)
dat1$graf=as.factor(dat1$graf)
dat1$subject=as.factor(dat1$subject)

#ロジスティック回帰
res1=glm(answer~graf+subject,data = dat1,family = "binomial")
summary(res1)
exp(res1$coefficients)

#決定係数
library(BaylorEdPsych)
PseudoR2(res1)

#ロジスティック回帰の適合度の検定
library(ResourceSelection)
hoslem.test(res1$y, res1$fitted.values)