funGP {DSWE}R Documentation

Function comparison using Gaussian Process and Hypothesis testing

Description

Function comparison using Gaussian Process and Hypothesis testing

Usage

funGP(
  datalist,
  xCol,
  yCol,
  confLevel = 0.95,
  testset,
  limitMemory = TRUE,
  opt_method = "nlminb",
  sampleSize = list(optimSize = 500, bandSize = 5000),
  rngSeed = 1
)

Arguments

datalist

A list of data sets to compute a function for each of them

xCol

A numeric or vector stating the column number of covariates

yCol

A numeric value stating the column number of target

confLevel

A single value representing the statistical significance level for constructing the band

testset

Test points at which the functions will be compared

limitMemory

A boolean (True/False) indicating whether to limit the memory use or not. Default is true. If set to true, 5000 datapoints are randomly sampled from each dataset under comparison for inference.

opt_method

A string specifying the optimization method to be used for hyperparameter estimation. Current options are: 'L-BFGS-B', 'BFGS', and 'nlminb'. Default is set to 'nlminb'.

sampleSize

A named list of two integer items: optimSize and bandSize, denoting the sample size for each dataset for hyperparameter optimization and confidence band computation, respectively, when limitMemory = TRUE. Default value is list(optimSize = 500, bandSize = 5000).

rngSeed

Random seed for sampling data when limitMemory = TRUE. Default is 1.

Value

a list containing :

References

Prakash, A., Tuo, R., & Ding, Y. (2020). "Gaussian process aided function comparison using noisy scattered data." arXiv preprint arXiv:2003.07899. <https://arxiv.org/abs/2003.07899>.

Examples


datalist = list(data1[1:100,], data2[1:100, ])
xCol = 2
yCol = 7
confLevel = 0.95
testset = seq(4,10,length.out = 20)
function_diff = funGP(datalist, xCol, yCol, confLevel, testset)


[Package DSWE version 1.5.1 Index]