Displaying 3 results from an estimated 3 matches for "betac".
Did you mean:
beta
2007 May 28
0
How to use filter function for this?
Hello,
I have been trying to compute the following:
for(t in 3:n){
h[t]<-(omegac/(1-betac))+alphac*Y[t-1]
b[t]<-Y[t-1]
bb[t]<-0
for(j in 2:(t-1)){
h[t]<-h[t]+alphac*((betac)^(j-1))*Y[t-j]
b[t]<-b[t]+((betac)^(j-1))*Y[t-j]
bb[t]<-bb[t]+((betac)^(j-2))*Y[t-j]*(j-1)
}
}
but it takes a long time. I have heard that "filter" fu...
2005 May 19
1
R 2.1.0 RH Linux Built from Source Segmentation Fault
...idates
#
oldbeta=beta0
mhess=diag(nvar)
candcov=chol2inv(chol(mhess))
root=s*chol(candcov)
priorcov=chol2inv(chol(A))
rootp=chol(priorcov)
rootpi=backsolve(rootp,diag(nvar))
#
# start main iteration loop
#
oldlpost=llmnl(y,X,beta0)+lmvn(beta0,betabar,rootpi)
naccept=0
for (rep in 1:R)
{
betac=oldbeta+t(root)%*%rnorm(nvar)
clpost=llmnl(y,X,betac)+lmvn(betac,betabar,rootpi)
ldiff=clpost-oldlpost
alpha=min(1,exp(ldiff))
if(alpha < 1) {unif=runif(1)} else {unif=0}
if (unif <= alpha)
{ beta0=betac
oldlpost=clpost
naccept=naccept+1}
oldbe...
2001 May 07
2
postscript problem
hello all,
when trying to produce a plot with the following code
postscript("g6-5000-0705-b.ps")
par(mfrow=c(4,1))
for (i in 1:4) {
somc<-cumsum(betaC[,i])/c(1:N)
plot(somc,ylim=c(min(somc),max(somc)));
}
dev.off()
I get a postscript file g6-5000-0705-b.ps, but if I try to open it
with ghostscript, or to convert it to pdf, I get the following error :
Error: /invalidfont in findfont
Operand stack:
Helvetica Font Helvetica 236270 Hel...