Displaying 1 result from an estimated 1 matches for "thetastart".
Did you mean:
thetastar
2008 May 23
0
Est. Component Size with AIC/BIC under Gamma Distribution
...smpl), theta_func=theta, samp=smpl)
# M-STEP
if (new_maxlogl$minimum > best) {
best <- new_maxlogl$minimum
}
}
best
}
# main program
# my samples
vsamples<- c(14.7, 18.8, 14, 15.9, 9.7, 12.8)
# initialize
start_ <- 300
error_ <- 0.001
thetastart <- 1
k <- 5
# compute AIC/BIC for k component
for (nofc in 1:k {
cat("k = ", nofc, "\n")
maxlogl <- -(mlogl_process(vsamples,error_,thetastart))
bic_k <- find_bic(maxlogl,vsamples,nofc)
aic_k <- find_aic(maxlogl,vsamples,nofc)
cat("BIC = ", b...