[R] How to fit a Hill model using braidrm package?
Luigi Marongiu
m@rong|u@|u|g| @end|ng |rom gm@||@com
Thu Sep 25 08:51:39 CEST 2025
Thank you but I am getting several error installing drc on Fedora...
On Wed, Sep 24, 2025 at 1:47 PM peter dalgaard <pdalgd using gmail.com> wrote:
>
> They're in basicdrm, not braidrm...
>
> -pd
>
> > On 24 Sep 2025, at 11:52 , Luigi Marongiu <marongiu.luigi using gmail.com> wrote:
> >
> > Hello,
> > I have a set of data coming from a dissociation experiment
> > (protein/ligand). Since the data is required to calculate the constant
> > of dissociation (Kd) of this pair, I am looking for a way of fitting a
> > Hill function to the data.
> > I have seen that the package braidrm
> > (https://cran.r-project.org/web/packages/braidrm/index.html) provides
> > this function, but when I launch the function `evalHillEqn`,
> > `findBestHill` and so forth I get the error of function not found.
> > Yet, the package is given as properly installed by the system.
> > How can I run this package?
> > Is there an alternative way to fit a Hill function to these data?
> > Thank you
> >
> > ```
> > df = data.frame(Response = c(890.72, 895.46, 895.63,
> > 894.72, 895.49, 893.59,
> > 892.53, 895.06, 897.21, 889.27, 876.05,
> > 857.96, 862.02, 858.36,
> > 890.94, 890.8, 887.22, 888.91, 890.83,
> > 889.92, 891.76, 890.32,
> > 886.35, 878.11, 866.57, 859.04, 863.64,
> > 880.16, 884.15, 879.57,
> > 878.89, 882.27, 881.59, 880.98, 881.45,
> > 876.19, 868.32, 859.16,
> > 850.53, 853.21, 859.34, 859.73, 861.19),
> > Dose = c(0.0000000015, 0.000000003,
> > 0.000000006, 0.000000012,
> > 0.000000024, 0.000000048, 0.000000095,
> > 0.00000018,
> > 0.00000038, 0.00000078, 0.0000015,
> > 0.000013, 0.000025,
> > 0.00005, 0.0000000015, 0.000000003,
> > 0.000000006,
> > 0.000000012, 0.000000024, 0.000000048,
> > 0.000000095,
> > 0.00000018, 0.00000038, 0.00000078,
> > 0.0000015, 0.000025,
> > 0.00005, 0.0000000015, 0.000000003,
> > 0.000000006,
> > 0.000000012, 0.000000024, 0.000000048,
> > 0.000000095,
> > 0.00000018, 0.00000038, 0.00000078,
> > 0.0000015, 0.000003,
> > 0.000006, 0.000013, 0.000025, 0.00005)
> > )
> > plot(Response~log10(Dose), df)
> > library(braidrm)
> > evalHillEqn(Response~Dose, df)
> > findBestHill(Response~Dose, df)
> > ```
> >
> > ______________________________________________
> > R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide https://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
>
> --
> Peter Dalgaard, Professor,
> Center for Statistics, Copenhagen Business School
> Solbjerg Plads 3, 2000 Frederiksberg, Denmark
> Phone: (+45)38153501
> Office: A 4.23
> Email: pd.mes using cbs.dk Priv: PDalgd using gmail.com
>
--
Best regards,
Luigi
More information about the R-help
mailing list