2020-01-14 12:09:27 +01:00
|
|
|
#' @title computeIndLogml
|
|
|
|
|
#' @description Palauttaa yksilön logml:n, kun oletetaan yksilön alkuperät
|
|
|
|
|
#' määritellyiksi kuten osuusTaulu:ssa.
|
2020-01-14 17:36:14 +01:00
|
|
|
#' @param omaFreqs own Freqs?
|
|
|
|
|
#' @param osuusTaulu Percentage table?
|
2020-01-14 12:09:27 +01:00
|
|
|
#' @export
|
|
|
|
|
computeIndLogml <- function (omaFreqs, osuusTaulu) {
|
|
|
|
|
|
|
|
|
|
apu <- repmat(t(osuusTaulu), c(1, dim(omaFreqs)[2]))
|
|
|
|
|
apu <- c(apu) * omaFreqs # c() avoids deprecation error re. matrix ops
|
|
|
|
|
if (length(apu) > 1) {
|
|
|
|
|
apu <- colSums(as.matrix(apu))
|
|
|
|
|
} else {
|
|
|
|
|
apu <- sum(apu)
|
|
|
|
|
}
|
|
|
|
|
|
2020-01-14 17:36:14 +01:00
|
|
|
apu <- log(apu)
|
2020-01-14 12:09:27 +01:00
|
|
|
|
|
|
|
|
loggis <- sum(apu)
|
|
|
|
|
return (loggis)
|
|
|
|
|
}
|