mrlplot(esk.rain$precip)

tcplot(esk.rain$precip, tlim=c(0.1,7), model="gpd", nt = 20)

### here we take a fixed threshold u=5, but you can choose u based on part (a)
thresh.fit_gpd<-fpot(esk.rain$precip, threshold=5, start=list(scale=1.2,shape=0.1))
### next, we look at the resulting diagnostic plots
par(mfrow=c(1,2))
plot(thresh.fit_gpd)