#Get an AMAX sample and derive GoTF score against the GenLog and the GEV distributions
AM <- GetAM(203018)
GoTF(AM$Flow, dist = "GenLog")
GoTF(AM$Flow, dist = "GEV")
#Derive the GF parameters for the ungauged pooled estimate for the AM and
#calculate a GoTF for GenLog (assuming the gauged QMED)
#For this assume 0.16 and 0.2 as the ungauged Lcv & LSkew pooled estimates
GoTF(AM$Flow, GF = c(0.16, 0.2, median(AM$Flow)))
#calculate the GoTF based on parameters of the GenLog estimated inadequately.
Loc <- mean(AM$Flow)
Scale <- sd(AM$Flow)
Skew <- 1-(median(AM$Flow)/mean(AM$Flow))
GoTF(AM$Flow, pars = c(Loc, Scale, Skew))
Run the code above in your browser using DataLab