I am trying to illustrate the three types of extreme value distributions in R.
For type one, the distribution function is as follows.
The problem is that the || operator only looks at the first element of z<b and z==b. This should work:
||
z<b
z==b
G2 <- function(z, a, b, alpha) ifelse(z<=b, 0,exp(-((z-b)/a)^(-alpha)))