A method of analyzing controlled source electromagnetic (CSEM) survey data to determine probability density functions (PDFs) for values of an electromagnetic parameter at locations in a subterranean region of interest is provided. Structural features in the subterranean strata are identified, e.g. from seismic survey data. An initial PDF for values of the electromagnetic parameter is then assigned to each feature. Models specifying values for the parameter in each structural feature are generated by sampling the PDFs. A subset of the models are deemed acceptable based on an acceptance criterion. The PDF for each feature is modified based on values for the parameter in the subset of accepted models to generate replacement PDFs for each feature. The process may be iterated a number of times to generate final PDFs for values of the electromagnetic parameter in the structural features identified in the subterranean region of interest.