使用simulate进行参数自举时遇到错误
使用simulate进行参数自举时遇到错误
我最近在给一堆简单模型实现参数自举功能,于是自己定义了一个bootstrap.para(model, data)函数——它接收一个模型以及该模型基于的对应数据,执行以下几个核心步骤:
- 提取模型的相关信息,方便后续的模型拟合操作以及创建用于存储结果的矩阵
- 创建空矩阵,为后续存储自举结果预留空间
- 循环执行
nboot次(也就是设定的自举迭代次数):每次模拟生成新的响应变量,重新拟合模型,再把新模型的相关信息存入结果矩阵中
下面是我写的这个函数内部的初始代码片段:
# fit the model (this would normally be provided to the function) model <- glm(survival ~ 1, data = current_data, family = "binomial") # get model family (poisson, binomial, gaussian) and response variable family_type <- family(model)$family response_variable <- as.character(best_model$formula[[2]]) # (原代码后续内容未完整提供)
备注:内容来源于stack exchange,提问作者Em Laskey
相关产品推荐
相关产品推荐

