Dear All, I am evaluating the value of loglikelihood and it ends up with the sum of tiny numbers. Below is an example: suppose I would like to calculate sum_i (log (sum_j x [i, j] )), the index of log (x) is in the range, say (-2000, 0). I am aware that exp(-744.5) will be expressed as 0 in 32 bit R and exp Is there a way to improve the result?
R example: powd <- sample(-2000:0, 100, replace=T) # the power of x [i, j] x <- matrix(exp(powd),10) # the value of x sum(log(rowSums(x))) # sum Thank you for your help. Best wishes, Jie [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.