topic Re: ED50 Calculations in SAS Health and Life Sciences

ED50 Calculations

CraigWW — Tue, 26 Jan 2010 00:20:46 GMT

I need to calculate the ED50 for a vasoconstriction pilot bioequivalence for corticosteroids as specified in the FDA Guidance "Guidance Topical Dermatologic Corticosteroids: In Vivo Bioequivalence" 2 June 1995. Please help me learn to do this in SAS.
Thanks

Re: ED50 Calculations

Cynthia_sas — Tue, 26 Jan 2010 17:00:04 GMT

Hi:
Example 71.1 in the PROC PROBIT documentation specifically mentions dosages and ED50:
http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/statug_probit_sect034.htm#

cynthia

Re: ED50 Calculations

CraigWW — Tue, 26 Jan 2010 17:16:58 GMT

Since the response for vasoconstriction pilot studies is AUC, not the number of subjects responding to a particular dose Proc Probit can not be used.

Re: ED50 Calculations

Doc_Duke — Fri, 29 Jan 2010 13:38:29 GMT

This article may help, it used something called an Emax model.

Clinical Pharmacology & Therapeutics 66, 346–357 (1 October 1999) | doi:10.1053/cp.1999.v66.a101209

Try a google search for
ed50 bioequivalence dermatologic

Re: ED50 Calculations

CraigWW — Fri, 29 Jan 2010 15:18:38 GMT

Thanks Doc,

An Emax model is what I am looking for.

Craig

Re: ED50 Calculations

k2k — Thu, 07 Feb 2013 05:05:21 GMT

dear craigWW.

i need the same for vasoconstriction study. if u have the sas code can u share me or guide me in this.

Re: ED50 Calculations

ssinojia — Fri, 11 Dec 2015 12:17:24 GMT

Hi All,

I want to calculate the ED50 and EMax using population model in SAS? Is there any built in procedure available in SAS to do this? I have read on some paper about SAS MIXLIN Macro, anyone have idea about this Macro?

Re: ED50 Calculations

Irla — Mon, 22 Aug 2022 14:12:22 GMT

Hi!

Have you found any code in SAS to calculate Emax with vasoconstriction data? I'm currently with the same issue. Thank you in advance

Re: ED50 Calculations

ssinojia — Mon, 22 Aug 2022 18:26:30 GMT

Hi Irla,

I didn’t find exact code as such however I was able to calculate ED50 using PROC NLMIXED by fitting the model. I don’t have that handy but you can try searching ED50 or Emax modeling with PROC NLMIXED. I am sure you will get some good results with Google.

Re: ED50 Calculations

Irla — Mon, 22 Aug 2022 18:32:44 GMT

Hey, thank you for answering me. I'll try to find it right now!

Re: ED50 Calculations

Irla — Tue, 23 Aug 2022 13:42:45 GMT

Hey! Since you told me you got to calculate it, can I ask you something? I've found a NLIN code, probably I'll be able to reproduce it in SAS.

But there some parameters to put on the code where I'm stuck. We use DOSE DURATION X AUEC, right? In the formula we have to type the inicial values, right? How do I find them? What about the slope, how to obtain this data?

that is the code:

proc nlin data = dr;
  parms top = 68.875 bottom = 20.575 EC50 = 1 hill = 1;
  model response = top + (bottom - top) / 
  (1 + (concentration / EC50)**hill);
run;

Re: ED50 Calculations

ssinojia — Tue, 23 Aug 2022 17:39:30 GMT

Hi,

I don’t remember them now, but I used to get that from my statistician. Then used to see model fitting. Again I used NLMIXED.

Re: ED50 Calculations

SteveDenham — Thu, 05 Jan 2023 13:32:40 GMT

First off, you may want to check the parameterization of the model that you are using. The example at the link below uses a more "interpretable" parameterization that you might find also has better properties for your goal.

If you don't have a good idea from prior knowledge as to what the best starting value(s) might be, you might try a grid search over a range of reasonable values. You don't need the exact starting values this way, and if several combinations of starting values lead to the same optimum solution, you get what you are looking for without knowing a precise set of values. This also can help you avoid converging to a saddle point rather than an extremum. See the documentation for the PARAMETER statement in NLIN https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/statug/statug_nlin_syntax12.htm and this example which goes into various parameterizations for an ED50 model https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/statug/statug_nlin_examples04.htm

Good luck.

SteveDenham

Re: ED50 Calculations

ssinojia — Thu, 05 Jan 2023 14:07:04 GMT

Hi Steve, thanks for your inputs. It’s helpful.