calcLikelihood: Calculate the likelihood for the samples in a MADproject...

Description Usage Arguments Details Value Methods (by class) Examples

Description

calcLikelihood returns an updated MADproject with the likelihood values based on the observation and realization data in the MADproject or, optionally, a subset thereof.

Usage

1
2
3
4
5
6
7
8
9
calcLikelihood(proj, dsubset, num_realz = max(proj@realizations$rid),
  samples = 1:proj@numSamples)

## S4 method for signature 'MADproject,numeric'
calcLikelihood(proj, dsubset,
  num_realz = max(proj@realizations$rid), samples = 1:proj@numSamples)

## S4 method for signature 'MADproject,ANY'
calcLikelihood(proj, num_realz, samples)

Arguments

proj

The MADproject object with data read from the MAD# databases.

dsubset

The subset of inversion data to use for the likelihood calculations.

num_realz

The number of realizations to use in the likelihood calculation (defaults to all in the realizations slot)

samples

A vector of sample IDs for which to calculate likelihood values (defaults to all available in the realizations slot)

Details

The likelihood calculation utilizes the np package for non- parametric density estimation with all inversion data as dependent (i.e. multivariate likelihood distributions are estimated).

Value

proj The updated MADproject object with a filled likelihood slot.

Methods (by class)

Examples

1
2
data(pumping)
pumping <- calcLikelihood(pumping, 100) #Inversion data as time step 100

anchoredDistr documentation built on June 20, 2017, 9:14 a.m.