param<-function(pars.est,dados){ SIGMA <- montaSigma(pars.est, dados$d) tchSIGMA <- t(chol(SIGMA)) ONE <- kronecker(diag(1, 2), rep(1, dados$n1)) ONEchSIGMA <- solve(tchSIGMA, ONE) mus <- solve(crossprod(ONEchSIGMA), crossprod(ONEchSIGMA, solve(tchSIGMA, dados$X))) dif<-dados$X-ONE%*%mus sig1.<-(1/(dados$n1+dados$n2))*(t(dif)%*%solve(tchSIGMA,dif)) sig1<-sig1.^2 sig2<-(pars.est$sig2)^2 tau1<-(pars.est$tau1)^2 tau2<-(pars.est$tau2)^2 nu1<-pars.est$nu1 nu2<-pars.est$nu2 a<-pars.est$a rho<-pars.est$rho mu1<-mus[1,] mu2<-mus[2,] parametros <- c(mu1,mu2,sig1,sig2,tau1,tau2,nu1,nu2,a,rho) names(parametros) = c('mu1','mu2','sig1','sig2','tau1','tau2','nu1','nu2','a','rho') return(parametros) }