This is part of a code of mine.I plot the density of beta.draws nad overlay the density of another model using lines()
You should be careful regarding the plotting regions,though!Code:plot(density(fix.m7.1n.rj$beta.draws[i,]),ylim=c(0,k),main="",col="steelblue") lines(density(fix.m7.1n.mc3$beta.draws[i,]),col="red");grid()





Reply With Quote


