[R] EM algorithm

[Elena 5/12]

l=∑_i∑_k log[ (Se(i))^(t_ik ) [1-Se(i) ]^(1-t_ik ) p+(Sp(i))^(1-t_ik)
[1-Sp(i) ]^(t_ik ) (1-p)]

This is my observed loglikelihood. It is a sum over i (there are I tests)
and over k (K patients). t_ik is the result for the kth patient on the ith
test.

Se(i)= sensitivity, Sp (i)= specificity and p=prevalence are my unknown
parameters.

Do you know how to code this loglikelihood and EM algorithm in R? I can't
get it to work.

Thanks!


-- 
View this message in context: http://www.nabble.com/EM-algorithm-tp24595358p24596434.html
Sent from the R help mailing list archive at Nabble.com.