# NOT RUN {
#Get an AMAX sample and derive GoTF score against the GenLog and the GEV distributions
# }
# NOT RUN {
AM <- GetAM(203018)
# }
# NOT RUN {
GoTF(AM$Flow, dist = "GenLog")
# }
# NOT RUN {
GoTF(AM$Flow, dist = "GEV")
# }
# NOT RUN {
#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
# }
# NOT RUN {
GoTF(AM$Flow, GF = c(0.16, 0.2, median(AM$Flow)))
# }
# NOT RUN {
#calculate the GoTF based on parameters of the GenLog estimated inadequately.
# }
# NOT RUN {
Loc <- mean(AM$Flow)
# }
# NOT RUN {
Scale <- sd(AM$Flow)
# }
# NOT RUN {
Skew <- 1-(median(AM$Flow)/mean(AM$Flow))
# }
# NOT RUN {
GoTF(AM$Flow, pars = c(Loc, Scale, Skew))
# }
Run the code above in your browser using DataLab