plotRetentionArea: Plots the convex increasing retention area

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

View source: R/plotRetentionArea.R

Description

For two evidence factors this function produces a grey-scale plot of the retention set, i.e. the set of sensitive bias levels.

Usage

1

Arguments

retentionBrd

Output of retentionBrd function for two evidence factors.

Gamlist

List of two vectors of gamma values for two evidence factors.

...

Further arguments to be passed to method image and legend used for plotting.

Details

For 2 evidence factors plotRetentionArea produces a grey-scale plot of the retention set, i.e. the set of sensitive bias levels.

Value

The output of the function is plot.

Author(s)

Bikram Karmakar

References

Karmakar, B., French, B., Sadakane, A. and Small, D. S. (2016) Integrating the evidence from evidence factors in observational studies.

See Also

plotRejDecbyAssm.

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
28
data(Plist)
data(Gamlist)
rbrd = retentionBrd(Plist, Gamlist)
plotRetentionArea(rbrd, Gamlist)

rbrd = retentionBrd(Plist, Gamlist, method = "Truncated")
plotRetentionArea(rbrd, Gamlist)

## mean tail moment data example
library(sensitivitymv)

data(mtm)
Gamseq <- seq(1, 15, by = 0.2)
Gamlist <- list(Gamseq, Gamseq)
Plist <- list(c(), c())

for(gam in Gamseq){
	Plist[[1]] = c(Plist[[1]], senmv(-mtm,gamma=gam,trim=1)$pval)
	Plist[[2]] = c(Plist[[2]], senmv(-mtm[,2:3],gamma=gam,trim=1)$pval)
}	

# Fisher's combination method
rbrd <- retentionBrd(Plist, Gamlist)
plotRetentionArea(rbrd, Gamlist)

# truncated product combination
rbrd <- retentionBrd(Plist, Gamlist, method = "Truncated", talpha = .5)
plotRetentionArea(rbrd, Gamlist)

evidenceFactors documentation built on Feb. 20, 2020, 5:07 p.m.