R语言在maxLik函数中，如何给待估计的参数设置条件

1. data<-read.xlsx("E:\\1.my construction\\vasicek\\mle.xlsx",sheet=1)
   
2. rt<-dataBMJP02Y(RA)`;rt<-rt[-length(rt)]/100;rT<-dataBMJP02Y(RA)`[-1]/100;n<-length(rt)
   
3. pai<-3.14159;e<-exp(1)
   
4. LL<-function(par)
   
5. {
   
6.   beta<-par[1];alphh<-par[2];sigma<-par[3]
   
7.   f<-0.5*(1-n)*(log(2*pai)+log((1-e^(-2*alphh))*sigma^2/(2*beta)))-sum(rT-beta-(rt-beta)*e^(-alphh))^2/((1-e^(-2*alphh))*sigma^2/(alphh))
   
8. 
   
9.   return(f)
   
10. }
    
11. library(maxLik)
    
12. mle<-maxLik(LL, start=c(beta=0.0015,alphh=0.035,sigma=0.00500275))
    
13. summary(mle)

复制代码

我想把sigma的值设定为大于等于0.0001，其他两个参数也都要大于0，这个在maxLik()里面怎么写呢？

我知道在optim中可以这么写，optim(par,LL,lower=c(-Inf, -Inf, 0.0001),method="L-BFGS-B",hessian=TRUE,control=list(fnscale=-1))，当然要先定义par<-c(0.0015,0.035,0.00500275)。 所以，想请问各位大大，在maxLik()中怎么写呢？ 我看它的help,看不太懂

请问楼主解决了么 我也遇到这个问题了

楼主解决问题了么 求解答

hhdxlll 发表于 2022-10-28 16:31
楼主解决问题了么 求解答

111111111