Skip to content

CensSpatial

Censored Spatial Models

v3.6 · Jan 24, 2023 · GPL (>= 2)
CRAN GitHub Mirror

Description

It fits linear regression models for censored spatial data. It provides different estimation methods as the SAEM (Stochastic Approximation of Expectation Maximization) algorithm and seminaive that uses Kriging prediction to estimate the response at censored locations and predict new values at unknown locations. It also offers graphical tools for assessing the fitted model. More details can be found in Ordonez et al. (2018) <doi:10.1016/j.spasta.2017.12.001>.

Downloads

308

Last 30 days

12962nd

814

Last 90 days

3.3K

Last year

Trend: +3.4% (30d vs prior 30d)

CRAN Check Status

2 ERROR
2 NOTE
10 OK
Show all 14 flavors
Flavor Status
r-devel-linux-x86_64-debian-clang NOTE
r-devel-linux-x86_64-debian-gcc NOTE
r-devel-linux-x86_64-fedora-clang ERROR
r-devel-linux-x86_64-fedora-gcc ERROR
r-devel-macos-arm64 OK
r-devel-windows-x86_64 OK
r-oldrel-macos-arm64 OK
r-oldrel-macos-x86_64 OK
r-oldrel-windows-x86_64 OK
r-patched-linux-x86_64 OK
r-release-linux-x86_64 OK
r-release-macos-arm64 OK
r-release-macos-x86_64 OK
r-release-windows-x86_64 OK
Check details (6 non-OK)
NOTE r-devel-linux-x86_64-debian-clang

CRAN incoming feasibility

Maintainer: ‘Alejandro Ordonez <ordonezjosealejandro@gmail.com>’

No Authors@R field in DESCRIPTION.
Please add one, modifying
  Authors@R: c(person(given = "Alejandro",
                      family = "Ordonez",
                      role = c("aut", "cre"),
                      email = "ordonezjosealejandro@gmail.com"),
               person(given = c("Christian", "E."),
                      family = "Galarza",
                      role = "aut"),
               person(given = c("Victor", "H."),
                      family = "Lachos",
                      role = "aut"))
as necessary.
NOTE r-devel-linux-x86_64-debian-gcc

CRAN incoming feasibility

Maintainer: ‘Alejandro Ordonez <ordonezjosealejandro@gmail.com>’

No Authors@R field in DESCRIPTION.
Please add one, modifying
  Authors@R: c(person(given = "Alejandro",
                      family = "Ordonez",
                      role = c("aut", "cre"),
                      email = "ordonezjosealejandro@gmail.com"),
               person(given = c("Christian", "E."),
                      family = "Galarza",
                      role = "aut"),
               person(given = c("Victor", "H."),
                      family = "Lachos",
                      role = "aut"))
as necessary.
NOTE r-devel-linux-x86_64-fedora-clang

dependencies in R code

Namespace in Imports field not imported from: ‘Rcpp’
  All declared Imports should be used.
ERROR r-devel-linux-x86_64-fedora-clang

examples

Running examples in ‘CensSpatial-Ex.R’ failed
The error most likely occurred in:

> ### Name: algnaive12
> ### Title: Naive 1 and Naive 2 method for spatial prediction.
> ### Aliases: algnaive12
> ### Keywords: Spatial Censored SAEM
> 
> ### ** Examples
> 
> ###simulated coordinates
> n<-200 ### sample size for estimation.
> n1=100 ### number of observation used in the prediction.
> 
> ###simulated coordinates
> r1=sample(seq(1,30,length=400),n+n1)
> r2=sample(seq(1,30,length=400),n+n1)
> coords=cbind(r1,r2)### total coordinates (used in estimation and prediction).
> 
> coords1=coords[1:n,]####coordinates used for estimation.
> 
> type="matern"### covariance structure.
> 
> xtot<-cbind(1,runif((n+n1)),runif((n+n1),2,3))## X matrix for estimation and prediction.
> xobs=xtot[1:n,]## X matrix for estimation.
> 
> ###simulated data
> obj=rspacens(cov.pars=c(3,.3,0),beta=c(5,3,1),x=xtot,coords=coords,kappa=1.2,
+ cens=0.25,n=(n+n1),n1=n1,cov.model=type,cens.type="left")
> 
> data2=obj$datare
> data2[,4:5]=xobs[,-1]
> 
> cc=obj$cc
> y=obj$datare[,3]
> cutoff=rep(obj$cutoff,length(y[cc==1]))
> 
> 
> aux2=algnaive12(data=data2,cc=obj$cc,covar=TRUE,covar.col=4:5,
+ copred=obj$coords1,thetaini=c(.1,.2),y.col=3,coords.col=1:2,
+ fix.nugget=TRUE,nugget=0,kappa=1.2,cutoff=cutoff,trend=~V4+V5,
+ cov.model=type)
OMP: Warning #96: Cannot form a team with 24 threads, using 2 instead.
OMP: Hint Consider unsetting KMP_DEVICE_THREAD_LIMIT (KMP_ALL_THREADS), KMP_TEAMS_THREAD_LIMIT, and OMP_THREAD_LIMIT (if any are set).
NOTE r-devel-linux-x86_64-fedora-gcc

