function 1 : log-likelihood to estimate logmu
f(logmu)
= function 2 : log-likelihood to estimate beta
f(beta)
= ![$M_1$](/wiki/pub/Main/WfccmAlgorithmHu/latex43a6778619ed77dfa7a456ba3bfc9384.png)
= data in group1 ![$M_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex203e6a91a1a89a7e76822f442244cc4b.png)
= data in group2
![$rchisq_1$](/wiki/pub/Main/WfccmAlgorithmHu/latex77b3ee29c1737eca25fa76e1c89fc3bf.png)
for ![$M_1$](/wiki/pub/Main/WfccmAlgorithmHu/latex43a6778619ed77dfa7a456ba3bfc9384.png)
, ![$rchisq_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex2da376b0a92f22942d40d5d6cbaf62d0.png)
for ![$M_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex203e6a91a1a89a7e76822f442244cc4b.png)
based on 10^6 time runs
![$Var_1$](/wiki/pub/Main/WfccmAlgorithmHu/latex2b7f0df083c7f4146465710245928456.png)
of ![$rchisq_1$](/wiki/pub/Main/WfccmAlgorithmHu/latex77b3ee29c1737eca25fa76e1c89fc3bf.png)
, ![$Var_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex451fe15d8b1dcca6dccfe1b50b51030f.png)
of ![$rchisq_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex2da376b0a92f22942d40d5d6cbaf62d0.png)
![$Quan_1$](/wiki/pub/Main/WfccmAlgorithmHu/latex56ae9462263e860b17f2364811c1423b.png)
, mean ![$Mean_1$](/wiki/pub/Main/WfccmAlgorithmHu/latex154908f221a0c8242070d5dd1b072315.png)
, variance ![$Variance_1$](/wiki/pub/Main/WfccmAlgorithmHu/latexe0339cd3273675e4dd35885cc4b1d66f.png)
of ![$Quan_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex8b222458cc296e17a09b38a48b15bb13.png)
, mean ![$Mean_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex650a673c98c3835840f1b7d435ba8261.png)
, variance ![$Variance_2$](/wiki/pub/Main/WfccmAlgorithmHu/latex14f0b4b81f723ba8d74945bf9d6a4c61.png)
of mean
= variance
= 1% quantile
= value at (pos-1) + diff_v*diff_pcoefficients
: residuals
= y - (intercept + slope * x) xd2hat
= variance of residuals - Var xinitial
= estimate value of nlm
lkdinitial
min
and estimate
of xinitiallkdinf
xbeta.mle
= xinitial[(m+1):(m+2)]xbar.mle
= exponent of xinitialxs2.mle
= exponent of xbeta.mle * xbar.mle ^ xbeta.mle