Questions tagged [drc]

Analysis of Dose-Response Curves.

Analysis of dose-response data is made available through a suite of flexible and versatile model fitting and after-fitting functions.

84 questions
1
vote
1 answer

Plotting the original x-axis from dose-response curve (drc package)

I have a dataset as follows. I want to plot a dose-response curve to show the IC50 value using the drc package. df <- structure(list(conc = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,…
MM1
  • 478
  • 15
1
vote
2 answers

How to find the X axis value (dose) at a desired Y axis (response) value in dose response curve plotted using R?

I am trying to plot a dose-response curve (for degradation) using R. I would like to get the interpolated X-axis value (dose concentration value) at the selected percent substrate remaining value (i.e., Y value, for example, Y = 50) as could be…
RanonKahn
  • 853
  • 10
  • 34
1
vote
0 answers

How to apply the drm function to subsets of a dataframe in R?

I want to do dose-response modelling with a large dataset. Is there an option to use e.g. the lapply or map function to apply the drm function to the individual subsets of the data? Or has this problem possibly to be solved by a loop? I successfully…
Anna
  • 11
  • 1
1
vote
1 answer

Calculate lethal dosis concentration calculation (LD50) on R

I am trying to calculate LD50 (lethal doses concentration that kills 50% of the organisms) for bees at different exposure concentrations (0, 3, 30 and 300ng) of a pesticide. I measure supervivency every 4 hours. Data: Control<-c(100, 100, 100, 96,…
1
vote
2 answers

How to get predictions from multiple fitted models

I have a data frame with calibration curves for multiple samples (signal s as a function of concentration c): cal <- data.frame(sample = c(rep("A", 8), rep("B", 8), rep("C", 8)), c_std = rep(c(0, 1, 5, 10, 25, 50, 100, 200), 3), …
empetrum
  • 37
  • 6
1
vote
1 answer

Plotting Dose Response Curve from Survival Data

I would like to make a dose response curve from the library(drc), and am stuck on how to prepare my dataset properly in order to make the plot. In particular, I’m struggling how to get my y-axis ready. I made up a dataframe (df) to help clarify what…
Andy
  • 413
  • 2
  • 15
1
vote
1 answer

ggplot for dose-response curve using drc package

I'm trying to plot dose-response curve in ggplot using drc package using below code and have two questions as follows. First: I need to include 0, 10, 100 etc and omit 4000 label on the x axis, how it can be done?. Second: Is it possible to squeeze…
Omi
  • 69
  • 8
1
vote
1 answer

In drc() package, drm fct = L.4 finds wrong intercept parameters, even though the graph is right

I have a problem with the following code. It calculates the drc curve correctly, but the ec50 wrongly, although the are closely related... x <- c(-1, -0.114074, 0.187521, 0.363612, 0.488551, 0.585461, 0.664642, 0.730782, 0.788875, 0.840106,…
EmaK
  • 154
  • 1
  • 9
1
vote
2 answers

coeftest() gives error with drc class object

I am trying to reproduce a supplementary example 2 Continuous response: one dose-response curve from the drc package. When trying to run coeftest(ryegrass.LL.4, vcov = sandwich), I get the following error: Error in UseMethod("estfun") : no…
Ana
  • 87
  • 9
1
vote
1 answer

Error in drc "Error in parmVec[3] - respl : non-numeric argument to binary operator"

I am pretty new to using R. I am currently trying to use drc in order to run a 4PL curve to analyze ELISA data. I have manually input the standard data and generated the curve, but I am having trouble returning the concentrations of my sample data…
wcarson
  • 13
  • 2
1
vote
0 answers

R drm power curve (DRC.powerCurve) second derivative

I am trying to calculate a power curve (Freundlich equation) in R using drm based on, https://www.statforbiology.com/nonlinearregression/usefulequations. Here is the code: model <- drm(numSpecies ~ Area, fct = DRC.powerCurve(), …
llsabang
  • 143
  • 1
  • 10
1
vote
0 answers

Is it possible to provide standard deviation to create dose response curve fits?

I am using the drc R package to fit four-parameter dose response curves to my data. I usually have replicates of four that I supply for one concentration. Would it be possible to provide the mean and standard deviation of these values instead? My…
jpquast
  • 333
  • 2
  • 8
1
vote
1 answer

Using purrr and modelr to make predictions from dose-response models (drm)

I have a data set that looks, in abbreviate form, like this: library(tidyverse) dat_s<-tibble( type=c(rep("A", 9), rep("B", 8), rep("C", 10)), ref=c("ref3", "ref3", "ref1", "ref2", "ref2", "ref1", "ref2", "ref2", "ref2", "ref2", …
aleoconn
  • 57
  • 8
1
vote
0 answers

Fitting dose-response models in python (R drc library equivalent)

I have been working with DRC library in R for quite a time. Now I would need a similar library to fit dose-response models in python, I've been googling with any discobering. Does anyone know some python package to do that? Many thanks! Xevi
Xevi
  • 379
  • 1
  • 2
  • 11
1
vote
0 answers

Wrong ED50 with Dose Response Model using DRM

When fitting a drm model to my dose response data the obtained ED50 is off log(ED50), log2(ED50), log10(ED50) # here is some actual data test_response <-…
Medibuntu
  • 11
  • 2