How do I create a Q-Q plot with Gamma distribution?

Question

I have continuous data that I determined had a gamma distribution. The key data are here:

data <- data.frame(yr = c("2008", "2009", "2