我在 R 中有我的代碼(模型):
library(ggplot2)
V <- function(c, E, HS, EC_50) E / (1 exp(HS * (c - EC_50)))
HS <- 1
log_EC50 <- log(1e-3)
E <- 0
log_c <- seq(-12, 0, by = 0.1)
df <- data.frame(log_c = log_c, response = V(log_c, E, HS, log_EC50))
ggplot(df, aes( response)) geom_histogram(binwidth=0.03)
我想創建一個函式,該函式將像在我的代碼中一樣繪制直方圖,但會考慮引數的變化,例如 HS、E。
請幫忙
uj5u.com熱心網友回復:
只需包裝ggplot到函式中,例如:
library(ggplot2)
log_c <- seq(-12, 0, by = 0.1)
E <- 1
HS <- 1
log_EC50 <- log(1e-3)
V <- function(c, E, HS, EC_50) {
E / (1 exp(HS * (c - EC_50)))
}
a <- function(log_c, E, HS, log_EC_50) {
df <- data.frame(log_c = {{log_c}}, response = V({{log_c}}, {{E}}, {{HS}}, {{log_EC50}}))
ggplot2::ggplot(df, aes( response)) geom_histogram(binwidth=0.03)
}
a(log_c, E=1, HS=2, log_EC50)

a(log_c, E=1, HS=5, log_EC50)

使用reprex v2.0.2創建于 2022-10-29
轉載請註明出處,本文鏈接:https://www.uj5u.com/net/523715.html
標籤:r功能阴谋
上一篇:C函式回傳字串長度
