Displaying 1 result from an estimated 1 matches for "spg1".
Did you mean:
sp1
2012 Oct 11
2
model selection with spg and AIC (or, convert list to fitted model object)
...rt((mu1 + mu2 + p1 + p2)^2 -
4*(mu2*p1 + mu1*(mu2 + p2))))*t)/(2*(mu2*p1 + mu1*(mu2 + p2))*
sqrt((mu1 + mu2 + p1 + p2)^2 - 4*(mu2*p1 + mu1*(mu2 + p2))))
P3 <- 1 - 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...