title: Predict ignores newdata argument in feglm estimation
body: <p>After estimating my model</p>
  <pre><code>train_data &lt;- df_est_hdfe[idx,]
  test_data &lt;- df_est_hdfe[-idx,]
  
  est &lt;- c(&quot;NATURAL&quot;, &quot;GDPsim&quot;, &quot;cap_lab_sim&quot;, &quot;cap_lab_sim_sq&quot;, &quot;cont&quot;)
  form &lt;- formula(paste0(&quot;rta ~ &quot;,paste(est, collapse = &quot; + &quot;), &quot;| exp_ind + imp_ind&quot;))
  
  train_data &lt;- as.data.frame(train_data)
  
  probithdfe &lt;- feglm(form, data = train_data[,c(&quot;rta&quot;, est, &quot;exp_ind&quot;, &quot;imp_ind&quot;)], family = binomial(link = &quot;probit&quot;))
  </code></pre>
  <p>I would like to use it to predict the output given some new data</p>
  <pre><code>pred_LM &lt;- predict(probithdfe, newdata = test_data[,c(est, &quot;exp_ind&quot;, &quot;imp_ind&quot;)], type = &quot;response&quot;)
  </code></pre>
  <p>However, I can't get my new data into the predict function as</p>
  <pre><code>dim(test_data)
  [1] 5910   13
  </code></pre>
  <p>while</p>
  <pre><code>length(pred_LM)
  [1] 11546
  </code></pre>
  <p>which is the size of the training data. How can I apply my estimated model to new data? I assume this might be specific to the feglm specification of the <code>alpaca</code> package.</p>