看看那个函数(ismev包)的定义:
<br />
gev.diag<br />
function (z)<br />
{<br />
n <- length(z$data)<br />
x <- (1:n)/(n + 1)<br />
if (z$trans) {<br />
oldpar <- par(mfrow = c(1, 2))<br />
plot(x, exp(-exp(-sort(z$data))), xlab = "Empirical",<br />
ylab = "Model")<br />
abline(0, 1, col = 4)<br />
title("Residual Probability Plot")<br />
plot(-log(-log(x)), sort(z$data), ylab = "Empirical",<br />
xlab = "Model")<br />
abline(0, 1, col = 4)<br />
title("Residual Quantile Plot (Gumbel Scale)")<br />
}<br />
else {<br />
oldpar <- par(mfrow = c(2, 2))<br />
gev.pp(z$mle, z$data)<br />
gev.qq(z$mle, z$data)<br />
gev.rl(z$mle, z$cov, z$data)<br />
gev.his(z$mle, z$data)<br />
}<br />
par(oldpar)<br />
invisible()<br />
}</p>
<p>
</p>
改对应的xlab,ylab就好了嘛。