From: <HIDDEN>
Newsgroups: comp.soft-sys.matlab
Subject: Re: How to combine two probability density functions (PDFs)?
Date: Fri, 27 Aug 2010 21:02:23 +0000 (UTC)
Organization: US Geological Survey
Lines: 44
Message-ID: <i5994v$koj$>
References: <i516ij$opj$> <i51dlr$5bf$> <i526t0$ruu$>
Reply-To: <HIDDEN>
Content-Type: text/plain; charset=UTF-8; format=flowed
Content-Transfer-Encoding: 8bit
X-Trace: 1282942943 21267 (27 Aug 2010 21:02:23 GMT)
NNTP-Posting-Date: Fri, 27 Aug 2010 21:02:23 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 2442676
Xref: comp.soft-sys.matlab:665621

"Ross W" <> wrote in message <i526t0$ruu$>...
> "Roger Stafford" <> wrote in message <i51dlr$5bf$>...
> > "Ryan " <> wrote in message <i516ij$opj$>...
> > > I am writing to seek help with modeling the probability of an event (age), which is bracketed by two ages, where each age is defined by a Gaussian probability density function (PDF).
> > > 
> > > To start simply, consider the case where an event is known to have occurred between two brackets. For example:
> > > -The upper age bracket is described by a Gaussian age with a mean value of 5000 (yr) and a 1-sigma value of 200 (yr).
> > > -The lower age bracket is also described by a Gaussian age, but with a mean value of 2000 (yr) and a 1-sigma value of 400 (yr).
> > > 
> > > My 1st question: how can I use the two age brackets to define a new PDF reflecting the probability of the event, which is known to have occurred between the two features? Note: my intuition is that the shape of the PDF should be a plateau between the brackets with Gaussian tails.
> > > 
> > > My 2nd question: how does the solution change where the distance between the end-member brackets decreases and the sigma values increase (e.g., when the tails of each age bracket overlap, especially given that the event defined by the old bracket has to have occurred prior to the event defined by the younger bracket)
> > > 
> > > Thanks for any help that you can provide!
> > > 
> > > Sincerely,
> > > Ryan
> > - - - - - - - - -
> >   I am struggling to understand what you are asking, Ryan.  It doesn't seem like a well-defined question.  You don't appear to have given the "event" (between the brackets) any a-priori probability distribution and therefore no question about its conditional probability given that it lies between the two Gaussians can be answered.  Also you haven't stated whether these Gaussians are mutually independent or at least jointly Gaussian.  I think you need to describe your problem a lot more carefully before anyone can give you meaningful help.
> > 
> > Roger Stafford
> Hi Ryan,
> Is it true that the age you seek is UNIFORMLY distributed between the upper age and the lower age? Or is it distributed in some other way between those two ages?
> Ross

Hi Ross,

Thanks for your response.

My expectation is that the age I seek is UNIFORMLY distributed.

Here's a description of a scenario: There is a stratigraphic horizon which I want to date (e.g., represent an earthquake in the geologic record). Below that horizon, I have dated deeper layer of sediment that provides a maximum constraint on the age of the earthquake. Similarly, above the earthquake horizon, I have dated a shallower layer of sediment that provides a minimum constraint on the age of the earthquake. The question: how do I describe the age of the earthquake based on these two boundaries?

Here is a modified scenario, from the one I previously described: Maximum age constraint = 10000 (yr) with a 1-sigma value of 400 (yr) [normally distributed] / Minimum age constraint = 5000 (yr) and a 1-sigma value of 200 (yr) [normally distributed].

My intuition is that I lack sufficient knowledge about sedimentation rates, erosion rates, etc., to state with confidence where the earthquake event that occurs between the boundaries is closer to the maximum or minimum bracket. Thus, I expect to see a uniformly distributed zone in the region between the constraints (e.g., flat PDF from 6000-8500 yr, in the example above), with tails that taper toward the boundaries.

Many thanks for additional input!