Surrogate (version 3.3.1)

plot MaxEnt ContCont: Plots the sensitivity-based and maximum entropy based Individual Causal Association when S and T are continuous outcomes in the single-trial setting

Description

This function provides a plot that displays the frequencies or densities of the individual causal association (ICA; \(rho[Delta]\)) as identified based on the sensitivity- (using the functions ICA.ContCont) and maximum entropy-based (using the function MaxEntContCont) approaches.

Usage

# S3 method for MaxEntContCont
plot(x, Type="Freq", Xlab, col, 
Main, Entropy.By.ICA=FALSE, ...)

Arguments

x

An object of class MaxEntContCont. See MaxEntContCont.

Type

The type of plot that is produced. When Type="Freq", the Y-axis shows frequencies of ICA. When Type="Density", the density is shown. Default Type="Freq".

Xlab

The legend of the X-axis of the plot.

col

The color of the bins (frequeny plot) or line (density plot). Default col <- c(8).

Main

The title of the plot.

Entropy.By.ICA

Plot with ICA on Y-axis and entropy on X-axis.

...

Other arguments to be passed to plot()

Author

Wim Van der Elst, Ariel Alonso, Paul Meyvisch, & Geert Molenberghs

References

Add

See Also

ICA.ContCont, MaxEntContCont

Examples

if (FALSE)  #time-consuming code parts
# Compute ICA for ARMD dataset, using the grid  
# G={-1, -.80, ..., 1} for the undidentifiable correlations

ICA <- ICA.ContCont(T0S0 = 0.769, T1S1 = 0.712, S0S0 = 188.926, 
S1S1 = 132.638, T0T0 = 264.797, T1T1 = 231.771, 
T0T1 = seq(-1, 1, by = 0.2), T0S1 = seq(-1, 1, by = 0.2), 
T1S0 = seq(-1, 1, by = 0.2), S0S1 = seq(-1, 1, by = 0.2))

# Identify the maximum entropy ICA
MaxEnt_ARMD <- MaxEntContCont(x = ICA, S0S0 = 188.926, 
S1S1 = 132.638, T0T0 = 264.797, T1T1 = 231.771)

  # Explore results using summary() and plot() functions
summary(MaxEnt_ARMD)
plot(MaxEnt_ARMD)
plot(MaxEnt_ARMD, Entropy.By.ICA = TRUE)