r - ゼロインフレ回帰モデルのna.actionを変更するには?

Question

パッケージの関数zeroinflを使用して、ゼロ膨張の負の二項回帰モデルを実行しています。pscl

分析の後半で従属変数に対して残差をプロットできるようにするために、モデルから NA を除外する必要があります。

したがって、設定したいと思いna.action="na.exclude"ます。たとえば、ゼロではない負の二項回帰モデル（パッケージglm.nbから使用）に対して問題なくこれを行うことができます。glm

fm_nbin <- glm.nb(DV ~ factor(IDV) + contr1
               +contr2 + contr3, data=df, 
               subset=(df$var<500), na.action="na.exclude")
fm_nbin.res = resid(fm_nbin) 
plot(fm_nbin.res~df$var)  

正常に動作します。ただし、ゼロ膨張モデルに対して同じことを行うと、機能しません。

zinfl <- zeroinfl(DV ~ factor(IDV) + contr1
               +contr2 + contr3 | factor(IDV) + contr1
               +contr2 + contr3, data=df, 
               subset=(df$var<500), na.action="na.exclude")
zinfl.res = resid(zinfl) 
plot(zinfl.res~df$var)

エラーを与える

Error in function (formula, data = NULL, subset = NULL, na.action = na.fail,  : 
  variable lengths differ (found for 'df$var')

回帰から NA を除外するために使用する必要がある他のコマンドはありますか?

編集： これは私が見つけた最も近い答えです。何らかの方法で私の問題に適用できますか? また、naresid何らかの形で適用できますか？