Dear List, I am trying to estimate a 3 dimensional density through the DPpackage. For example # model sigma <- matrix(c(0.1,0.05,0.05,0.05,0.1,0.05,0.05,0.05,0.1), ncol=3) rnormm<- rmvnorm(n=100, mean=c(5,100,150), sigma=sigma) sigma2 <- matrix(c(10,0.05,0.05,0.05,10,0.05,0.05,0.05,10), ncol=3) rnormm2<- rmvnorm(n=100, mean=c(20,1,110), sigma=sigma) rnormm<-rbind(rnormm,rnormm2) #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 time: 32.125 s) MCMC scan 2000 of 5000 (CPU time: 61.797 s) MCMC scan 3000 of 5000 (CPU time: 91.281 s) MCMC scan 4000 of 5000 (CPU time: 120.750 s) MCMC scan 5000 of 5000 (CPU time: 150.047 s) Error in matrix(foo$f, nrow = ngrid, ncol = ngrid) : Versuch ein Attribut von NULL zu setzen.... Some idea? Best regards, rodrigo -- Besuchen Sie meine Web Seite rodrigoherrera.blogspot.com Viel Spass beim Surfen ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ Habe nun, ach! Philosophie, Juristerei, Theologie und Medizin und leider auch Volkswirtschaft durchaus studiert, mit hei?em Bem?hn. Da steh' ich nun, ich armer Tor, und bin so klug als wie zuvor! Hei?e Magister, hei?e Doktor gar, und ziehe schon and die zehen Jahr' herauf, herab und quer und krumm meinen Sch?ler an der Nase herum - Und sehe, da? wir nichts wissen k?nnen! Das will mir schier das Herz verbrennen. Aus dem Faust Johann Wolfgang von Goethe 1749 - 1832 Jetzt 1 Monat kostenlos! GMX FreeDSL - Telefonanschluss + DSL f?r nur 17,95 Euro/mtl.!* dsl.gmx.de/?ac=OM.AD.PD003K11308T4569a