Displaying 1 result from an estimated 1 matches for "tot_mll".
Did you mean:
tot_mles
2008 May 23
2
About Passing Arguments to Function
Hi,
Below I have a function mlogl_k,
later it's called with "nlm" .
__BEGIN__
vsamples<- c(14.7, 18.8, 14, 15.9, 9.7, 12.8)
mlogl_k <- function( k_func, x_func, theta_func, samp) {
tot_mll <- 0
for (comp in 1:k_func) {
curr_mll <- (- sum(dgamma(samp, shape = x_func,
scale=theta_func, log = TRUE)))
tot_mll <- tot_mll + curr_mll
}
tot_mll
}
# Calling the function above
mlogl_out <- nlm(mlogl_k, mean(vsamples), k_func =2, x_func = 1,
theta_func...