Method to create a.wght component from the LKinfo object.

Description

This method takes a vector or more complicated object and based on the geometry creates a list with the a.wght information.

Usage

LKrigSetupAwght(object, ...)
## Default S3 method:
LKrigSetupAwght(object, ...)
## S3 method for class 'LKRectangle'
LKrigSetupAwght(object, ...)

Arguments

Details

The simplest function of this method is to convert the sequence of a.wght values into a list. If only a scalar a.wght value is supplied then the default method just repeats this for each level. The attribute fastNormalize (either TRUE or FALSE) is attached to this list to indicate how the marginal variance of the process should be found.

 	 	
  LKinfo<- LKrigSetup( x,LKGeometry="LKInterval", alpha=c( 1,.2,.01),
                   nlevel=3, a.wght=4.5, NC= 3)
  LKrigSetupAwght( LKinfo)
  
[[1]]
[1] 4.5

[[2]]
[1] 4.5

[[3]]
[1] 4.5

attr(,"fastNormalize")
[1] FALSE

Currently the only geometry with fastNormalization being TRUE is for a rectangular domain.

For the LKRectangle geometry, however, more complicated anisotropic and nonstationary a.wght specifications are possible. See LKrig for details. Also in the case that the fastNormalization is TRUE for rectangles several more attributes are added to the a.wght list that precompute some matrices of the SAR.

Value

A list with nlevel components. The attribute fastNormalize is added to this list. In the case that the geometry is LKRectangle several more attributes are added indicating the type of covariance model and possibly an eigen decomposition of the SAR matrix exploiting Kronecker products.

Author(s)

Doug Nychka

See Also

LKrigSetup, LKrigSetupAlpha, LKrigSAR, LKrig

Examples

  x<- cbind( c(0,1))
  LKinfo<- LKrigSetup( x,LKGeometry="LKInterval", alpha=c( 1,.2,.01),
                   nlevel=3, a.wght=4.5, NC= 3)
  a.wghtList<- LKrigSetupAwght( LKinfo)
  
  x<- cbind( c(0,1), c(0,1))
  LKinfo<- LKrigSetup( x, alpha=c( 1,.2,.01),
                   nlevel=3, a.wght=4.5, NC= 3)
  a.wghtList<- LKrigSetupAwght( LKinfo)
# see   
  names(attributes( a.wghtList))
 
  

Results

R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(LatticeKrig)
Loading required package: spam
Loading required package: grid
Spam version 1.3-0 (2015-10-24) is loaded.
Type 'help( Spam)' or 'demo( spam)' for a short introduction 
and overview of this package.
Help for individual functions is also obtained by adding the
suffix '.spam' to the function name, e.g. 'help( chol.spam)'.

Attaching package: 'spam'

The following objects are masked from 'package:base':

    backsolve, forwardsolve

Loading required package: fields
Loading required package: maps

 # maps v3.1: updated 'world': all lakes moved to separate new #
 # 'lakes' database. Type '?world' or 'news(package="maps")'.  #


> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/LatticeKrig/LKrigSetupAwght.Rd_%03d_medium.png", width=480, height=480)
> ### Name: LKrigSetupAwght
> ### Title: Method to create a.wght component from the 'LKinfo' object.
> ### Aliases: LKrigSetupAwght LKrigSetupAwght.default
> ###   LKrigSetupAwght.LKRectangle
> ### Keywords: spatial
> 
> ### ** Examples
> 
>   x<- cbind( c(0,1))
>   LKinfo<- LKrigSetup( x,LKGeometry="LKInterval", alpha=c( 1,.2,.01),
+                    nlevel=3, a.wght=4.5, NC= 3)
>   a.wghtList<- LKrigSetupAwght( LKinfo)
>   
>   x<- cbind( c(0,1), c(0,1))
>   LKinfo<- LKrigSetup( x, alpha=c( 1,.2,.01),
+                    nlevel=3, a.wght=4.5, NC= 3)
>   a.wghtList<- LKrigSetupAwght( LKinfo)
> # see   
>   names(attributes( a.wghtList))
[1] "fastNormDecomp" "fastNormalize"  "first.order"    "stationary"    
>  
>   
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>