In reply to Rick Wicklin.

Thank you for your response. The needed link function of GEV has the following form: {{[-ln(cdf of GEV)]^(-xi)}-1}/xi=(transpose(beta))(x_i). How can it be implemented in generalized linear model?