我目前正在研究超出极端降雨概率的图表。但是,我正在努力解决这些问题。
这是我到目前为止所得到的:
PE2010 <- read.table(file="C:/Users/PE2001-2010.txt",header=TRUE)
plot(PE2010$Pre ~ PE2010$PE,
col = "red",
xlab = "Overschrijdingskans",
ylab = "Neerslagintensiteit (mm/d)",
main = "Overschrijdingskans voor de 3 simulaties")
结果如图所示:
我想做什么:
总结一下,我希望结果是这样的:
答案 0 :(得分:0)
如果没有您的实际数据,很难让它完全正确,但这可能是您想要的很好的近似值:
# an index that will put the data in order from left to right
indx <- order(PE2010$PE)
plot(y =PE2010$Pre[indx],
x = PE2010$PE[indx],
col = "red",
xlab = "Overschrijdingskans",
ylab = "Neerslagintensiteit (mm/d)",
main = "Overschrijdingskans voor de 3 simulaties",
type = "l",
# log transformt the x-axis
log = 'x',
# reverse the x axis
xlim = range(PE2010$PE,na.rm=T)[2:1],
# don't include the default axes
axes=F
)
# make the axes
box()
axis(2)
axis(1,at=c(1, 0.1, 0.01,0.001 ))