Hi, I am doing Bayesian Calibration for a model which has 24 parameters using MCMC. I get a sample of 100,000 points from the posterior pdf. Let's say, It is a dataframe named as "BC_output", with 24 parameters (V1,V2,....V24)
I would like to show the Maximum a Posteriori (MAP)parameter value of the model. Can I do it like this: For V1: density(BC_output$V1)$x[which.max(density(BC_output$V1)$y)] If not, what is the correct function to estimate a MAP value in R? Thanks a lot, -- _________________________________________________ Mr.Bin You Abt.Biometrie und Informatik Wonnhaldestraße 4 79100 Freiburg Telefon: 0761 / 4018-365 Fax: 0761 / 4018-333 bin....@forst.bwl.de www.fva-bw.de www.manfredproject.eu www.waldwissen.net ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
