Displaying 1 result from an estimated 1 matches for "nofc".
Did you mean:
noc
2008 May 23
0
Est. Component Size with AIC/BIC under Gamma Distribution
...w_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 = ", bic_k, " - AIC = ", aic_k, "- MLL= ", maxlogl...