This is the part of the function in my likelihood function
prob<-function(t1,t2,m1,m2,s_r,s_p,s_e,cor_m){
if ((t1==0) & (t2==0))
log_lik<-log(pmvnorm(lower=rep(-Inf, 2),
upper=c(m1/sqrt(s_r+s_p+s_e+1),m2/sqrt(s_r+s_p+s_e+1)),
mean=rep(0,2), corr=cor_m))
else if ((t1==0) & (t2==1))
log_lik<-log(pmvnorm(lower=c(-Inf,m2/sqrt(s_r+s_p+s_e+1)),
upper=c(m1/sqrt(s_r+s_p+s_e+1),Inf),
mean=rep(0,2), corr=cor_m))
else if ((t1==1) & (t2==0))
log_lik<-log(pmvnorm(lower=c(m1/sqrt(s_r+s_p+s_e+1),-Inf),
upper=c(Inf,m2/sqrt(s_r+s_p+s_e+1)),
mean=rep(0,2), corr=cor_m))
else if ((t1==1) & (t2==1))
log_lik<-log(pmvnorm(lower=c(m1/sqrt(s_r+s_p+s_e+1),m2/sqrt(s_r+s_p+s_e+1)),
upper=c(Inf,Inf),
mean=rep(0,2), corr=cor_m))
return(log_lik)
}
--
View this message in context:
http://r.789695.n4.nabble.com/Maximization-problem-in-the-optim-function-tp4473408p4473602.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
