code: def compare_linreg_to_regenie(
      spark, output_prefix, missing=[], single_chr=True
  ):
      (
          genotype_df,
          phenotype_df,
          covariate_df,
          offset_df,
      ) = fx.get_input_dfs(
          spark, binary=False, missing=missing, single_chr=single_chr
      )
      glowgr_df = lr.linear_regression(
          genotype_df,
          phenotype_df,
          covariate_df,
          offset_df,
          values_column="values",
      ).toPandas()
      fx.compare_to_regenie(output_prefix, glowgr_df)
  
imports: import functions
  import numpy
  import statsmodels.api
  import pandas
  import glow.gwas.functions
  import glow.gwas.lin_reg
  import pytest
label_desc: Use image node_id 40 for calling a global function with example usage: compare_linreg_to_regenie(spark, output_prefix, missing, single_chr) without return types