Page 1 of 1

Anamorphose

PostPosted: Fri Feb 10, 2017 10:59 am
by bez
Pour la version dite « empirique » (linéaire par morceaux ?), la help mentionne des arguments (disc.val et min.ratio) qui n’existent plus. Seul l’argument sigma2e est disponible.

Pourquoi le résultat est-il affiché avec un nombre de discrétisation qui ne correspond pas à celui demandé ? E.g. 33 quand j’en demande 100 ; 17 quand j’en demande 50…soit, il semblerait, un tiers de ce qui est demandé.

Pourquoi l’intervalle de fluctuation de Z n’est-il pas respecté et décalé vers la gauche ?

Pourquoi le premier intervalle des valeurs de Y est-il atomique (calé sur -10) ? C’est le côté atomique qui me pose question.

L’appel à anam.plot() en fin de fonction anam.fit() n’est pas alimenté avec les arguments d’entrée comme ndisc, azmin/max, aymin/max, sigma2e. Or anam.plot() appelle anam.eval() qui utilise donc les valeurs par défaut de ces arguments. Il me semble que c’est une erreur et qu’il faudrait passer les arguments de anam.fit() dans anam.plot().

Re: Anamorphose

PostPosted: Mon Mar 06, 2017 10:35 am
by bez
The corrections are not fully satisfactory.

arguments aymin-max and azmin.max are now compulsory for the graph to be represented.
Running anam.fit() without specifying them is OK for the computation but not for the representation.
In addition, the x and y limits of the graph are not consistent by default.
a <- db.create(x1=runif(1000),z1=rlnorm(1000))
aa <- anam.fit(a) # ==> pas de représentation graphique
aa <- anam.fit(a,aymin=-5,ayma=5,azmin=0,azmax=30) # graphe mais non calé sur les valeurs fournies.
aa <- anam.fit(a,aymin=-5,ayma=5,azmin=0,azmax=30,xlim=c(-7,7),ylim=c(0,35)) # graphe complet visible

In addition ndisc is used for type="emp" AND for anam.plot. The ndisc argument should thus also be passed to anam.plot.
So the graph is the same whatever the ndisc value.

Gaussian anamorphosis
the specification of aymin-max and azmin.max are not honoured strictly in the computations. The reasons are propbably the way intervals are handled. However when asking for azmin=0 and azmax=30, I think these these values should be honoured
aa <- anam.fit(a,aymin=-5,ayma=5,azmin=0,azmax=30,xlim=c(-7,7),ylim=c(0,35))
aaa <- anam.y2z(aa,db.create(x1=runif(100),z1=seq(-10,10,length=100)))
aaa[,4]
The values range from 0.0188 and 38.24 instead of ranging from 0 to 30.

Empirical anamorphosis
aa <- anam.fit(a,type="emp",ndisc=100,sigma2e=0,xlim=c(-7,7),ylim=c(0,35))
==> no graph is produced
aa[,4]
==> wrong max for z even if no additional variance is provided (sigma2e = 0)
the input variable ranges from 0.07 to 19.6
range(a[,3])
[1] 0.07119225 19.56747995
the anamorphose ranges from 0.07 to 18.6
range(aaa[,4])
[1] 0.07118513 18.57276428