fcenLRinv: Inverse of fcenLR transformations (functional)

In matthias-da/robCompositions: Compositional Data Analysis

Inverse of fcenLR transformations (functional)

Description

Inverse of fcenLR transformations

Usage

fcenLRinv(z, z_step, fcenLR, k = 1)

Arguments

z	grid of points defining the abscissa
z_step	step of the grid of the abscissa
fcenLR	grid evaluation of (i) fcenLR[lambda] transformed lambda-density, (ii) fcenLR[u] transformed P-density, (iii) fcenLR[P] transformed P-density
k	value of the integral of density; if k=1 it returns a unit-integral representation of density

Details

By default, it returns a unit-integral representation of density.

Value

out ... grid evaluation of (i) lambda-density in B2(lambda), (ii) P-density in unweighted B2(lambda), (iii) P-density in B2(P)

Author(s)

R. [e-mail], A. Menafoglio, K. [e-mail], J. J. Egozcue, J. Palarea-Albaladejo

Examples

# Example (normal density)
t = seq(-4.7,4.7, length = 1000)
t_step = diff(t[1:2])

mean = 0; sd = 1.5
f = dnorm(t, mean, sd)
f1 = f/trapzc(t_step,f)

f.fcenLR = fcenLR(t,t_step,f) 
f.fcenLRinv = fcenLRinv(t.fine,t_step,f.fcenLR)

plot(t,f.fcenLR, type="l",las=1, ylab="fcenLR(density)", 
  cex.lab=1.2,cex.axis=1.2, col="darkblue",lwd=2)
abline(h=0, col="red")

plot(t,f.fcenLRinv, type="l",las=1, 
  ylab="density",cex.lab=1.2,cex.axis=1.2, col="darkblue",lwd=2,lty=1)
lines(t,f1,lty=2,lwd=2,col="gold")  

matthias-da/robCompositions documentation built on Jan. 15, 2024, 11:24 p.m.