30 lines
1 KiB
R
30 lines
1 KiB
R
#' @title Etsi Paras
|
|
#' @description Search for the best?
|
|
#' @param osuus Percentages?
|
|
#' @param omaFreqs own Freqs?
|
|
#' @param osuusTaulu Percentage table?
|
|
#' @param logml log maximum likelihood
|
|
etsiParas <- function(osuus, osuusTaulu, omaFreqs, logml) {
|
|
ready <- 0
|
|
while (ready != 1) {
|
|
muutokset <- admix1_muutokset$new()
|
|
muutokset <- muutokset$laskeMuutokset4(osuus, osuusTaulu, omaFreqs, logml)
|
|
|
|
# Work around R's base::max() limitation on complex numbers
|
|
if (any(vapply(muutokset, class, vector("character", 1)) == "complex")) {
|
|
maxRe <- base::max(Re(as.vector(muutokset)))
|
|
maxIm <- base::max(Im(as.vector(muutokset)))
|
|
maxMuutos <- complex(real = maxRe, imaginary = maxIm)
|
|
} else {
|
|
maxMuutos <- base::max(as.vector(muutokset))
|
|
}
|
|
indeksi <- which(muutokset == maxMuutos)
|
|
if (Re(maxMuutos) > 0) {
|
|
osuusTaulu <- suoritaMuutos(osuusTaulu, osuus, indeksi)
|
|
logml <- logml + maxMuutos
|
|
} else {
|
|
ready <- 1
|
|
}
|
|
}
|
|
return(c(osuusTaulu, logml))
|
|
}
|