■

3本サブミット．関係各位おつかれさまです．

LL<-colSums(Ctakashi*log(Ptakashi) + Cmituko*log(Pmituko) + Cf_matu*log(Pf_matu) + Cjyasuko*log(Pjyasuko) + Cjoepla*log(Pjoepla))
res<-optim(b0,fr, method = "BFGS", hessian = TRUE, control=list(fnscale=-1))

BFGSって，L-BFGS-Bじゃないよね．

optim(par, fn, gr = NULL,
method = c("Nelder-Mead", "BFGS", "CG", "L-BFGS-B", "SANN"),
lower = -Inf, upper = Inf,
control = list(), hessian = FALSE, ...)