Rを用いて「多重代入法で欠損値補完をしたのちに、GLM/GEEを行いたい」

R Studio version 1.3.1093を用いて統計解析を行なっています。
Rの使用については初心者で、参考書やネットで調べたのですが調べきれていなかったらごめんなさい。

前提・実現したいこと

データベース解析をするにあたって
①欠損値補完（多重代入法）
②GLM、可能であればGEEを行なって解析する
という２ステップをしたいです。

多重代入法をmice を使って行い、withで解析し、poolで統合するという手段が最適かと思って試みたのですが、lmは可能でもglmだとうまく解析できません。そして、本当はgeeがしたいです。

多重代入法のパッケージがmiceではないのを使った方がいいとか、glmのコードが間違っているのかとか解決法がわからなかったので質問させていただきました。

ここに質問の内容を詳しく書いてください。

library(mice)
imppt<- mice(selection_pt,
m=10,
maxit=50,
method="pmm",
printFlag=FALSE,
seed=500)
fit<- with(data = imppt,lm(motality~ISS+male))
pool(fit)
summary(pool(fit))

これは可能です。

発生している問題・エラーメッセージ

エラーメッセージ

lmではなくて、

fit1 <- with(data=imppt, glm(motality~ISS+male, family = binomial)

とすると

Error: Incomplete expression: fit1 <- with(data=imppt, glm(motality~ISS+male, family = binomial)

というエラーメッセージが出ます。

該当のソースコード

library(mice)
imppt<- mice(selection_pt,
m=10,
maxit=50,
method="pmm",
printFlag=FALSE,
seed=500)
fit1 <- with(data=imppt, glm(motality~ISS+male, family = binomial)