Re: [R] What to do with positive likelihoods

Michael, It's not entirely clear what you are doing from the text (and the code is incompplete), so just a few remarks. I assume you are after the log likelihood rather than the likelihood (otherwise the issue of positive/negative would not show up, dnorm always returns positive numbers). Taking

[R] What to do with positive likelihoods

Hi all, I'll preface this with saying I've gone through the archives, and am still in need of some help. I've been using this likelihood model with mean = 0 and s.d. = sqrt( (c + ( 1 / N1 ) + ( 1 / N2 ) ) * x * ( 1 - x )), where c is a genetic drift parameter (usually very small, like between