search for: klambda

...t;-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/sigmathelda)^2)) klambda1<-besselM3(lambdastar+1, x=2, logvalue=FALSE) klambda<-besselM3(lambdastar,x=2,logvalue=FALSE) klambda2<-besselM3(lambdastar-1,x=2,logvalue=FALSE) zi[i]<-((theldastar[i]*klambda1*(theldastar[i]*gammastar))/(gammastar*klambda*theldastar[i]*gammastar)) ti[i]<-((gammastar*klambda2*(...