code: def _compare_full_firth_beta(test_data, golden_firth_beta):
      beta_init = np.zeros(test_data.covariates.shape[1])
      X = test_data.covariates
      y = test_data.phenotypes
      offset = test_data.offset
  
      test_firth_fit = af._fit_firth(
          beta_init=beta_init, X=X, y=y, offset=offset
      )
      test_firth_beta = test_firth_fit.beta
      assert np.allclose(golden_firth_beta, test_firth_beta)
  
imports: from dataclasses import dataclass
  import functions
  import glow.gwas.log_reg
  import glow.gwas.approx_firth
  import pandas
  from nptyping import Float, NDArray
  import numpy
  import pytest
  from typing import Any
label_desc: Use image node_id 2 for calling a global function with example usage: _compare_full_firth_beta(test_data, golden_firth_beta) without return types