code: def skat_dataset():
      ds2 = hl.import_vcf(resource("sample2.vcf"))
  
      covariates = hl.import_table(
          resource("skat.cov"), impute=True
      ).key_by("Sample")
  
      phenotypes = hl.import_table(
          resource("skat.pheno"),
          types={"Pheno": hl.tfloat64},
          missing="0",
      ).key_by("Sample")
  
      intervals = hl.import_locus_intervals(
          resource("skat.interval_list")
      )
  
      weights = hl.import_table(
          resource("skat.weights"),
          types={"locus": hl.tlocus(), "weight": hl.tfloat64},
      ).key_by("locus")
  
      ds = hl.split_multi_hts(ds2)
      ds = ds.annotate_rows(
          gene=intervals[ds.locus], weight=weights[ds.locus].weight
      )
      ds = ds.annotate_cols(
          pheno=phenotypes[ds.s].Pheno, cov=covariates[ds.s]
      )
      ds = ds.annotate_cols(
          pheno=hl.if_else(
              ds.pheno == 1.0,
              False,
              hl.if_else(ds.pheno == 2.0, True, hl.missing(hl.tbool)),
          )
      )
      return ds
  
imports: import hail
  import pytest
  from hail.utils import FatalError, HailUserError
  from ..helpers import resource, test_timeout, qobtest
label_desc: Use image node_id 10 for calling a global function with example usage: skat_dataset() and returns: ds