Install any skill in seconds. Free to start, no credit card required.
Get Started Free →mims-harvard/tooluniverse-dose-response
Dose-response / concentration-response curve fitting — IC50, EC50, Hill slope, Emax/Emin efficacy, and relative potency from paired concentration vs response data (enzyme/cell assays, drug screening, agonist/antagonist pharmacology). Fits the 4-parameter logistic (Hill sigmoidal) model.
Passed 18 of 22 test cases with the skill, and 13 of 22 without it.
and 2 more in the full report
5 fixed13 passed either way4 still failing0 broke
Tokens, turns and time describe resource use, not a dollar price.
One of 5 test cases the skill fixed. Same model, same prompt, with and without the skill.
I want to extract only the full 4-parameter logistic curve parameters for a single dataset without calling the summary helper wrapper. Which ToolUniverse function provides this?
| Checked | Without skill | With skill |
|---|---|---|
| Identifies DoseResponse_fit_curve as the tool function to call. | Missed | Passed |
The instructions your model reads when the skill loads.
Turn paired concentration vs response measurements into a potency (IC50/EC50), a Hill slope, an efficacy (Emax), and a quality judgment — and compare potency between compounds.
The model is the 4-parameter logistic (4PL) / Hill sigmoidal:
f(x) = Emin + (Emax − Emin) / (1 + (EC50/x)^n) — where n is the Hill slope.
| Issue | What to do |
|---|---|
| Concentration units | Pick ONE unit (µM, nM, M) and use it for every point. The IC50 comes back in that unit. Don't mix. |
| Log vs linear concentrations | Pass concentrations on the linear scale (e.g. 0.01, 0.1, 1, 10), not log10. The fitter logs internally. |
| Zero/control concentration | Drop a literal 0 concentration (log(0) is undefined). Keep it only as the Emin/Emax reference if normalizing. |
| Direction | Inhibition curves go high→low (IC50); activation curves go low→high (EC50). The tools handle both; just be consistent. |
| Normalization | Convert raw signal to % of control if you want Emax/Emin near 100/0: % = 100 × (raw − blank)/(control − blank). Raw values also fit, but plateaus are then in raw units. |
| Replicates | Average technical replicates per concentration before fitting, or pass all points (the fit weights them equally). |
Coverage requirement: you need ≥4 points (the tools require it) and ideally 6–8 spanning both plateaus — points clearly above and clearly below the inflection. A curve that never plateaus gives an unreliable, extrapolated IC50 (see Step 4).
Single curve → IC50 or EC50 (same math; "IC50" for inhibition, "EC50" for activation):
tu run DoseResponse_calculate_ic50 '{"operation":"calculate_ic50",
"concentrations":[0.001,0.01,0.1,1,10,100],
"responses":[98,95,80,45,12,3]}'
Returns ic50, ic50_95_confidence_interval, hill_slope, emax, emin, r_squared, log_ic50.
Full 4PL parameters only → DoseResponse_fit_curve (same inputs). Two compounds → DoseResponse_compare_potency with conc_a/resp_a/conc_b/resp_b (returns each IC50 + ic50_fold_shift_b_over_a + more_potent).
test cases passed, vs 13 / 22 without
Published skills whose descriptions closely match this one, with their headline benchmark lift.
More skills from the same source and nearby in the registry, with their headline benchmark lift.