我有代碼。
V <- function(C, E, HS, EC_50) {
response <- E (1 - E) / (1 exp(HS * (C - EC_50)))
list(response = response,
mean_to_hist = mean(response),
sd_to_hist = 1)
}
example1 <- V(seq(-0.1, 0, by = 0.01),0,1,log(1e-3))
example2 <-V(seq(-2500, 0, by = 0.01),0,1,log(1e-3))
Fluorescence_Intensity <- function(my_V, METHOD, ...){
fun <- METHOD
fun(mean = my_V$mean_to_hist, sd = my_V$sd_to_hist, ...)
}
x <- Fluorescence_Intensity(example1, rnorm, n = 1000)
y <- Fluorescence_Intensity(example2, rnorm, n = 1000)
x.bar <- mean(x)
y.bar <- mean(y)
x.q95 <- quantile(x, 0.95)
hist(x, xlim = c(0, 40), col = scales::alpha('gray70', 0.4), border = FALSE)
hist(y, add = TRUE, col = scales::alpha('gray70'), border = FALSE)
abline(v = c(x.bar, y.bar, x.q95), col = c("green", "green", "blue"), lwd = 2)
hx <- hist(x, plot = FALSE)
hy <- hist(y, plot = FALSE)
xlim <- range(c(x, y))
ylim <- c(0, max(hx$counts, hy$counts))
plot(1, type = "n", xlim = xlim, ylim = ylim,
xlab = "Fluorescence Intensity", ylab = "Frequency")
grid()
rect(xlim[1] - 10, -10, x.q95, ylim[2] 10, col = scales::alpha('lightblue', 0.4), border = FALSE)
hist(x, add = TRUE, col = scales::alpha('gray70', 0.4), border = FALSE)
hist(y, add = TRUE, col = 'gray70', border = FALSE)
abline(v = c(x.bar, y.bar, x.q95), col = c("green", "green", "blue"), lwd = 2)
我想更改它，以便只有較暗的直方圖和較綠的均值條可見，但藍色條和背景位于完全相同的位置。（只要不可見，就可以在代碼中創建光直方圖）。
uj5u.com熱心網友回復：
它解決了你的問題嗎？（移除光直方圖及其相關的均值作為斜線）
example1 <- V(seq(-0.1, 0, by = 0.01),0,1,log(1e-3))
example2 <-V(seq(-2500, 0, by = 0.01),0,1,log(1e-3))
Fluorescence_Intensity <- function(my_V, METHOD, ...){
fun <- METHOD
fun(mean = my_V$mean_to_hist, sd = my_V$sd_to_hist, ...)
}
x <- Fluorescence_Intensity(example1, rnorm, n = 1000)
y <- Fluorescence_Intensity(example2, rnorm, n = 1000)
x.bar <- mean(x)
y.bar <- mean(y)
x.q95 <- quantile(x, 0.95)
hist(x, xlim = c(0, 40), col = scales::alpha('gray70', 0.4), border = FALSE)
hist(y, add = TRUE, col = scales::alpha('gray70'), border = FALSE)
abline(v = c(x.bar, y.bar, x.q95), col = c("green", "green", "blue"), lwd = 2)
hx <- hist(x, plot = FALSE)
hy <- hist(y, plot = FALSE)
xlim <- range(c(x, y))
ylim <- c(0, max(hx$counts, hy$counts))
plot(1, type = "n", xlim = xlim, ylim = ylim,
xlab = "Fluorescence Intensity", ylab = "Frequency")
grid()
rect(xlim[1] - 10, -10, x.q95, ylim[2] 10, col = scales::alpha('lightblue', 0.4), border = FALSE)
hist(y, add = TRUE, col = 'gray70', border = FALSE)
abline(v = c(y.bar, x.q95), col = c("green", "blue"), lwd = 2)
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/527689.html
標籤：r阴谋