Normal Limit Check

1library(yamaa)23yamaa_spec("spec.yaml")$run()$output

1import yamaa23yamaa.spec("spec.yaml").run().output

1# Reference solution for the yamaa benchmark adam-adlb-reported-precision (R track).2#3# Written from the benchmark's full prompt and its inputs alone, as an agent4# would write it. Harbor's oracle agent runs it to check that the prompt can5# be solved and that the grader scores a correct answer 1.67library(dplyr, warn.conflicts = FALSE)8library(readr)910lb <- read_csv(11  "/app/input/lb.csv",12  col_types = cols(13    AVAL = col_double(),14    ANRLO = col_double(),15    .default = col_character()16  )17)1819# The result as a multiple of the lower limit, rounded once to four places20# with an exact half away from zero.21adlb <- lb |>22  mutate(23    ratio = AVAL / ANRLO,24    R2ANRLO = case_when(25      is.na(AVAL) | is.na(ANRLO) | ANRLO == 0 ~ NA_real_,26      ratio >= 0 ~ floor(ratio * 10000 + 0.5) / 10000,27      .default = ceiling(ratio * 10000 - 0.5) / 1000028    )29  ) |>30  arrange(USUBJID) |>31  select(STUDYID, USUBJID, PARAMCD, AVAL, ANRLO, R2ANRLO)3233dir.create("/app/output", showWarnings = FALSE)34write_csv(adlb, "/app/output/adlb.csv", na = "")

1# Reference solution for the yamaa benchmark adam-adlb-reported-precision (Python track).2#3# Written from the benchmark's full prompt and its inputs alone, as an agent4# would write it. Harbor's oracle agent runs it to check that the prompt can5# be solved and that the grader scores a correct answer 1.67import csv8from decimal import Decimal, ROUND_HALF_UP, getcontext910getcontext().prec = 501112FOUR_PLACES = Decimal("0.0001")131415def parse_number(text):16    text = (text or "").strip()17    if text == "":18        return None19    return Decimal(text)202122def report(value):23    """Format a Decimal to exactly four decimal places, rounding an exact24    half away from zero. A missing value is reported as absent."""25    if value is None:26        return ""27    return str(value.quantize(FOUR_PLACES, rounding=ROUND_HALF_UP))282930def main():31    with open("/app/input/lb.csv", newline="") as f:32        rows = list(csv.DictReader(f))3334    out_rows = []35    for row in rows:36        aval = parse_number(row["AVAL"])37        anrlo = parse_number(row["ANRLO"])38        if aval is None or anrlo is None or anrlo == 0:39            ratio = None40        else:41            # Nothing is rounded before the report: the ratio keeps every42            # digit it was calculated with and is rounded once, in report().43            ratio = aval / anrlo44        out_rows.append(45            {46                "STUDYID": row["STUDYID"],47                "USUBJID": row["USUBJID"],48                "PARAMCD": row["PARAMCD"],49                "AVAL": report(aval),50                "ANRLO": report(anrlo),51                "R2ANRLO": report(ratio),52            }53        )5455    fieldnames = ["STUDYID", "USUBJID", "PARAMCD", "AVAL", "ANRLO", "R2ANRLO"]56    with open("/app/output/adlb.csv", "w", newline="") as f:57        writer = csv.DictWriter(f, fieldnames=fieldnames)58        writer.writeheader()59        writer.writerows(out_rows)606162if __name__ == "__main__":63    main()

1/* opensas reference solution for the yamaa benchmark adam-adlb-reported-precision.2   Reads the benchmark inputs and derives the requested output.3   Execute with the existing opensas runtime: opensas /app/output/result.sas. */45data lb;6  length STUDYID USUBJID PARAMCD $1024;7  infile "/app/input/lb.csv" dsd dlm="," firstobs=2 truncover lrecl=32767;8  input STUDYID : $char1024. USUBJID : $char1024. PARAMCD : $char1024. AVAL ANRLO;9run;1011data result;12  set lb;13  if not missing(AVAL) and not missing(ANRLO) and ANRLO ne 0 then R2ANRLO=round(AVAL/ANRLO,0.0001);14run;1516proc export data=result(keep=STUDYID USUBJID PARAMCD AVAL ANRLO R2ANRLO) outfile="/app/output/adlb.csv" dbms=csv replace; run;