Need help NIntegrating stepwise 3D probability density functions
- To: mathgroup at smc.vnet.net
- Subject: [mg114071] Need help NIntegrating stepwise 3D probability density functions
- From: Jason Neuswanger <jrneuswanger at alaska.edu>
- Date: Mon, 22 Nov 2010 07:40:36 -0500 (EST)
I need to integrate 3D probability density functions over the region
in which the probability density exceeds some threshold (the end goal
being to calculate isopleths in 3D animal territories fit by kernel
distributions).
The simplest test case I could come up with is this:
testDistribution =
MultinormalDistribution[{1, 2,
4}, {{2, -1/4, 1/3}, {-1/4, 2/3, 1/5}, {1/3, 1/5, 1/2}}];
thresholdDensity = 0.05;
NIntegrate[
PDF[testDistribution, {x, y, z}]*
UnitStep[PDF[testDistribution, {x, y, z}] - thresholdDensity]
, {x, -Infinity, Infinity}, {y, -Infinity, Infinity}, {z, -Infinity,
Infinity}
, Method -> {"LocalAdaptive", "SymbolicProcessing" -> 0}
]
I've tried quite a few things and just can't get remotely sensible
results. The analogous case in 1 dimension works just fine, shown
below:
testDistribution2 = NormalDistribution[];
thresholdDensity2 = 0.35;
NIntegrate[
PDF[testDistribution2, {x}]*
UnitStep[PDF[testDistribution2, {x}] - thresholdDensity2]
, {x, -Infinity, Infinity}
, Method -> {"LocalAdaptive", "SymbolicProcessing" -> 0}
]
(* check answer using symbolic methods *)
thresholdRoot =
x /. FindRoot[
PDF[testDistribution2, {x}] - thresholdDensity2, {x, 1}];
Integrate[
PDF[testDistribution2, {x}], {x, -thresholdRoot, thresholdRoot}]
Any help here would be VERY much appreciated!