Displaying 2 results from an estimated 2 matches for "nburn".
Did you mean:
burn
2009 Feb 21
0
density estimation for d>2 for the DPpackage
...2)
#with the following priors
# Prior information
s2<-matrix(c(2,0,0,0,2,0,0,0,2),ncol=3)
m2<-c(5,10,15)
psiinv2<-solve(matrix(c(10000,0,0,0,1,0,0,0,2),ncol=3))
prior<-list(a0=1,b0=1/5,nu1=4,nu2=4,s2=s2,
m2=m2,psiinv2=psiinv2,tau1=0.01,tau2=0.01)
state <- NULL
# MCMC parameters
nburn<-1000
nsave<-5000
nskip<-10
ndisplay<-1000
mcmc <- list(nburn=nburn,nsave=nsave,nskip=nskip,ndisplay=ndisplay)
# Fit the model
fit3<-DPdensity(y=rnormm,mcmc=mcmc, prior=prior,
state=state,status=TRUE,na.action=na.omit)
All is fine untill I bacame...
MCMC scan 1000 of 5000 (CPU...
2009 Mar 31
1
error during DPpackage compilation
...*(p+1)/2+3)
? ? ? ? ? ? ? ? ? ? ? ? ? 2
Invalid declaration of or reference to symbol `nsave' at (2) [initially seen
at (1)]
LDDPdensity.f:12:
? ? ? ? & ? ? ? ? ? ? ? ? ? ? ? ?mcmc,nsave,seed,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 1
LDDPdensity.f:236: (continued):
? ? ? ? ?integer mcmc(3),nburn,nskip,nsave,ndisplay
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?2
Invalid declaration of or reference to symbol `nsave' at (2) [initially seen
at (1)]
LDDPdensity.f:9:
? ? ? ? & ? ? ? ? ? ? ? ? ? ? ? ?cpo,thetasave,randsave,
? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?1
LDDPdensity.f:291: (continued):...