R - рассчитать вероятность и перевернуть ось x/y кумулятивной кривой (ECDF)
В R I построите кумулятивную кривую, используя функцию ecdf для отображения площади и высоты. По умолчанию высота отображается на оси x, в области на оси y, где высота указывается в общих значениях (например, 1000–3000 м), и в области вероятности (0,0–1,0).
Я хотел бы показать повышение в процентах / значения вероятности тоже.
Я хотел бы перевернуть ось и визуализировать область на оси х и высоту на оси у.
Как я могу это сделать? Я не мог найти решение.
Это мой код:
library(rgdal)
setwd("G:/data")
raster <- readGDAL("elevation.tif")
plot(ecdf(raster$band1), main="Hypsometric curve")
1 ответ
Практически любое решение начнется с создания массива, представляющего все вершины в графике. Они получены путем извлечения координат из ecdf
объект. Для иллюстрации давайте сгенерируем некоторые данные.
n <- 32
k <- floor(n/2)
x <- 1 - runif(n)^2
x <- c(x, sample(x, k, replace=TRUE))
Вот стандартный расчет ECDF и график, для справки:
f <- ecdf(x)
plot(f, main="Default Plot")
R
предпочитает помещать данные в data.frame
объекты. Вот один из способов. Обратите внимание, как диапазон данных немного расширяется, чтобы показать естественные конечные точки кривой, которая начинается на уровне y=0 и заканчивается на уровне y=1:
expand <- function(x, f=1.05, ...) {# Expand a data range
r <- range(x, ...)
(r - mean(r)) * f + mean(r)
}
X <- with(environment(f), {
x.range <- expand(x)
data.frame(x=c(x.range[1], x),
xend=c(x, x.range[2]),
y=c(0, y),
yend=c(0, y))
})
С этой удобной структурой данных X
на месте, все что осталось, это описать сюжет. я использую ggplot2
за это:
library(ggplot2)
g.ECDF <- ggplot(X, aes(x=x, y=y, xend=xend, yend=yend)) +
geom_hline(yintercept=0:1, linetype=2, size=1.0, color="Gray") +
geom_segment(color="#b0b0b0", size=0.8) +
geom_point(data=X[-1,]) +
xlab("Elevation") + ylab("Proportion of Total") +
ggtitle("Hypsometric Curve")
print(g.ECDF)
Теперь для государственного переворота: просто скажи ggplot2
для переключения осей:
print(g.ECDF + coord_flip() + ggtitle("Axes Flipped"))