基于R的统计功效计算-SIMR包学习笔记2

---

title: “SIMR-2”
output: html_document

SIMR包可以进行各种假设检验的功效分析。函数doTest()可以用来对一个模型进行检验，以确定在进行功效分析之前，该test是合理的。使用**？tests**可以查看可以进行哪些假设检验。

?tests
library(simr)

1分类变量作为预测因子的二项GLMM

数据cbpp是有关牛传染性胸膜肺炎的数据。加入一个观测变量（1-56）以允许过度分散。响应变量可以使用cbind来定义。

data("cbpp")
str(cbpp)
cbpp$obs <- 1:nrow(cbpp) 
gm1 <- glmer(cbind(incidence, size - incidence) ~ period + (1|herd) + (1|obs), data=cbpp, family=binomial)
summary(gm1)$coef

变量period有4个水平，所以有3个哑变量，如果要同时检验3个哑变量，可以使用似然率检验lr

doTest(gm1, fixed("period", "lr"))

如果您只对变量period2感兴趣，您可以使用z检验。

doTest(gm1, fixed("period2", "z"))

如果您的模型加入了一个连续预测变量如size，您可以使用fixed选择检验哪一个固定效应。

gm2 <- glmer(cbind(incidence, size - incidence) ~ period + size + (1|herd), data=cbpp,family=binomial)
summary(