Lme4：循环 glmer 时随机效应的错误

如何解决Lme4：循环 glmer 时随机效应的错误

我正在尝试对大量变量运行混合效应模型循环。

下面，我在一个空模型和一个具有一个预测变量的模型上测试了模型结构，并且这些模型没有错误或警告运行

## Null model 
df_m0 <- glmer(pred_recent_pf_bin ~ +   (1 | houseID),data = df,family = binomial,weights = weight,control = glmerControl(optimizer = "bobyqa"),nAGQ = 10)
#test one variable
df_mtest <- glmer(pred_recent_pf_bin ~  Motorsikal + (1 | houseID),nAGQ = 10)

但是，当我尝试循环模型时，我收到以下警告：

# set list of variables for looped model 
df_var <- df[2:253]

## Run model for all variables
for(i in 1:length(df_var)){
  
  # Set formula
  fmla <- paste0("pred_recent_pf_bin ~ + ",df_var[i]," + (1|houseID)")
  
  # Run model 
  m1 <- glmer(fmla,nAGQ = 10)
}

Error: couldn't evaluate grouping factor houseID within model frame: try adding grouping factor to data frame explicitly if possible

HouseID 设置为一个因子，有 2849 个级别

从我从类似问题 here 和 here 中可以看出，这可能是由于模型缺少数据参数，但我已确保将其包含在我的模型中。此错误的其他报告似乎具有嵌套随机效应，但我只有一个。

数据的最小可重现答案，给出相同的错误：

library(lme4)

#take subset of df dataset
rep <- subset(df,select = c("houseID","pred_recent_pf_bin","weight","smoke","Motorsikal","X100m__slopemean","ac200pa"))

#set list of variables for looped models
rep_vars <- rep[4:7]

#run model for all variables 
for(i in 1:length(rep_vars)){
  
  # Set formula
  fmla <- paste0("pred_recent_pf_bin ~ +  ",rep_vars[i],data = rep,nAGQ = 10)
}
Error: couldn't evaluate grouping factor houseID within model frame: try adding grouping factor to data frame explicitly if possible

