Displaying 1 result from an estimated 1 matches for "spg2".
Did you mean:
sp2
2012 Oct 11
2
model selection with spg and AIC (or, convert list to fitted model object)
...P1 - P2
p.all <- c(P1, P2, P3)
-sum(dmnom2(c(n1, n2, n3), prob = p.all, log = TRUE))
}
# Optimization with spg
spg1 <- spg(par = c(2, 0.3, 0.001, 0.001), fn = NLL1, y = y,
lower = c(0.0001, 0.0001, 0.0001, 0.0001),
control = list(maxit = 5000, trace = FALSE))
spg2 <- spg(par = c(2, 0.001), fn = NLL2, y = y,
lower = c(0.0001, 0.0001),
control = list(maxit = 5000, trace = FALSE))
# Model selection with ICtab; currently gives an error
ICtab(spg1, spg2)
##########################################################################
Any...