Displaying 1 result from an estimated 1 matches for "betatheldaplus".
2010 Sep 21
2
Need help for EM algorithm ASAP !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
...dastar[i]*gammastar))/(gammastar*klambda*theldastar[i]*gammastar))
ti[i]<-((gammastar*klambda2*(theldastar[i]*gammastar))/(theldastar[i]*klambda*theldastar[i]*gammastar))
zimean<-sum(zi)/n
timean<-sum(ti)/n
mutheldaplus<-(zimean*(1/n)* sum((ti[i]*y[i])-mean(y)))/((zimean*timean)-1)
betatheldaplus<- sum(y[i]- mutheldaplus)/(n*zimean)
sigmatheldaplus<-((1/n)*sum((ti[i]*((y[i]-mutheldaplus)^2))-(2*betatheldaplus*(y[i]-mutheldaplus))-((betatheldaplus^2)*zi[i])))
print(muthelda)
print(mutheldaplus)
print(betathelda)
print(betatheldaplus)
print(sigmathelda)
print(sigmatheldaplus)
return(...