当前位置:网站首页>[R language] interaction test data
[R language] interaction test data
2022-08-09 06:01:00 【The bionic bug】
#载入数据
data<-read.csv("Labor analgesia and7Weekly postpartum prognostic score.csv")
str(data)
data$Score <- as.numeric(data$Score)
data$LaborAnalgesia <- as.factor(data$LaborAnalgesia)
data$Delivery <- as.factor(data$Delivery)
data$ID <- as.factor(data$ID)
#Simple chi-square test
data$Score_group <- ifelse(data$Score<=13,0,1)
table(data$Score_group,data$LaborAnalgesia)
s=chisq.test(data$Score_group,data$LaborAnalgesia,correct = TRUE)
s
s$expected
#Hierarchical Chi-square analysis 【按Delivery分层】
table(data$Score_group,data$LaborAnalgesia,data$Delivery)
mantelhaen.test(data$Score_group,data$LaborAnalgesia,data$Delivery)
#【广义线性模型】 【as a dependent variable】
fit <- glm(Score_group ~ LaborAnalgesia+Delivery,family = binomial(link = "logit"), data)
summary(fit)
#Consider interactions
fit2 <- glm(Score_group ~ LaborAnalgesia*Delivery,family = binomial(link = "logit"), data)
summary(fit2)
# 广义线性模型fit的R-Square
library(rsq)
rsq(fit, #模型
adj=FALSE #Whether to correctR2
,type=c("lr"))
rsq(fit2,adj=FALSE,type=c("lr"))
#【线性混合模型】【分娩方式Deliveryas a random effect】
library(lme4)
lmer1 <- lmer(Score_group ~ LaborAnalgesia+(1|Delivery),data )
summary(lmer1)
rsq(lmer1,adj=FALSE,type=c("lr"))
边栏推荐
猜你喜欢
随机推荐
正则表达式-判断字符串是否匹配“AABB”模式
list 字符串的输出方法 print(*a) print(““.join(str(c) for c in a) )
SiO2-NH2/Fe3O4表面氨基化的Fe3O4/SiO2复合磁性纳米粒子(齐岳试剂)
二硫化钼/二氧化铪的复合纳米材料(MoS2/HfO2)|钽掺杂二氧化铪纳米颗粒(齐岳bio)
RT201 Domestic PA RF Power Amplifier Compatible with RFX2401C
GNNExplainer应用于节点分类任务
三剑客进阶
四氧化三铁/硫化铋纳米复合材料([email protected]@BSABiS纳米颗粒)|树状大分子稳定的硫化铋纳米颗粒|科研试剂
22 high mid term paper topics forecast
Shell 文本三剑客 (grep、sed、awk)
51串口通信(上)
deploy发布的jar包中文乱码
Selection of MOS tube
untiy 倒计时
Kubernetes apparmor profile
kubernetes apparmor 简介
se2clam
AIOT介绍
tidb crash test
phpstudy install flarum forum