0

我正在尝试在 GLM 模型上运行引导重采样,但我一直遇到错误,并且在网上找不到任何解决方案。

当我尝试使用归一化变量(sub.size_z)和转换为因子(小时)的连续变量时,我遇到的第一个问题是“变量长度不同”。我检查了我的数据集,发现任何地方都没有 NA,所以我不知道如何解决这个问题。 在此处输入图像描述

第二个错误是“错误:从基础到 coefs 列表的命名不匹配”,当我尝试使用 GLM 运行引导程序而没有转换我之前提到的变量时,就会发生这种情况。

在此处输入图像描述

谁能指出我的解决方案?非常感谢!

这是我使用的代码!

data <-read.delim2("C:/Users/BRIZ_/Desktop/analisi scratching/Hypothesis 3/Hypothesis 3/database.txt")


library(car)
library(lme4)
library(MASS)
library(dplyr)
library(blmeco)
library(MuMIn)
library(lmerTest)
library(AER)
library(DHARMa)

sub.size_z<-scale(sub.size, center=TRUE,scale=TRUE)
hours<- as.factor(hour)


modelboot<-glmmTMB(sc.count~ offset(log(DURATION)) + hours + general.activity + sub.size_z + approach5+ id,data=data, family=genpois(link = "log"))
summary(modelboot)

bootstrap_model(modelboot,base_data=data,resamples =100)
4

0 回答 0