View License

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video

Highlights from
Trimmed Spearman-Karber method.

5.0 | 1 rating Rate this file 8 Downloads (last 30 days) File Size: 6.52 KB File ID: #28479 Version: 1.5
image thumbnail

Trimmed Spearman-Karber method.



13 Aug 2010 (Updated )

This function analyzes a dose-response curve using the Trimmed Spearman-Karber method.

| Watch this File

File Information

The MATLAB function in TSK.m analyzes a dose-response curve and calculates the median dose (for instance, the lethal dose for 50% of the population (LD50)) using the Trimmed Spearman-Karber method. This is an alternative to dose-response analysis methods such as probit and logit methods. It does not require glmfit or the stats toolbox.

All calculations in this code are based on the method from Hamilton, M.A.; Russo, R.C.; Thurston, R.V. Trimmed Spearman-Karber Method for Estimating Median Lethal Concentrations in Toxicity Bioassays. Enviro. Sci. Tech. 1977, 11 (7), 714-719: See also ibid, 1978, 12 (4), 417:

This code was validated using the data given in Hamilton et al. (1977) and the example given in the documentation for the EPA's DOS program to perform the Spearman-Karber method. (Note that the confidence intervals given in Hamilton et al. (1977) are incorrect, as per Hamilton et al. (1978).) This code replicates the EPA program's results. Hamilton's test cases are provided in TSKTest.m.


Herrorbar inspired this file.

MATLAB release MATLAB 7.5 (R2007b)
Tags for This File   Please login to tag files.
Please login to add a comment or rating.
Comments and Ratings (2)
08 Nov 2016 mopper

mopper (view profile)


28 Sep 2010 Brenton Stone

Someone asked me via e-mail if it would be better to use the Student's distribution instead of the normal distribution to calculate the confidence intervals. This isn't what was done in Hamilton, but if you want to modify the code and have an install of the stats toolbox, change line 248 to:

v = tinv((1-conf)/2,s+1);

Note: I don't have a copy of the stats toolbox and have not tested this.

Comment only
13 Aug 2010 1.3

text editing, adding URLs: no changes to code

16 Aug 2010 1.4

adding BSD organization clause

25 Mar 2011 1.5

bug fixes

Contact us