Displaying 1 result from an estimated 1 matches for "parlpe3".
Did you mean:
carpe3
2010 Feb 22
1
lmom: plotting log Pearson Type III
...29900, 28000, 25600,
19700, 25700, 29500, 26800, 30000, 29500)
# estimate moments
moments = samlmu(mackenzieRiver, sort.data = TRUE)
log.moments <- samlmu( log(mackenzieRiver), sort.data = TRUE )
# estimate parameters
parGEV <- pelgev(moments) # GEV
parPE3 <- pelpe3(moments) # Pearson
parLPE3 <- pelpe3(log.moments) # log Pearson
# plot result
evplot(mackenzieRiver, rp.axis = TRUE)
evdistq(quagev, parGEV, col = 'black')
evdistq(quape3, parPE3, col = 'blue')
# estimate 1:100 yr event
flood.est <- list(
GEV = quagev(0.99, parGEV),...