robCompositions (version 2.4.1)

fcenLRinv: Inverse of fcenLR transformations (functional)

Description

Inverse of fcenLR transformations

Usage

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

Value

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

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

Author

R. Talska, A. Menafoglio, K. Hron, J. J. Egozcue, J. Palarea-Albaladejo

Details

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

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")