Anamorphose

If you have problems using RGeostats package or if you simply need some help on an unlisted subject.

Anamorphose

Postby bez » Fri Feb 10, 2017 10:59 am

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().
bez
 
Posts: 13
Joined: Tue Oct 23, 2012 10:31 am

Re: Anamorphose

Postby bez » Mon Mar 06, 2017 10:35 am

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
bez
 
Posts: 13
Joined: Tue Oct 23, 2012 10:31 am


Return to Other Troubleshooting

Who is online

Users browsing this forum: No registered users and 1 guest

cron