首先,識別V(c, E, HS, EC_50)具有3-parameter log-logistic{dose-response函數的形式是有用的,其中 EC_50是log-transformed EC50劑量, HS是希爾斜率或斜率因子, E是漸近響應下限, c是log-transformed劑量。 請注意exp^(...)不是有效的R語法;解決這個問題,我們得到以下V(c, E, HS, EC_50)的表單 V <- function(c, E, HS, EC_50) E + (1 - E) / (1 + exp(HS * (c - EC_50))) 我們現在可以為這三個參數定義一些(合理的)值 HS <- 1log_EC50 <- log(1e-3)E <- 0 然后計算(log-transformed劑量范圍內的響應)值c。我們將值存儲在data.frame中。 log_c <- seq(-12, 0, by = 0.1)df <- data.frame(log_c = log_c, response = V(log_c, E, HS, log_EC50)) 我們現在可以將response繪制為log_c的函數;我在這里使用ggplot2進行繪圖。 library(ggplot2)ggplot(df, aes(log_c, response