mathematical optimization - Optimize a function several times R code -
i have function. can estimate parameters easily.
sex <- c("f","m","f","m","f") age <- rnorm(5,28,1.2) dat <- data.frame(sex,age) myfun <- function(par, x1,x2){ -sum(log(exp(par[1]*x1+par[2]*x2))) } myf <- optim(myfun, par=c(0.1,0.4), x1=dat$age,x2=as.numeric(dat$sex))$par
i want optimize function 10 times. if use replicate(10,myf)
gives same values 10 times. guess give different parameter estimates every time due age <- rnorm(5,28,1.2)
. want using loop, how can proceed?
i'm not sure need, work you?
myfun <- function(par, x1,x2){ -sum(log(exp(par[1]*x1+par[2]*x2))) } sex <- c("f","m","f","m","f") (i in 1:10){ age <- rnorm(5,28,1.2) dat <- data.frame(sex,age) print(optim(myfun, par=c(0.1,0.4), x1=dat$age,x2=as.numeric(dat$sex))$par) }
Comments
Post a Comment