Hi all, in the function DPdensity( ) example -- ?DPdensity --, we can extract the density estimate by: # Extracting the density estimate cbind(fit1.1$x1,fit1.1$dens) cbind(fit1.2$x1,fit1.2$dens) How is density calculated? I can imagine it is a mixture of normal densities, but where is the information about the density stored? Can we extract that so that we can write out the density function? My second question is: what is the last three columns of -- fit1.1$save.state$randsave---? The plot of the last column of fit1.1$save.state$randsave doesn't look like the original data? Thanks a lot. David [[alternative HTML version deleted]]