# How determine The most possibility probable with programming

i have a scatter chart of probable of happen some thing like below chart

i what determine the probable The most point in below chart , can use a mathematic method to determine for programming

=================

Do you just have this chart as a picture, or do you also have the values for the data points?
– Karsten 7.
Aug 29 ’14 at 8:59

have valuse of points , you can see the caluse of google chart link goo.gl/wY4OXg
– MOB
Aug 29 ’14 at 9:00

want single value of most likely value
– MOB
Aug 29 ’14 at 10:21

of course it’s better find some single value that are most likely values
– MOB
Aug 29 ’14 at 10:24

=================

1

=================

With

xValues={-1.27,-3.65,-0.64,8.65,-1.52,-1.65,-3.9,-2.21,-0.07,-1.62,1.5,-2.71,3.37,-8.29,16.16,-31.27,3.07,2.97,9,-2.63,16.03,-24.49,-0.35,-13.36,-2.57,5.08,4.53,-4.47,0.34,1.66,2.32,0.23,-1.77,-0.86,5.46,-0.84,-1.6,1.54,-3.69,-5.87,-2.09,-10.86,-1.33,0.64,0.84,3.96,1.57,2.07,5.87,-2.09,-1.57,-3.65,1.67,-10.23,-1.94,-0.38,1.45,-0.85,1.87,-3.71,0.4,-2.26,-1.66,0.77,1.9,-0.56,7.53,31.38,0.23,0.06,-5.05,9.04}
yValues={2.18,1.14,1.37,2.18,2.16,1.27,2.2,1.78,1.81,0.86,1.32,1.67,1.63,1.55,1.57,2,2.14,1.45,1.7,2.09,2.11,2.09,2.06,2.08,1.79,1.95,1.93,2.03,2.11,1.04,2.12,2.19,1.32,1.88,2.15,2.1,2.16,1.97,1.47,2.02,1.4,1.65,2.15,1.57,2.14,2.12,2.14,2.15,2.06,1.98,2.18,2.14,1.67,1.82,1.71,2.18,2.15,1.93,2.01,1.75,2.11,1.86,2.16,2.11,1.91,2.2,1.79,2.16,2.08,1.82,2.07,2.14}

There a two possible interpretations for these data points.

Interpretation as univariate distribution

Interpreting the yValues as weights and finding the most probable xValue as the objective.
One can create a WeightedData object using

wData = WeightedData[xValues, yValues]

and then calculate the mean with

Mean[wData]

-0.153426

and other statistics, like the standard deviation with

StandardDeviation[wData]

It is even possible to plot a PDF for the data

With[{dist = SmoothKernelDistribution[wData]},
Plot[Evaluate[PDF[dist, \[FormalT]]], {\[FormalT],
Min[wData[“InputData”]] – 2 dist[“Bandwidth”],
Max[wData[“InputData”]] + 2 dist[“Bandwidth”]}, PlotRange -> All]]

Interpretation as multivariate distribution

Interpreting the xValues and yValues as 2D data points and with “finding the most probable {x,y} data point” as the objective.
One can create an EmpiricalDistribution object

ed2D = EmpiricalDistribution[Transpose[{xValues, yValues}]]

and then calculate its mean with

Mean[ed2D]

{-0.227361, 1.88819}

A visualization of the PDF together with a red line at the mean point and a blue line at the maximum point

Module[{dist = SmoothKernelDistribution[Transpose[{xValues, yValues}]],
meanPoint = Mean[ed2D], maxPoint},
maxPoint = ArgMax[PDF[dist, {x, y}], {x, y}];
Show[Plot3D[Evaluate[PDF[dist, {x, y}]],
{x, Min[xValues] – 2 First@dist[“Bandwidth”], Max[xValues] + 2 First@dist[“Bandwidth”]},
{y, Min[yValues] – 2 Last@dist[“Bandwidth”], Max[yValues] + 2 Last@dist[“Bandwidth”]},
PlotRange -> All],
Graphics3D[{Darker@Red, Thickness[.02],
Line[{Append[meanPoint, 0], Append[meanPoint, 0.13]}]}],
Graphics3D[{Darker@Blue, Thickness[.02],
Line[{Append[maxPoint, 0], Append[maxPoint, 0.18]}]}] ]]

– MOB
Aug 29 ’14 at 11:10

Just curious:Why are you interpreting the y values as weights and not the whole thing as 2D points?
– Dr. belisarius
Aug 29 ’14 at 11:22

@belisarius It’s just that I “saw” the PDF of a univariate distribution when I looked at the image. Nevertheless, a multivariate distribution is reasonable too, based on the limited information given about the context of these data points. E.g. something like Plot3D[Evaluate@CDF[EmpiricalDistribution[Transpose[{xValues‌​,yValues}]],{x,y}],{‌​x,-35,35},{y,0,2.5},‌​Exclusions->None] makes sense too.
– Karsten 7.
Aug 29 ’14 at 11:55

@belisarius And of course Mean[EmpiricalDistribution[Transpose[{xValues, yValues}]]]
– Karsten 7.
Aug 29 ’14 at 12:08

1

Perhaps you could add that “another” POV
– Dr. belisarius
Aug 29 ’14 at 12:10