#Get some CDs, form an ungauged pooling group and apply EVPlot.
CDs.28015 <- GetCDs(28015)
Pool.28015 <- Pool(CDs.28015, exclude = 28015)
EVPool(Pool.28015)
#Do the same for the gauged case, change the title, and convert with a QMED of 105.5.
PoolG.28015 <- Pool(CDs.28015)
EVPool(PoolG.28015, gauged = TRUE, Title = "Gauged frequency curve - Site 28015", QMED = 9.8)
#Pretend we have an extra AMAX for the gauge. Amend the pooling group Lcv and LSkew
#for the site accordingly then apply EVPool with the updated AMAX.
#Firstly, get the AMAX sample
AM.28015 <- GetAM(28015)
#Add an extra AMAX flow of 15m3/s
Append28015 <- append(AM.28015$Flow, 15)
#Amend the Lcv and Lskew in the pooling group
PoolG.28015[1, c(16, 17)] <- c(Lcv(Append28015), LSkew(Append28015))
#Now plot gauged with the updated AMAX
EVPool(PoolG.28015, AMAX = Append28015, gauged = TRUE)
Run the code above in your browser using DataLab