> dput(rep[1:6,])
structure(list(houseID = structure(c(644L,2191L,823L,821L,368L,1274L),.Label = c("AB006","AB009","AB018","AB027","AB038","AB042","AB043","AB044","AB045","AB046","AB050","AB051","AB052","AB057","AB059","AB060","AB061","AB064","AB068","AB072","AB075","AJ001","AJ003","AJ004","AJ005","AJ006","AJ009","AJ010","AJ011","AJ012","AJ013","AJ014","AJ015","AJ016","AJ018","AJ020","BA002","BA003","BA007","BA009","BA010","BA012","BA014","BA015","BA016","BA017","BA018","BA020","BA026","BA027","BA029","BA032","BA033","BA034","BA039","BA041","BA042","BB002","BB003","BB004","BB005","BB008","BB009","BB010","BB011","BB012","BB014","BB017","BB019","BB020","BB021","BB022","BB025","BB026","BB027","BB028","BB029","BB030","BC004","BC009","BC010","BC011","BC012","BC013","BC014","BC015","BC016","BC022","BC023","BC025","BC026","BC027","BC028","BC029","BC030","BC034","BC037","BC040","BC041","BD001","BD003","BD004","BD008","BD009","BD013","BD015","BD016","BD017","BD020","BD022","BD023","BD024","BD026","BD027","BD029","BD030","BD031","BE001","BE002","BE004","BE005","BE006","BE007","BE008","BE009","BE010","BE014","BE015","BE018","BE019","BE020","BE021","BE023","BE024","BE025","BE026","BE028","BE029","BE032","BE033","BE034","BE035","BE037","BE038","BG002","BG003","BG004","BG005","BG008","BG009","BG010","BG011","BG012","BG013","BG015","BG016","BG018","BG019","BG020","BG022","BG023","BG024","BG026","BH001","BH002","BH003","BH004","BH006","BH007","BH008","BH009","BH010","BH012","BH016","BH018","BH019","BH020","BH026","BH028","BH030","BH031","BH032","BH033","BI001","BI005","BI006","BI009","BI013","BI022","BI024","BI025","BI031","BI033","BI036","BI043","BI046","BI050","BI052","BI059","BI061","BI066","BI067","BI073","BJ002","BJ003","BJ006","BJ007","BJ009","BJ010","BJ014","BJ015","BJ016","BJ019","BJ037","BJ039","BJ044","BJ045","BJ046","BJ047","BJ048","BJ049","BJ051","BK001","BK005","BK006","BK007","BK008","BK009","BK010","BK011","BK012","BM001","BM002","BM003","BM004","BM005","BM007","BM008","BM009","BM011","BM013","BM014","BM015","BM016","BM017","BM018","BM019","BM020","BM022","BM023","BM024","BO002","BO004","BO005","BO006","BO007","BO008","BO009","BO012","BO013","BO014","BO015","BO016","BO017","BP003","BP004","BP005","BP006","BP008","BP009","BP010","BP012","BP014","BP016","BP017","BP019","BP020","BP021","BP022","BP024","BP025","BP026","BP028","BP030","BP031","BR001","BR002","BR004","BR005","BR007","BR008","BR009","BR012","BR015","BR017","BR023","BR026","BR028","BR030","BR031","BR032","BR033","BR034","BR036","BR037","BS001","BS002","BS004","BS005","BS007","BS011","BS012","BS013","BS020","BS021","BS023","BS024","BS034","BS035","BS036","BS037","BS039","BS040","BS044","BS046","BS048","BT001","BT002","BT003","BT004","BT005","BT006","BU001","BU002","BU003","BU004","BU005","BU009","BU012","BU013","BW001","BW002","BW005","BW007","BW008","BW009","BW010","BW011","BW012","BW013","BW014","BW015","BW016","BW017","BW018","BX014","BX015","BX022","BX055","BX057","BX082","BX083","BX097","BX099","BX110","BX111","BX113","BX118","BX119","BX121","BX122","BZ001","BZ002","BZ004","BZ005","BZ006","BZ007","BZ008","BZ009","BZ010","BZ011","BZ012","BZ013","BZ014","DA008","DA012","DA027","DA034","DA035","DA039","DA040","DA055","DA064","DA065","DA066","DA071","DA078","DA089","DA103","DA105","DA113","DA123","DA130","DA139","DB001","DB002","DB004","DB006","DB012","DB017","DB025","DB028","DB039","DB041","DB045","DB048","DB049","DB053","DL003","DL004","DL005","DL007","DL008","DL012","DL013","DL014","DL016","DL017","DL018","DL019","DL020","DL021","DL022","DU005","DU008","DU010","DU011","DU013","DU018","DU020","DU021","DU022","DU029","DU032","DU033","DU036","DU043","DU045","DU048","DU049","GA002","GA005","GA006","GA007","GA008","GA009","GA010","GA011","GA013","GA015","GA016","GA017","GA023","GA024","GA026","GA027","GL001","GL002","GL003","GL019","GL020","GL021","GL022","GL023","GL024","GL026","GL028","GL029","GL030","GO001","GO005","GO007","GO008","GO009","GO013","GO015","GO016","GO018","GO019","GO023","GO025","GO028","GO029","GO030","GO032","GO034","GO037","GO038","GO039","GP001","GP002","GP003","GP004","GP005","GP006","GP007","GP008","GP009","GP010","GP011","GP014","GP015","GP016","GP017","GP019","GR001","GR002","GR003","GR004","GR007","GR008","GR009","GR010","GR011","GR012","GR013","GR014","GR015","GR016","GR017","GR018","GR020","GR021","GX001","GX002","GX003","GX004","KA005","KA007","KA008","KA009","KA010","KA017","KA018","KA019","KA020","KA021","KA022","KA026","KA028","KA031","KA033","KA034","KA036","KA040","KA041","KA045","KB001","KB002","KB003","KB004","KB005","KB006","KB007","KB008","KB009","KB010","KB012","KB013","KB014","KB015","KB016","KB017","KB018","KB019","KC002","KC003","KC004","KC005","KC006","KC008","KC010","KC013","KC015","KC020","KC021","KC022","KC023","KC024","KC025","KC026","KC027","KC042","KC045","KC049","KC053","KC055","KC056","KC058","KD001","KD002","KD003","KD005","KD007","KD008","KD020","KD021","KD022","KD027","KD029","KD030","KD032","KD038","KD040","KD041","KD042","KD044","KD045","KD046","KE001","KE002","KE003","KE004","KE005","KE006","KE007","KE008","KE009","KE010","KE011","KE012","KE013","KE014","KE015","KE016","KE017","KE018","KF001","KF002","KF003","KF004","KF005","KF006","KF007","KF008","KF009","KF010","KF011","KF012","KF014","KF015","KF016","KF017","KF018","KF019","KF020","KF021","KF023","KG001","KG003","KG004","KG005","KG006","KG007","KG009","KG010","KG011","KG012","KG013","KG014","KG015","KG017","KG018","KG019","KG020","KG021","KG022","KI001","KI002","KI004","KI005","KI007","KI010","KI011","KI012","KI013","KI014","KI015","KK001","KK004","KK006","KK009","KK010","KK011","KK013","KK014","KK015","KK016","KK017","KK018","KK020","KK022","KK026","KK028","KK033","KK034","KK035","KK036","KK037","KL001","KL003","KL005","KL008","KL010","KL014","KL015","KL017","KL018","KL021","KL024","KL026","KL027","KL029","KL036","KL040","KL043","KL045","KL047","KM007","KM008","KM009","KM012","KM013","KM015","KM017","KM019","KM020","KM021","KM024","KM025","KM027","KM029","KM031","KM034","KM036","KM037","KM039","KM040","KN001","KN004","KN005","KN007","KN008","KN009","KN011","KN013","KN014","KN015","KN016","KN017","KN018","KN020","KP002","KP003","KP012","KP020","KP042","KP043","KP047","KP054","KP064","KP065","KP075","KP078","KP081","KR001","KR002","KR003","KR004","KR005","KR007","KR008","KR009","KR010","KR011","KR013","KR014","KR016","KR018","KR022","KR023","KR024","KR029","KR030","KR032","KT001","KT002","KT003","KT004","KT005","KT006","KT007","KT008","KT010","KT011","KT012","KT013","KT014","KT015","KT016","KT017","KT018","KT019","KT020","KU001","KU002","KU004","KU007","KU011","KU012","KU014","KU015","KU016","KU018","KU019","KU027","KU028","KU030","KU031","KU034","KU035","KW004","KW005","KW006","KW007","KW008","KW009","KW010","KW012","KW013","KW014","KW015","KW017","KW018","KW019","KW020","KW022","KW024","KW026","KW028","KX004","KX006","KX007","KX009","KX010","KX011","KX012","KX022","KX024","KX027","KX029","KX031","KX033","KX034","KX036","KX037","KX038","KX043","KY001","KY002","KY004","KY005","KY006","KY007","KY008","KY009","KY010","KY011","KY012","KY013","KY014","KY015","LA001","LA003","LA012","LA013","LA016","LA017","LA027","LA028","LA029","LA030","LA032","LA037","LA039","LA043","LB001","LB003","LB004","LB007","LB008","LB010","LB011","LB012","LB013","LB014","LB015","LB016","LB017","LB018","LB019","LC002","LC004","LC005","LC006","LC007","LC008","LC009","LC010","LC011","LC012","LC013","LC015","LC016","LD006","LD012","LD014","LD016","LD023","LD026","LD027","LD029","LD031","LD034","LD036","LD037","LD038","LD039","LD041","LD042","LD051","LD052","LE003","LE004","LE005","LE006","LE007","LE008","LE012","LE013","LE016","LE017","LE019","LE021","LE022","LE023","LE024","LE027","LE028","LE029","LE030","LE031","LG004","LG008","LG012","LG016","LG017","LG018","LG019","LG023","LG024","LG026","LG027","LG029","LG030","LG035","LG036","LG037","LG044","LG047","LI002","LI005","LI006","LI009","LI010","LI012","LI014","LI016","LI019","LI023","LI024","LI026","LI027","LI029","LI037","LI038","LI040","LI043","LI044","LI046","LI054","LI057","LJ001","LJ002","LJ003","LJ012","LJ015","LJ016","LJ018","LJ020","LJ025","LJ028","LJ039","LJ041","LJ044","LJ047","LJ048","LJ051","LJ056","LJ061","LJ066","LK001","LK002","LL003","LL004","LL006","LL008","LL010","LL011","LL012","LL015","LL016","LL022","LL024","LL025","LL027","LL029","LL030","LL035","LL036","LL038","LM001","LM002","LM003","LM005","LM006","LM007","LM010","LM011","LM013","LM015","LM016","LM017","LM018","LN001","LN004","LN006","LN008","LN010","LN011","LN012","LN013","LN014","LN016","LN018","LN019","LN022","LN023","LN026","LN028","LN030","LN032","LN033","LN034","LN035","LO003","LO005","LO008","LO010","LO012","LO014","LO015","LO017","LO018","LO019","LO020","LO022","LO023","LO024","LO026","LO028","LO030","LO031","LP001","LP002","LP006","LP007","LP008","LP013","LP014","LP016","LP020","LP021","LP024","LP026","LP029","LP030","LP033","LP034","LP035","LP037","LP043","LR002","LR005","LR007","LR010","LR012","LR013","LR014","LR015","LR017","LR020","LR021","LR022","LR024","LR025","LR026","LR030","LR031","LR033","LR034","LR036","LS002","LS004","LS005","LS006","LS007","LS008","LS009","LS010","LS012","LS013","LS014","LS015","LS016","LS018","LS020","LS021","LS022","LS023","LS027","LS028","LT001","LT003","LT005","LT008","LT009","LT011","LT012","LT022","LT023","LT024","LT025","LT027","LX002","LX004","LX005","LX006","LX007","LX008","LX009","LX010","LX011","LX013","LX015","LX016","LX017","LX018","LX019","LX020","LX021","LX022","LX024","LX028","LY001","LY002","LY005","LY006","LY008","LY009","LY010","LY015","LY018","LY019","LY020","LY021","LY024","LY025","LY026","LY027","LY028","LY029","LZ001","LZ002","LZ003","LZ004","LZ005","LZ006","LZ007","LZ008","LZ009","LZ010","LZ011","LZ012","MA004","MA006","MA007","MA008","MA011","MA012","MA014","MA017","MA018","MA020","MA021","MA022","MA023","MA024","MA025","MA026","MA027","MA028","MA031","MA032","MB001","MB002","MB004","MB006","MB007","MB009","MB010","MB011","MB014","MB016","MB017","MB018","MB019","MB020","MB021","MB022","MB023","MC002","MC005","MC006","MC007","MC010","MC012","MC013","MC014","MC015","MC016","MC017","MC018","MC021","MC027","MC029","MC030","MC035","MC036","MD004","MD007","MD009","MD016","MD017","MD024","MD025","MD028","MD032","MD034","MD037","MD038","MD044","MD047","MD050","MD052","MD053","MD054","MD056","MD058","ME005","ME016","ME017","ME027","ME029","ME032","ME039","ME041","ME044","ME048","ME051","ME055","ME056","ME059","ME060","ME061","ME062","ME068","ME069","MF001","MF002","MF005","MF008","MF009","MF013","MF014","MF015","MF016","MF017","MF018","MF021","MF022","MF023","MF024","MF026","MF027","MF028","MF029","MF030","MF033","MF037","MF041","MF043","MF044","MF047","MF048","MF049","MF050","MF111","MF222","MF333","MG001","MG005","MG008","MG015","MG019","MG020","MG021","MG023","MG024","MG029","MG031","MG033","MG036","MG037","MG039","MG040","MG041","MG047","MG048","MG050","MG051","MH001","MH003","MH004","MH005","MI001","MI002","MI003","MI004","MI005","MI006","MI008","MI009","MI010","MI011","MI012","MI013","MI014","MI015","MI016","MI017","MI018","MI019","MI020","MK002","MK003","MK006","MK007","MK008","MK009","MK010","MK011","MK020","MK022","MK023","MK024","MK025","MK026","MK027","MK028","MK029","MK030","MK031","MK034","ML001","ML007","ML010","ML012","ML015","ML016","ML019","ML021","ML022","ML024","ML027","ML032","ML034","ML036","ML037","ML043","ML044","ML047","ML052","ML053","MM004","MM005","MM007","MM008","MM011","MM013","MM014","MM015","MM016","MM018","MM019","MM020","MM021","MM022","MM023","MM027","MM028","MM029","MN001","MN003","MN004","MN005","MN008","MN010","MN011","MN012","MN014","MN016","MN017","MN018","MN020","MN022","MN024","MN031","MN033","MN034","MN035","MN039","MN044","MN045","MN048","MN051","MN057","MN059","MP003","MP006","MP008","MP015","MP018","MP029","MP043","MP044","MP046","MP047","MP049","MP050","MP067","MP089","MP091","MP097","MP105","MP109","MQ004","MQ005","MQ007","MQ010","MQ013","MQ015","MQ018","MQ022","MQ025","MQ027","MQ028","MQ030","MQ032","MQ035","MQ040","MQ044","MQ045","MQ048","MQ049","MQ051","MR006","MR008","MR016","MR017","MR018","MR019","MR025","MR026","MR035","MR038","MR039","MR041","MU005","MU007","MU008","MU010","MU012","MU014","MU016","MU018","MU024","MU025","MU028","MU029","MU033","MU038","MU046","MX002","MX003","MX004","MX005","MX006","MX007","MX008","MX009","MX010","MX011","MX012","MX013","MX014","NA001","NA004","NA006","NA010","NA012","NA014","NA016","NA017","NA019","NA020","NA022","NA025","NA026","NA028","NA029","NA030","NA032","NA033","NA036","NA041","NB002","NB008","NB011","NB012","NB013","NB015","NB022","NB025","NB029","NB032","NB037","ND001","ND002","ND003","ND006","ND009","ND010","ND011","ND012","ND013","ND014","ND015","ND019","ND020","ND021","ND022","ND024","ND025","ND026","ND027","ND028","ND030","NP002","NP005","NP007","NP011","NP012","NP013","NP017","NP022","NP023","NP029","NP030","NP031","NP033","NP034","NP035","NP037","NU001","NU004","NU006","NU009","NU010","NU011","NU012","NU014","NU016","NU017","PA002","PA005","PA006","PA007","PA008","PA009","PA010","PA011","PA012","PA014","PA016","PA017","PA018","PA019","PB001","PB002","PB003","PB004","PB005","PB006","PB008","PB009","PB012","PB013","PB014","PB015","PB016","PB017","PB019","PB021","PB022","PB023","PB025","PB026","PC001","PC016","PC017","PC018","PC019","PC020","PC021","PC022","PC023","PC025","PC027","PC035","PC039","PC041","PC048","PC052","PC056","PC061","PC063","PD001","PD002","PD003","PD004","PD005","PD006","PD007","PD010","PD013","PD014","PD016","PD017","PD020","PD022","PD023","PD024","PD025","PD026","PE001","PE002","PE003","PE004","PE005","PE007","PE009","PE010","PE014","PE015","PE016","PE017","PE019","PE020","PE022","PE023","PE026","PE027","PG001","PG004","PG017","PG020","PG023","PG026","PG040","PG041","PG043","PG044","PG046","PG049","PH004","PH005","PI001","PI002","PI003","PI005","PI014","PI016","PI023","PI027","PI029","PI035","PI036","PK001","PK003","PK004","PK005","PK006","PK009","PK010","PK011","PK012","PK013","PK015","PK016","PK017","PK018","PK019","PK021","PK022","PL002","PL003","PL005","PL006","PL007","PL009","PL010","PL011","PL012","PL013","PM007","PM009","PM012","PM013","PM014","PM018","PM021","PM024","PM025","PM027","PM029","PM031","PM032","PM038","PM041","PM043","PM044","PM045","PM046","PN001","PN002","PN004","PN007","PN008","PN009","PN012","PN015","PN016","PN018","PN019","PN020","PN023","PN024","PN034","PN037","PN038","PN042","PN043","PO001","PO004","PO006","PO010","PO012","PO021","PO029","PO033","PO036","PO038","PO043","PO044","PO045","PO049","PO052","PO055","PO058","PO061","PO070","PP002","PP003","PP004","PP005","PP006","PP007","PP008","PP009","PP010","PP011","PP012","PP013","PP014","PP015","PP016","PP017","PP020","PP022","PP023","PR002","PR004","PR005","PR009","PR012","PR014","PR017","PR018","PR020","PR021","PR027","PR028","PR034","PR035","PR038","PR039","PR043","PR044","PR046","PS001","PS002","PS003","PS008","PS009","PS010","PS011","PS013","PS016","PS019","PS020","PS021","PS023","PS024","PS026","PS027","PS028","PS030","PS031","PS032","PT001","PT004","PT007","PT010","PT011","PT013","PT014","PT015","PT016","PT017","PT022","PT028","PT029","PT031","PT032","PT033","PT035","PT036","PT037","PU001","PU002","PU003","PU005","PU006","PU007","PU008","PU009","PU010","PW002","PW003","PW007","PW011","PW012","PW017","PW018","PW019","PW020","PW021","PW024","PW026","PW031","PW033","PW038","PW043","PW050","PW062","PW068","PX001","PX003","PX005","PX006","PX007","PX008","PX009","PX010","PX012","PX014","PX015","PX016","PX017","RB002","RB004","RB006","RB007","RB008","RB012","RB013","RB014","RB015","RB017","RB018","RB021","RB022","RB026","RB027","RB032","RB033","RB034","RB035","RB037","RG001","RG003","RG007","RG009","RG010","RG012","RG013","RG014","RG015","RG016","RG017","RG018","RG021","RG022","RG023","RG024","RG025","RG026","RG027","RG028","RG029","RK002","RK003","RK004","RK005","RK012","RK013","RK018","RK019","RK020","RK022","RK024","RK026","RK027","RK028","RK029","RK030","RK032","RK036","RK038","RK039","RK041","RM001","RM003","RM005","RM007","RM009","RM010","RM011","RM012","RM013","RM015","RM016","RM017","RM018","RM019","RM020","RM021","RU001","RU002","RU004","RU007","RU008","RU011","RU014","RU015","RU016","RU017","RU020","RU021","RU022","RU024","RU025","RU026","RU027","RU028","SA001","SA004","SA005","SA006","SA007","SA009","SA014","SA016","SA017","SA021","SA024","SA027","SA028","SA029","SA031","SA045","SA047","SA052","SA058","SA059","SA069","SA070","SB001","SB004","SB012","SB013","SB015","SB016","SB018","SB019","SB021","SB022","SB023","SB024","SB029","SB030","SB033","SB034","SB035","SB037","SB040","SC001","SC002","SC005","SC006","SC008","SC015","SC017","SC018","SC022","SD004","SD005","SD010","SD015","SD016","SD018","SD023","SE006","SE013","SE014","SE015","SE016","SE017","SE018","SE019","SE020","SE021","SF002","SF003","SF005","SF007","SF009","SF010","SF011","SF012","SF013","SF018","SF020","SF024","SF026","SF029","SF031","SF034","SF035","SG001","SG002","SG005","SG006","SG007","SG011","SG012","SG013","SG014","SG016","SG018","SG019","SG024","SG025","SG027","SG031","SG033","SG035","SG036","SG038","SH004","SH005","SH006","SH008","SH010","SH014","SH016","SH019","SH023","SH024","SH025","SH026","SH027","SH028","SH030","SH031","SI007","SI015","SI016","SI046","SI047","SI053","SI059","SI060","SI076","SI079","SI080","SI085","SI093","SI094","SI098","SI109","SI124","SJ002","SJ003","SJ009","SJ014","SJ015","SJ018","SJ021","SJ022","SJ024","SJ028","SJ029","SJ033","SJ038","SJ039","SJ040","SJ044","SJ045","SK008","SK011","SK013","SK018","SK024","SK028","SK029","SK033","SK034","SK036","SK039","SK041","SK043","SK046","SK047","SK049","SK052","SK058","SK059","SK060","SL001","SL003","SL005","SM002","SM004","SM007","SM008","SM010","SM011","SM016","SM017","SM018","SM019","SM020","SM021","SM022","SM023","SM025","SM026","SN001","SN002","SN005","SN006","SN008","SN009","SN010","SN013","SN014","SN015","SN016","SN017","SN018","SN021","SN025","SN026","SN027","SN028","SN029","SO001","SO010","SO011","SO013","SO014","SO016","SO017","SO022","SO024","SO027","SO028","SO030","SO033","SO046","SO047","SP001","SP003","SP004","SP006","SP007","SQ001","SQ002","SQ003","SQ004","SQ005","SQ006","SQ007","SQ008","SQ009","SQ010","SQ011","SQ012","SR001","SR002","SR003","SR004","SR006","SR007","SR008","SR009","SR010","SR011","SR013","SR015","SR016","SR017","SR018","SR019","SR021","SR022","SR024","SS001","SS003","SS012","SS013","SS014","SS016","SS023","SS024","SS026","SS027","SS030","SS031","SS035","SS036","SS041","ST004","ST005","ST010","ST012","ST013","ST020","ST022","ST025","ST026","SU001","SU002","SU003","SU004","SU006","SU007","SU008","SU009","SU010","SU012","SU013","SU014","SU015","SV001","SV003","SV004","SV005","SV006","SV008","SV010","SV011","SV012","SV014","SV015","SV016","SV017","SV019","SV020","SV021","SV022","SV023","SV024","SV025","SW001","SW002","SW003","SW005","SW006","SW007","SW008","SW009","SW010","SW011","SW012","SW013","SW014","SW015","SW016","SW018","SW019","SW021","SW022","SX002","SX003","SX004","SX005","SX006","SX007","SX008","SX009","SX010","SX012","SX013","SX015","SX016","SX017","SX018","SX019","SX020","SX023","SY001","SY002","SY011","SY013","SY014","SY015","SY018","SY024","SY025","SY029","SY031","SY034","SY035","SY036","SY040","SY041","SY046","SY048","SY053","SY059","SZ001","SZ002","SZ004","SZ006","SZ009","SZ011","SZ012","SZ013","SZ014","SZ015","SZ016","SZ017","SZ018","SZ019","SZ020","SZ021","SZ025","TA002","TA003","TA004","TA005","TA007","TA008","TA009","TA011","TA012","TA013","TA014","TA015","TA016","TA017","TA018","TA019","TA021","TA022","TA023","TB001","TB002","TB003","TB004","TB006","TB010","TB012","TB013","TB014","TB015","TB016","TB017","TB018","TB019","TB021","TB023","TC009","TC010","TC011","TC012","TC013","TC016","TC017","TD003","TD006","TD012","TD014","TD015","TD016","TD017","TD019","TD020","TD024","TD026","TD027","TD029","TD032","TD033","TD035","TD040","TD041","TD042","TE001","TE003","TE007","TE008","TE010","TE021","TE022","TE026","TE027","TE031","TE033","TE039","TE040","TE047","TE052","TE059","TE069","TE072","TE074","TF002","TF006","TF012","TF018","TF020","TF025","TF028","TF032","TF034","TF043","TF049","TF052","TF053","TF054","TF055","TF058","TF060","TG005","TG006","TG007","TG012","TG018","TG019","TG020","TG023","TG025","TG033","TG034","TG035","TG048","TG049","TG055","TG056","TH003","TH005","TH006","TH007","TH008","TH009","TH011","TH013","TH014","TH015","TH017","TH020","TH024","TI001","TI002","TI003","TI005","TI007","TI008","TI011","TI013","TI015","TI017","TI018","TI019","TI020","TI021","TI022","TI025","TI028","TI030","TI031","TI034","TI035","TI036","TI037","TI039","TI040","TI041","TI042","TI043","TJ002","TJ003","TJ004","TJ009","TJ010","TJ013","TJ019","TJ020","TJ023","TJ026","TJ029","TJ038","TJ040","TJ045","TJ047","TJ051","TJ052","TL002","TL003","TL004","TL006","TL007","TL008","TL009","TL010","TL011","TL013","TL014","TL015","TL016","TL018","TL019","TL020","TL021","TM002","TM003","TM004","TM006","TM007","TM011","TM012","TM014","TM018","TM019","TM020","TM021","TM022","TM023","TM027","TM028","TM029","TM030","TM031","TN006","TN008","TN009","TN011","TN013","TN017","TN020","TN023","TN024","TN028","TN031","TN033","TN035","TN039","TN053","TN054","TN059","TN062","TO001","TO002","TO003","TO005","TO006","TO008","TO009","TO010","TO011","TO012","TO013","TO014","TO015","TO016","TO018","TO019","TP001","TP003","TP004","TP005","TP006","TP008","TP010","TP011","TP012","TP015","TP016","TP017","TP019","TP020","TP025","TP041","TP044","TP055","TR002","TR003","TR008","TR011","TR012","TR013","TR018","TR020","TR023","TR024","TR027","TR028","TR029","TR031","TR032","TR033","TR040","TR043","TR044","TR047","TS001","TS002","TS003","TS005","TS007","TS015","TS017","TT001","TT002","TT003","TT005","TT006","TT007","TT008","TT009","TT010","TT011","TU002","TU003","TU004","TU006","TU007","TU010","TU012","TU015","TU017","TU018","TU021","TU024","TU026","TU033","TU036","TU037","TU038","TU039","TU041","TU044","TW001","TW002","TW003","TW005","TW007","TW009","TW010","TX001","TX002","TX005","TX012","TX018","TX019","TX025","TX027","TX028","TX033","TX035","TX036","TX040","TX041","TX042","TX044","TX047","US002","US003","US009","US011","US015","US016","US018","US020","US021","US023","US027","US029","US030","US031","US034","US035","US036","US038","WA008","WA009","WA018","WA019","WA021","WA026","WA028","WA029","WA039","WA043","WA049","WA053","WA054","WA058","WA060","WA062","WA068","WA069","WA072","WA076","WA077","WA079","WG001","WG002","WG007","WG008","WG009","WG013","WG016","WG017","WG019","WG020","WG021","WG022","WG023","WG026","WG028","WG029","WG032","WG034","WG036","WG037","WG043","WL001","WL002","WL003","WL005","WL006","WL007","WL008","WL009","WL010","WL011","WL012","WL013","WL014","WL015","WL016","WL017","WL018","WL019"),class = "factor"),pred_recent_pf_bin = structure(c(1L,1L,2L,1L),.Label = c("0","1"),weight = c(1,3,2),smoke = structure(c(1L,2L),.Label = c("FALSE","TRUE"),Motorsikal = structure(c(1L,X100m__slopemean = c(-0.0255895725057277,-1.08949221424853,-0.809496920233557,-1.0486043943469,-0.799390887160726,0.255158983311231),ac200pa = c(-0.121035991288976,-0.121035991288976,-0.121035991288976)),row.names = c(NA,6L),class = "data.frame")

解决方法

在您的公式中，“~”后不需要“+”，而 df_var[i] 会引入您只需要列名的整个向量。如果您使用以下命令构建公式，您的代码将起作用：

fmla <- paste0("pred_recent_pf_bin ~ ",names(df_var[i])," + (1|houseID)")