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
2 answers

mselect function of drc package does not work inside a function

I wanted to use mselec function from drc package for selecting the best model of a dose-response dataset. However, mselect does not work if you use it inside a function. The following code works: library(drc) ryegrass.m1 <- drm(rootl~conc, data =…
1
vote
1 answer

Why single confidence interval when plotting drm model with multiple curves?

I'm doing nonlinear regression in R using drm function from drc package. The drm function takes a curveid argument, which makes drm to fit multiple curves and store the results in single model object. Next, I plot the curves and confidence…
Vallo Varik
  • 137
  • 1
  • 10
1
vote
1 answer

Fitting steep dose response curves

Hi I'm trying to fit a rather steep dose-response curve in the drc package in R. Moles is my x function and pct_control is my y function. I tried fitting it on the Log-logistic 3 parameter function and the fit was terrible. I tried other model fits…
VChu
  • 11
  • 4
1
vote
2 answers

Plotting drc model in ggplot2; issue with seq( )

My model does not continue towards the asymptotes when plotted in ggplot2, though it does in R base graphics. In ggplot2, it stops at certain points on the X axis (images included below), I am 90% certain this is related to seq(). I am using…
Arch
  • 192
  • 2
  • 16
1
vote
1 answer

drc: Error in drm when used interaction terms

I want to fit log-logistic regression for the following data in drc R package for the combination of Temp and Variety. However, my code throws the following error Error in Temp:Variety : NA/NaN argument Code: df2 <- structure(list(Temp = c(15L,…
MYaseen208
  • 22,666
  • 37
  • 165
  • 309
1
vote
0 answers

Fitting Multiple Logistic Regression with Interaction between Quantitative and Qualitative Explanatory Variables with drm function from drc package

As a follow up to this question answered by @EDi. I wonder how to fit the following glm model with drc function from drc package. GLM Code Type <- rep(x=LETTERS[1:3], each=5, times=2) Conc <- rep(rep(x=seq(from=0, to=40, by=10), times=3), 2) Rep …
MYaseen208
  • 22,666
  • 37
  • 165
  • 309
0
votes
0 answers

Dose-response curve, geom_segment

When use the function of scale_x_continuous, geom_segment of y axis dose not working. (x=coefs["e"], y=y50, xend=0, yend=y50) Wheon use the function of scale_x_log10, it dose work but x axis labels are very weird. So I want to use…
0
votes
0 answers

Dose-response curve, geom_point

I plotted the dose-response curve with the code below. However, the points on the graph don't seem to be correct. I would like to know how to plot the y-values corresponding to the x-values with the correct…
0
votes
1 answer

Dose-response curve, ED50

I used the following code to plot the dose-response curve and display the ED50 value. However, the ED50 value was displayed at a value other than 50 on the y-axis. I used the ED function to find the ED50 value, but why am I getting these results? Is…
0
votes
1 answer

How to store model coefficients to be used as variables in a function in R Shiny?

I'm creating an app that takes a .csv file and spits out concentration values after doing some analysis on the back end. The analysis involves running a 5-parameter logistic curve using the drc() package and then storing the model coefficients to…
0
votes
1 answer

Non finite value error when trying to use drm() to fit logistic regression model to the data

I have the dataframe below d50<-structure(list(VOT = c("VOT 0", "VOT 5", "VOT 10", "VOT 15", "VOT 20", "VOT -10", "VOT -20", "VOT -30", "VOT -40", "VOT -50", "VOT -60", "VOT 0", "VOT 5", "VOT 10", "VOT 15", "VOT 20", "VOT -10", "VOT -20", "VOT…
firmo23
  • 7,490
  • 2
  • 38
  • 114
0
votes
0 answers

Why does datafreme error after using GRfit for grouping?

I saw error in R-script. library(GRmetrics) df <- read.csv("E:/Result 2023-1-19/KKU100.csv") head(mydata) df1 <- na.omit(df) drc_output = GRfit(df1, groupingVariables = c('cell_line','treatment')) drc_output and the error is : Error in…
0
votes
1 answer

R drc Plot y-axis scale issues

Trying to make a simple plot and there is a large gap between the first and second data points (0 and 10) hatch <- drm(ha_.~leachate, data=data, fct = LL.3()) plot(hatch, xlim=c(0,500),ylim=c(0,100), xlab="Leachate Concentration…
0
votes
0 answers

DRC multiple curves plotting issue

So I saw some question that are similar to my problem, the thing is that the answers are not working for me. drc_ll.4_ppo <- drm(DWP~ai, pop, data = DRC_PPO, fct = LL.4(fixed = c(NA, NA, 100, NA))) upperl = c(NA, NA, 100, NA) lowerl = c(NA, 0, NA,…
0
votes
0 answers

'ungroup error' when using summary table for a drc model

I have a drc model in R, exponential decay with 3 parameters. m1<- drm(df$Time ~ df$DeltaYear, fct=EXD.3()) modelFit(m1) summary(m1) summary(m1) gives me the model results but I need to print a table of this in R markdown. I have used the…
bee761
  • 25
  • 3