got_evd <- requireNamespace("evd", quietly = TRUE)
if (got_evd) {
library(evd)
small <- fgev(ow$temp, nsloc = ow[, "loc"])
adj_small <- alogLik(small, cluster = ow$year)
tiny <- fgev(ow$temp)
adj_tiny <- alogLik(tiny, cluster = ow$year)
anova(adj_small, adj_tiny)
set.seed(4082019)
uvdata <- evd::rgev(100, loc = 0.13, scale = 1.1, shape = 0.2)
M0 <- fgev(uvdata)
M1 <- fgev(uvdata, nsloc = (-49:50)/100)
adj0 <- alogLik(M0)
adj1 <- alogLik(M1)
anova(adj1, adj0)
}
got_extRemes <- requireNamespace("extRemes", quietly = TRUE)
if (got_extRemes) {
library(extRemes)
large <- fevd(temp, ow, location.fun = ~ loc, scale.fun = ~ loc,
shape.fun = ~ loc)
medium <- fevd(temp, ow, location.fun = ~ loc, scale.fun = ~ loc)
small <- fevd(temp, ow, location.fun = ~ loc)
tiny <- fevd(temp, ow)
adj_large <- alogLik(large, cluster = ow$year)
adj_medium <- alogLik(medium, cluster = ow$year)
adj_small <- alogLik(small, cluster = ow$year)
adj_tiny <- alogLik(tiny, cluster = ow$year)
anova(adj_large, adj_medium, adj_small, adj_tiny)
}
Run the code above in your browser using DataLab