qadqd.QMC: ~~function to do ... ~~

Description Usage Arguments Details Value Note Author(s) References See Also Examples

View source: R/quaddisc.R

Description

~~ A concise (1-5 lines) description of what the function does. ~~

Usage

1
qadqd.QMC(dP, fnH, mZ, sForm = "v", lExtraArgsKernel, lArgsUniroot = list(interval = c(0, 100)), ...)

Arguments

dP

~~Describe dP here~~

fnH

~~Describe fnH here~~

mZ

~~Describe mZ here~~

sForm

~~Describe sForm here~~

lExtraArgsKernel

~~Describe lExtraArgsKernel here~~

lArgsUniroot

~~Describe lArgsUniroot here~~

...

~~Describe ... here~~

Details

~~ If necessary, more details than the description above ~~

Value

~Describe the value returned If it is a LIST, use

comp1

Description of 'comp1'

comp2

Description of 'comp2'

...

Note

~~further notes~~

~Make other sections like Warning with Warning .... ~

Author(s)

~~who you are~~

References

~put references to the literature/web site here ~

See Also

~~objects to See Also as help, ~~~

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function (dP, fnH, mZ, sForm = "v", lExtraArgsKernel, lArgsUniroot = list(interval = c(0, 
    100)), ...) 
{
    if ((dP <= 0) | (dP >= 1)) 
        return("wrong argument (should be in (0, 1))")
    if (!is.null(lExtraArgsKernel)) 
        lArgs <- c(list(mZ = mZ, fnH = fnH), lExtraArgsKernel)
    else lArgs <- list(mZ = mZ, fnH = fnH)
    if (sForm == "u") 
        mH <- do.call("matHu.QMC", lArgs)
    if (sForm == "v") 
        mH <- do.call("matH.QMC", lArgs)
    vLambda <- eigen(mH, symmetric = TRUE)$values
    fnObj <- function(dC) {
        lP <- pgausschaos(dC, vLambda, ...)
        if (lP$iFault != 0) 
            print("warning: error in computation of pgausschaos")
        return(lP$dResult - dP)
    }
    lArgsUniroot <- c(f = fnObj, lArgsUniroot)
    return(do.call("uniroot", lArgsUniroot))
  }

quaddisc documentation built on May 2, 2019, 4:50 p.m.