data(voting)
p1 <- Lindsey(voting,voting_tally)
p2 <- optimizer(voting,voting_tally,start=p1)
theta(p1) - theta(p2) # Should be zero
if (FALSE) {
data(pollen)
p1 <- optimizer(pollen)
p2 <- Lindsey(pollen)
theta(p1) - theta(p2) # Isn't zero...numerical scruff...
}
Run the code above in your browser using DataLab