Mercurial > hg > ltpda
view m-toolbox/classes/+utils/@math/fitPrior.m @ 22:b11e88004fca database-connection-manager
Update collection.fromRepository
author | Daniele Nicolodi <nicolodi@science.unitn.it> |
---|---|
date | Mon, 05 Dec 2011 16:20:06 +0100 |
parents | f0afece42f48 |
children |
line wrap: on
line source
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Function fitpriors fits the pdf computed from a pest object. % % It returns matrix (# of params x 3) with the mean, sigma and % % normalization constant for each parameter at each column. -Nikos- % %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% function params = fitPrior(prior,nparam,chain,bins) params = []; for ii=1:nparam x = linspace(min(prior(:,2*ii-1)),max(prior(:,2*ii-1)),bins); y = normpdf(x,mean(chain(:,ii)),std(chain(:,ii))); s=sum(y); params(ii,:) = [mean(chain(:,ii)) std(chain(:,ii)) s]; end end