search

becarioprecario / slm

R code with examples of the slm latent effect in R-INLA
3 stars 1 forks source link

katrina-slm-impacts.R #3

Open rsbivand opened 3 years ago

rsbivand commented 3 years ago 
> 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
rsbivand commented 3 years ago

katrina-slm.R and katrina-slm-neigh.R run OK, it's just katrina-slm-impacts.R that doesn't complete.

becarioprecario commented 3 years ago

Are you using the ‘update' branch? This is the one that is going to replace the current files on master.

El 4 feb 2021, a las 10:38, Roger Bivand notifications@github.com<mailto:notifications@github.com> escribió:

katrina-slm.R and katrina-slm-neigh.R run OK, it's just katrina-slm-impacts.R that doesn't complete.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://github.com/becarioprecario/slm/issues/3#issuecomment-773170667, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABYD6WVMGN5OUPTM2OTIVHLS5JTKNANCNFSM4XCPLEIA.

rsbivand commented 3 years ago

No, I usually ignore branches and never use them. I always assume that branches were just for local work, and are on the central repo for reference only. I've been working only in the main branch of my fork - how should I merge your update with my fork main, or update my fork of your update to my main and your update ????

rsbivand commented 3 years ago

I'm trying to merge manually locally. Do not use your update branch, I will make a PR against your main branch.

becarioprecario commented 3 years ago

I'm trying to merge manually locally. Do not use your update branch, I will make a PR against your main branch.

Thanks. I removed all the code related to the approximation that we had in the first version of the paper. I could try to make the merge if you are unsure about what to keep/remove.

rsbivand commented 3 years ago

I'm some way through. For boston, the use of eigenvalues assumes that n is small enough, but for house, the new impacts approach is not feasible. I'll PR to main branch only when completed, forget any other branches. I have no sarp.RData, will try to reconstruct from the spatialprobit package.

becarioprecario commented 3 years ago

OK. Thanks for checking this.

I'm some way through. For boston, the use of eigenvalues assumes that n is small enough, but for house, the new impacts approach is not feasible. I'll PR to main branch only when completed, forget any other branches. I have no sarp.RData, will try to reconstruct from the spatialprobit package.

— You are receiving this because you commented. Reply to this email directly, view it on GitHubhttps://github.com/becarioprecario/slm/issues/3#issuecomment-773340150, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABYD6WSLWOJV2EKCVYVO62DS5KUFVANCNFSM4XCPLEIA.