The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.

geoSAE

This function is an extension of the Small Area Estimation (SAE) model. Geoadditive Small Area Model is a combination of the geoadditive model with the Small Area Estimation (SAE) model, by adding geospatial information to the SAE model.

Authors

Ketut Karang Pradnyadika, Ika Yuni Wulansari

Maintainer

Ketut Karang Pradnyadika 221709776@stis.ac.id

Installation

You can install the released version of geoSAE from CRAN or find my github repository Github

Example

#Load the dataset for unit level
library(geoSAE)
data(dataUnit)
#Load the dataset for spline-2
data(zspline)
#Load the dataset for area level
data(dataArea)
#Construct the data frame
y       <- dataUnit$y
x1      <- dataUnit$x1
x2      <- dataUnit$x2
x3      <- dataUnit$x3
formula <- y~x1+x2+x3
zspline <- as.matrix(zspline[,1:6])
dom     <- dataUnit$area
xmean   <- cbind(1,dataArea[,3:5])
zmean   <- dataArea[,7:12]
number  <- dataUnit$number
area    <- dataUnit$area
data    <- data.frame(number, area, y, x1, x2, x3)
#Estimate EBLUP
eblup_geosae <- eblupgeo(formula, zspline, dom, xmean, zmean, data)
eblup_geosae$eblup
#>           [,1]
#>  [1,] 29.04625
#>  [2,] 33.43651
#>  [3,] 34.66706
#>  [4,] 33.81857
#>  [5,] 23.52744
#>  [6,] 22.89752
#>  [7,] 21.86852
#>  [8,] 21.26004
#>  [9,] 33.73404
#> [10,] 38.43505
#> [11,] 33.77393
#> [12,] 28.98660
#> [13,] 32.29918
#> [14,] 24.31817
#> [15,] 31.23797
 
#Estimate MSE
mse_geosae <- pbmsegeo(formula,zspline,dom,xmean,zmean,data,B=100)
#> 
#> Bootstrap procedure with B = 100 iterations starts.
mse_geosae$mse
#>  [1]  2.052566  1.978709  2.231913 10.926587  1.480916  4.157471  3.172412
#>  [8]  1.839984  2.466619  1.563998  3.051762 16.937056 16.238457  2.648152
#> [15]  5.538344
 
## eblup_geosae$eblup        #to see the result of EBLUPs with Geoadditive Small Area Model each area
## mse_geosae$mse            #to see the result of MSE with Geoadditive Small Area Model each area

References

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.