The slm from becarioprecario

slm's People

Contributors

Stargazers

Watchers

slm's Issues

katrina-slm-impacts.R

> tabdi$INLASDM<-(inla.zmarginal(trIrhoWbsdm, TRUE)$mean * sdmm1$summary.random$idx[idxbeta,"mean"]+inla.zmarginal(trIrhoWgsdm, TRUE)$mean * sdmm1$sum .... [TRUNCATED] 
Error in splinefun(xnew, xx, method = "hyman") : zero non-NA points
In addition: Warning message:
In sqrt(dx) : NaNs produced
> traceback()
5: stop("zero non-NA points")
4: splinefun(xnew, xx, method = "hyman")
3: inla.smarginal(marginal)
2: inla.emarginal(function(xx) c(xx, xx^2), marginal)
1: inla.zmarginal(trIrhoWbsdm, TRUE)

Did something change in the internals?

> o <- inla.zmarginal(trIrhoWbsdm, TRUE)$mean
Error in splinefun(xnew, xx, method = "hyman") : zero non-NA points
In addition: Warning message:
In sqrt(dx) : NaNs produced
> o <- inla.zmarginal(trIrhoWgsdm, TRUE)$mean
Error in splinefun(xnew, xx, method = "hyman") : zero non-NA points
In addition: Warning message:
In sqrt(dx) : NaNs produced

First INLA fit in katrina-slm.R fails

Any ideas?

> semm1<-inla(
+    update(f1, ~.+f(idx, model="slm", args.slm=args.slm, hyper=hyper.slm)),
+    data=Katrina, family="binomial",
+    control.family = list(link="probit", hyper=zero.variance),
+    control.compute=list(dic=TRUE, cpo=TRUE)
+ )
Error in inla.set.hyper(family[i.family], "likelihood", cont.family[[i.family]]$hyper,  : 
  Model binomial [ likelihood ], has none hyperparameters, but 'hyper' is  list(prec = list(initial = 25, fixed = TRUE))

Recommend Projects