max/ourMELONS
1
0
Fork 0
Code Issues 1 Pull requests Projects Releases Packages Wiki Activity Actions

Translated function into R

Browse source
This commit is contained in:
Waldir Leoncio 2020-01-14 13:06:36 +01:00
parent f8ca66f827
commit c72d2b6896
2 changed files with 15 additions and 19 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

20
R/admix1.R
View file

@ -536,22 +536,4 @@ admix1 <- function(tietue) {
# i2 = ceil(indeksi / npops);
# osuusTaulu(i1) = osuusTaulu(i1)-osuus;
# osuusTaulu(i2) = osuusTaulu(i2)+osuus;
# %-------------------------------------------------------------------------
# function [osuusTaulu, logml] = etsiParas(osuus, osuusTaulu, omaFreqs, logml)
# ready = 0;
# while ready ~= 1
# muutokset = laskeMuutokset4(osuus, osuusTaulu, omaFreqs, logml);
# [maxMuutos, indeksi] = max(muutokset(1:end));
# if maxMuutos>0
# osuusTaulu = suoritaMuutos(osuusTaulu, osuus, indeksi);
# logml = logml + maxMuutos;
# else
# ready = 1;
# end
# end
# osuusTaulu(i2) = osuusTaulu(i2)+osuus;

14
R/etsiParas.R Normal file
View file

@ -0,0 +1,14 @@
etsiParas <- function = (osuus, osuusTaulu, omaFreqs, logml) {
ready <- 0;
while (ready != 1) {
muutokset <- laskeMuutokset4(osuus, osuusTaulu, omaFreqs, logml)
[maxMuutos, indeksi] = max(muutokset[1:end]) # TODO: how does this work on Matlab?
if (maxMuutos > 0) {
osuusTaulu <- suoritaMuutos(osuusTaulu, osuus, indeksi)
logml <- logml + maxMuutos
} else {
ready <- 1
}
}
return (c(osuusTaulu, logml))
}