weltby Ergodic Sign inGet started
All use cases

Is This Jar of Spread Contaminated?

From 30 sensor measurements taken on the line, does this jar of cocoa-hazelnut spread contain a foreign body? A missed jar means a recall, and a false alarm wastes product.

Use caseBiologyPredictionModel evaluation

The table

2,400 jars scanned on a production line with 30 sensor measurements each, labelled clean or contaminated. It's in Welt's research catalogue as Contaminant Detection in Cocoa Hazelnut Spread Jars, 2,400 rows by 31 columns, licensed CC BY-SA. The catalogue measures its signal at 0.81. At 2,400 rows and 30 features, TabICLv2 and Kumo Tabular Medium take the table as it is. The others can't take a table this size, so the walkthrough leaves them out.

Walkthrough

  1. Load it from the catalogue

    Find the table, read its recorded target and print the column names.

    import pandas as pd
    from welt import Client
    
    with Client() as client:
        item = client.research_datasets(q="contaminant detection cocoa hazelnut", limit=1)["items"][0]
        info = client.research_dataset(item["id"])
        data = client.load_research_dataset(item["id"])
    
    target = info["targets"][0]
    target = target["name"] if isinstance(target, dict) else target
    print([column["name"] for column in info["schema"]], target)
  2. Pick the rows

    Set aside 20 rows to predict at the end, and keep the rest for the models to learn from.

    known = data.dropna(subset=[target])
    new_rows = known.sample(20, random_state=1)
    table = known.drop(new_rows.index)
  3. Compare the models and predict

    The models compete on the same folds against a simple guess, then the top-ranked one is fitted on the whole table. ROC AUC runs from 0.5, a coin toss, to 1, a perfect ranking, and a lower Brier score means the probabilities can be trusted more. The last table puts each prediction next to the real answer.

    with Client() as client:
        evaluation = client.evaluate(table, target, task="classification", seed=7,
                                     models=["tabicl-v2", "kumo-medium"])
        model = evaluation.fit()
    
    print(evaluation.summary()[["name", "rank", "roc_auc_mean", "brier_mean"]])
    features = new_rows.drop(columns=target)
    print(pd.DataFrame({"actual": new_rows[target].values, "predicted": model.predict(features),
                        "probability": model.predict_proba(features).max(axis=1)}))

What you get

A probability per jarReject the jars most likely to be contaminated.
Fewer false alarmsSet the threshold from the out-of-fold scores.
A strong signalThe catalogue measures a high signal.

Try it

Use Open with to run the research datasets notebook in Colab, Kaggle or Jupyter and search for "contaminant detection cocoa hazelnut", or run the code above after pip install "welt-client[parquet]".