Displaying 1 result from an estimated 1 matches for "theldastar".
2010 Sep 21
2
Need help for EM algorithm ASAP !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
...library(QRMlib)
library(ghyp)
############ simulation part
simulation<-function(n,lambda,mu,thelda,gamma,sigma,beta){
set.seed(235)
chi<-thelda^2
psi<-gamma^2
W <- rGIG(n, lambda, chi, psi);
Z <- rnorm(n,0,1);
y<-mu + beta * W + sqrt(W) * Z *gamma;
for (i in 1:n){
theldastar<-rep(0,n)
zi<-rep(0,n)
ti<-rep(0,n)
muthelda<-mu
gammathelda<-thelda*gamma
sigmathelda<-(thelda^2)*sigma
betathelda<-(thelda^2)*sigma*beta
lambdastar<-lambda-0.5
theldastar[i]<-sqrt(1+((y[i]-muthelda)/sigmathelda)^2)
gammastar<-sqrt((gammathelda^2)+((betathelda/...