Displaying 2 results from an estimated 2 matches for "casedata".
2012 Oct 30
2
issues with krige function
...strange problem with the krige function from geoR. The
problem that I am having is that while the krige function seems to
work well, the resulting predicted values are all NAs. Given the size
of the datasets I am working with can't attach it, but I can provide
snippets of the datasets.
> casedata
station year month day obs mpe bias type lat lon
140 319147 2011 8 28 0.00 0.000000000 0.00000000 COOP 35.48667 -82.96833
141 319354 2011 8 28 0.02 0.001305799 -0.01869420 COOP 34.25722 -78.68722
142 319357 2011 8 28 0.00 0.045194085 0.04519409 COOP...
2011 Jun 27
0
OpenBUGS and glmmBUGS package
Dear All,
I am trying to use the glmmBUGS package with OpenBUGS under Linux.
Here is a short list of code I modified from the diseasemapping
package document:
----------------------------------------
library(diseasemapping)
data(popdata)
data(casedata)
model = getRates(casedata, popdata, ~age*sex, breaks=seq(0, 90, by=10) )
ontario = getSMR(popdata,model, casedata)
spplot(ontario, 'SMR')
library(spdep)
popDataAdjMat = poly2nb(ontario, ontario[["CSDUID"]])
library(glmmBUGS)
forBugs = glmmBUGS(formula=observed + logExpected...