dependencies in R code

Namespace in Imports field not imported from: ‘Rcpp’
  All declared Imports should be used.
ERROR r-devel-linux-x86_64-fedora-gcc

examples

Running examples in ‘CensSpatial-Ex.R’ failed
The error most likely occurred in:

> ### Name: algnaive12
> ### Title: Naive 1 and Naive 2 method for spatial prediction.
> ### Aliases: algnaive12
> ### Keywords: Spatial Censored SAEM
> 
> ### ** Examples
> 
> ###simulated coordinates
> n<-200 ### sample size for estimation.
> n1=100 ### number of observation used in the prediction.
> 
> ###simulated coordinates
> r1=sample(seq(1,30,length=400),n+n1)
> r2=sample(seq(1,30,length=400),n+n1)
> coords=cbind(r1,r2)### total coordinates (used in estimation and prediction).
> 
> coords1=coords[1:n,]####coordinates used for estimation.
> 
> type="matern"### covariance structure.
> 
> xtot<-cbind(1,runif((n+n1)),runif((n+n1),2,3))## X matrix for estimation and prediction.
> xobs=xtot[1:n,]## X matrix for estimation.
> 
> ###simulated data
> obj=rspacens(cov.pars=c(3,.3,0),beta=c(5,3,1),x=xtot,coords=coords,kappa=1.2,
+ cens=0.25,n=(n+n1),n1=n1,cov.model=type,cens.type="left")
> 
> data2=obj$datare
> data2[,4:5]=xobs[,-1]
> 
> cc=obj$cc
> y=obj$datare[,3]
> cutoff=rep(obj$cutoff,length(y[cc==1]))
> 
> 
> aux2=algnaive12(data=data2,cc=obj$cc,covar=TRUE,covar.col=4:5,
+ copred=obj$coords1,thetaini=c(.1,.2),y.col=3,coords.col=1:2,
+ fix.nugget=TRUE,nugget=0,kappa=1.2,cutoff=cutoff,trend=~V4+V5,
+ cov.model=type)

Check History

ERROR 10 OK · 2 NOTE · 0 WARNING · 2 ERROR · 0 FAILURE Mar 10, 2026
NOTE r-devel-linux-x86_64-debian-clang

CRAN incoming feasibility

Maintainer: ‘Alejandro Ordonez <ordonezjosealejandro@gmail.com>’

No Authors@R field in DESCRIPTION.
Please add one, modifying
  Authors@R: c(person(given = "Alejandro",
                      family = "Ordonez",
                      role = c("aut", 
...[truncated]...
(given = c("Christian", "E."),
                      family = "Galarza",
                      role = "aut"),
               person(given = c("Victor", "H."),
                      family = "Lachos",
                      role = "aut"))
as necessary.
NOTE r-devel-linux-x86_64-debian-gcc

CRAN incoming feasibility

Maintainer: ‘Alejandro Ordonez <ordonezjosealejandro@gmail.com>’

No Authors@R field in DESCRIPTION.
Please add one, modifying
  Authors@R: c(person(given = "Alejandro",
                      family = "Ordonez",
                      role = c("aut", 
...[truncated]...
(given = c("Christian", "E."),
                      family = "Galarza",
                      role = "aut"),
               person(given = c("Victor", "H."),
                      family = "Lachos",
                      role = "aut"))
as necessary.
ERROR r-devel-linux-x86_64-fedora-clang

dependencies in R code

Namespace in Imports field not imported from: ‘Rcpp’
  All declared Imports should be used.
ERROR r-devel-linux-x86_64-fedora-gcc

dependencies in R code

Namespace in Imports field not imported from: ‘Rcpp’
  All declared Imports should be used.

Dependency Network

Dependencies Reverse dependencies geoR Rcpp mvtnorm optimx (>= 2021.10-12) tmvtnorm msm psych numDeriv (>= 2.11.1) raster moments lattice tlrmvnmvt CensSpatial

Version History

new 3.6 Mar 10, 2026
updated 3.6 ← 3.5 diff Jan 23, 2023
updated 3.5 ← 3.0 diff Jan 16, 2023
updated 3.0 ← 2.59 diff Aug 14, 2022
updated 2.59 ← 2.58 diff Oct 15, 2021
updated 2.58 ← 2.1 diff May 3, 2020
updated 2.1 ← 1.3 diff Jan 8, 2020
updated 1.3 ← 1.2 diff Mar 2, 2017
updated 1.2 ← 1.1 diff Dec 3, 2016
updated 1.1 ← 1.0 diff Oct 16, 2016
new 1.0 Oct 7, 2016