Hi y'all!
I 'm doing a special Q-Q plot and I'm trying to control the xlab, but it doesnt work.
Here is my code
monthly$quantile=exp(-(1-0.808*((monthly$MonthlyMax-82.119)/8.085))^(-1/-0.808))
monthly$returnperiod=1/(1-monthly$quantile)
model=vector("numeric",length(sort(monthly$MonthlyMax)))
for (i in 1:length(model)){
model=qgev(cdf,-0.808,82.119,8.085,lower.tail=TRUE)
}
plot(model~sort(monthly$MonthlyMax),xlab='Return Periods',main='Q-Q plot',xaxt='n')
Do you have any idea why?
MonthlyMax MonthlyMin Date
1: 70.4 55.2 2015-01-01
2: 74.7 56.3 2015-02-01
3: 77.5 61.8 2015-03-01
4: 82.9 63.1 2015-04-01
5: 88.5 77.8 2015-05-01
6: 89.5 83.9 2015-06-01
7: 89.5 81.5 2015-07-01
8: 91.9 81.0 2015-08-01
9: 87.6 79.7 2015-09-01
10: 85.1 69.6 2015-10-01
11: 81.9 64.7 2015-11-01
12: 75.6 -99.0 2015-12-01
13: 72.0 40.2 2016-01-01
14: 75.5 48.9 2016-02-01
15: 76.2 -99.0 2016-03-01
16: 82.9 71.0 2016-04-01
17: 87.8 -99.0 2016-05-01
18: 89.9 80.2 2016-06-01
19: 91.7 83.0 2016-07-01
20: 89.1 80.6 2016-08-01