% This subroutine accomplishes reweighting of posterior simulator output % and/or construction of new functions of interest. The code given here % leaves the prior and functions of interest unchanged, and therefore % accomplishes nothing. This code must be replaced by code specific to % the client's interests. % Inputs: % npar Number of functions of interest in the investigator's posterior % simulation file. % thetai Vector of dimension npar, containing the parameter values read from % the investigator's simulation file in the current record. % pri Logarithm of the investigator's prior evaluated at thetai % Outputs: % nc Number of functions of interest in the client's posterior simulation % file. % gc Vector of dimension nc, containing the values of the client's % function of interest % priorc Logarithm of the client's prior evaluated at thetai function [nc, gc, priorc]=client(npar,thetai,pri); nc=npar; priorc=pri; gc=thetai(1:nc); end;