How to write ri, which includes 76 parameters? Because with this syntax r gives an error.
logLikFun <- function(param) { ri <- param[1:76] sigma1 <- param[77] sigma2 <- param[78] s1 <- sum(dnorm(ri-1/(1-z2)*(z5*(1+z3)+0.02-1)+ri-1/(1-z2)*(0.1*log(z4)+0.02-log(0.99)), mean = 0, sd = sigma1, log = TRUE)) s2 <- sum(dnorm(ri-1/(1-z2)*(z5*(1+z3)+0.02-1)+ri-1/(1-z2)*(0.1*log(z4)+0.02-log(0.99)), mean = 0, sd = sigma2, log = TRUE)) s <- s1 + s2 } mle <- maxLik(logLik = logLikFun, start = c(ri = 0.01, sigma1 = 0.1, sigma2 = 0.1)) summary(mle)
param
is a list , then$ r -e 'param<-c(1:78);ri<-param[1:76];print(ri)'
works fine. Please give the minimum reproducible example (by clicking edit ) - aleksandr barakin