Python: Binning data with a weight










0















I have got a dataset with strong noise, e.g.



import numpy as np
import matplotlib.pyplot as plt

N = 1000
x = np.linspace(0,10,N)
y = x + 20 * np.random.rand(N)


I want to bin the data for a given binsize (or binnumber). By that I basically just mean a Δx. The binned data should be weighted by a gaussian function that you can think of as a gaussian that is extended over the y-axis weighting the data depending on the distance of the expectation value µ. Also, I want the data to give me the 1σ-error.



I know about numpy.digitize and scipy.stats.binned_statistic but I am failing to apply any of the two to get my desired binning. Maybe the latter should be the easiest to use for this case as it offers the parameter statistic=<function> but I am open to suggestions.










share|improve this question

















  • 1





    Did you take a look at numpy.histogram? It supports binning with weights.

    – MisterMiyagi
    Nov 13 '18 at 17:28











  • Yes but I don't quite get how this would look like... I need a specific example to start off but unfortunately I cannot find (m)any

    – famfop
    Nov 13 '18 at 17:34















0















I have got a dataset with strong noise, e.g.



import numpy as np
import matplotlib.pyplot as plt

N = 1000
x = np.linspace(0,10,N)
y = x + 20 * np.random.rand(N)


I want to bin the data for a given binsize (or binnumber). By that I basically just mean a Δx. The binned data should be weighted by a gaussian function that you can think of as a gaussian that is extended over the y-axis weighting the data depending on the distance of the expectation value µ. Also, I want the data to give me the 1σ-error.



I know about numpy.digitize and scipy.stats.binned_statistic but I am failing to apply any of the two to get my desired binning. Maybe the latter should be the easiest to use for this case as it offers the parameter statistic=<function> but I am open to suggestions.










share|improve this question

















  • 1





    Did you take a look at numpy.histogram? It supports binning with weights.

    – MisterMiyagi
    Nov 13 '18 at 17:28











  • Yes but I don't quite get how this would look like... I need a specific example to start off but unfortunately I cannot find (m)any

    – famfop
    Nov 13 '18 at 17:34













0












0








0








I have got a dataset with strong noise, e.g.



import numpy as np
import matplotlib.pyplot as plt

N = 1000
x = np.linspace(0,10,N)
y = x + 20 * np.random.rand(N)


I want to bin the data for a given binsize (or binnumber). By that I basically just mean a Δx. The binned data should be weighted by a gaussian function that you can think of as a gaussian that is extended over the y-axis weighting the data depending on the distance of the expectation value µ. Also, I want the data to give me the 1σ-error.



I know about numpy.digitize and scipy.stats.binned_statistic but I am failing to apply any of the two to get my desired binning. Maybe the latter should be the easiest to use for this case as it offers the parameter statistic=<function> but I am open to suggestions.










share|improve this question














I have got a dataset with strong noise, e.g.



import numpy as np
import matplotlib.pyplot as plt

N = 1000
x = np.linspace(0,10,N)
y = x + 20 * np.random.rand(N)


I want to bin the data for a given binsize (or binnumber). By that I basically just mean a Δx. The binned data should be weighted by a gaussian function that you can think of as a gaussian that is extended over the y-axis weighting the data depending on the distance of the expectation value µ. Also, I want the data to give me the 1σ-error.



I know about numpy.digitize and scipy.stats.binned_statistic but I am failing to apply any of the two to get my desired binning. Maybe the latter should be the easiest to use for this case as it offers the parameter statistic=<function> but I am open to suggestions.







python binning weighting






share|improve this question













share|improve this question











share|improve this question




share|improve this question










asked Nov 13 '18 at 17:24









famfopfamfop

28118




28118







  • 1





    Did you take a look at numpy.histogram? It supports binning with weights.

    – MisterMiyagi
    Nov 13 '18 at 17:28











  • Yes but I don't quite get how this would look like... I need a specific example to start off but unfortunately I cannot find (m)any

    – famfop
    Nov 13 '18 at 17:34












  • 1





    Did you take a look at numpy.histogram? It supports binning with weights.

    – MisterMiyagi
    Nov 13 '18 at 17:28











  • Yes but I don't quite get how this would look like... I need a specific example to start off but unfortunately I cannot find (m)any

    – famfop
    Nov 13 '18 at 17:34







1




1





Did you take a look at numpy.histogram? It supports binning with weights.

– MisterMiyagi
Nov 13 '18 at 17:28





Did you take a look at numpy.histogram? It supports binning with weights.

– MisterMiyagi
Nov 13 '18 at 17:28













Yes but I don't quite get how this would look like... I need a specific example to start off but unfortunately I cannot find (m)any

– famfop
Nov 13 '18 at 17:34





Yes but I don't quite get how this would look like... I need a specific example to start off but unfortunately I cannot find (m)any

– famfop
Nov 13 '18 at 17:34












0






active

oldest

votes











Your Answer






StackExchange.ifUsing("editor", function ()
StackExchange.using("externalEditor", function ()
StackExchange.using("snippets", function ()
StackExchange.snippets.init();
);
);
, "code-snippets");

StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "1"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);

StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);

else
createEditor();

);

function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: true,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: 10,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);



);













draft saved

draft discarded


















StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53286468%2fpython-binning-data-with-a-weight%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown

























0






active

oldest

votes








0






active

oldest

votes









active

oldest

votes






active

oldest

votes















draft saved

draft discarded
















































Thanks for contributing an answer to Stack Overflow!


  • Please be sure to answer the question. Provide details and share your research!

But avoid


  • Asking for help, clarification, or responding to other answers.

  • Making statements based on opinion; back them up with references or personal experience.

To learn more, see our tips on writing great answers.




draft saved


draft discarded














StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fstackoverflow.com%2fquestions%2f53286468%2fpython-binning-data-with-a-weight%23new-answer', 'question_page');

);

Post as a guest















Required, but never shown





















































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown

































Required, but never shown














Required, but never shown












Required, but never shown







Required, but never shown







這個網誌中的熱門文章

How to read a connectionString WITH PROVIDER in .NET Core?

In R, how to develop a multiplot heatmap.2 figure showing key labels successfully

Museum of Modern and Contemporary Art of Trento and Rovereto