andybega / spduration

Split-population duration model in R
https://andybeger.com/spduration/
4 stars 0 forks source link

spdur loglog version #1

Closed andybega closed 11 years ago

andybega commented 11 years ago

library(spduration)

data(insurgency) duration.ins <- buildDuration(insurgency, 'insurgency', unitID='ccode', tID='date')

model_loglog <- spdur( duration ~ low_intensity + high_neighbors + exclpop.l1, atrisk ~ excl_groups_count.l1 + high_neighborhood + high_intensity + exclpop.l1 + lgdppc.l1, last='end.spell', data=duration.ins, distr="loglog", max.iter=300)

Fitting base loglog... Error in optim(base.inits, loglog.lik, method = "BFGS", control = list(maxit = max.iter), : initial value in 'vmmin' is not finite