AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • Início
  • system&network
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • Início
  • system&network
    • Recentes
    • Highest score
    • tags
  • Ubuntu
    • Recentes
    • Highest score
    • tags
  • Unix
    • Recentes
    • tags
  • DBA
    • Recentes
    • tags
  • Computer
    • Recentes
    • tags
  • Coding
    • Recentes
    • tags
Início / coding / Perguntas / 79250609
Accepted
BlackPhoenix
BlackPhoenix
Asked: 2024-12-04 17:48:43 +0800 CST2024-12-04 17:48:43 +0800 CST 2024-12-04 17:48:43 +0800 CST

Média e mediana da distribuição dada pelo histograma numpy

  • 772

Suponha que você tenha um numpy histogramcalculado a partir de alguns dados (aos quais você não tem acesso), então você só conhece bins e contagens. Existe uma maneira eficiente de calcular a média e a mediana da distribuição descrita pelo histograma?

numpy
  • 1 1 respostas
  • 26 Views

1 respostas

  • Voted
  1. Best Answer
    mozway
    2024-12-04T18:14:18+08:002024-12-04T18:14:18+08:00

    Não, você não pode. Após agregação como histograma, a informação inicial é parcialmente perdida. Você não pode calcular exatamente a média/mediana da população original.

    Como demonstração, aqui estão duas matrizes diferentes (com médias/medianas diferentes) que fornecem as mesmas contagens e compartimentos:

    a1 = np.array([10, 20, 100, 300, 310])
    np.mean(a1), np.median(a1)
    # (148.0, 100.0)
    
    a2 = np.array([10, 10, 130, 300, 310])
    np.mean(a2), np.median(a2)
    # (152.0, 130.0)
    
    np.histogram(a1, bins=2)
    # (array([3, 2]), array([ 10., 160., 310.])
    
    np.histogram(a2, bins=2)
    # (array([3, 2]), array([ 10., 160., 310.])
    

    Aproximação

    No entanto, você pode determinar os limites da média:

    low = np.average(bins[:-1], weights=cnt)
    high = np.average(bins[1:], weights=cnt)
    print(f'The average is in the {low}-{high} range.')
    # The average is in the 70.0-220.0 range.
    

    E para a mediana:

    cnt_cumsum = np.add.accumulate(cnt)
    idx = np.searchsorted(cnt_cumsum, half)
    low = bins[idx]
    high = bins[idx+1]
    print(f'The median is in the {low}-{high} range.')
    # The median is in the 10.0-160.0 range.
    

    Exemplo com 1000 valores aleatórios e 20 caixas:

    True data mean: 0.496, median: 0.481
    The average is in the 0.471-0.521 range.
    The median is in the 0.45-0.5 range.
    
    • 2

relate perguntas

  • Como classificar o tensor em "lote" por valor de chave específico?

  • Aviso de descontinuação do notebook Jupyter ao encontrar a raiz do determinante de uma matriz

  • Como você concatena matrizes internas do tensor ao longo do eixo?

  • Digite regra de promoção para i4 e S8 no documento numpy

  • Transmitindo uma matriz numpy em uma matriz de tamanho maior usando uma matriz de índice

Sidebar

Stats

  • Perguntas 205573
  • respostas 270741
  • best respostas 135370
  • utilizador 68524
  • Highest score
  • respostas
  • Marko Smith

    Vue 3: Erro na criação "Identificador esperado, mas encontrado 'import'" [duplicado]

    • 1 respostas
  • Marko Smith

    Por que esse código Java simples e pequeno roda 30x mais rápido em todas as JVMs Graal, mas não em nenhuma JVM Oracle?

    • 1 respostas
  • Marko Smith

    Qual é o propósito de `enum class` com um tipo subjacente especificado, mas sem enumeradores?

    • 1 respostas
  • Marko Smith

    Como faço para corrigir um erro MODULE_NOT_FOUND para um módulo que não importei manualmente?

    • 6 respostas
  • Marko Smith

    `(expression, lvalue) = rvalue` é uma atribuição válida em C ou C++? Por que alguns compiladores aceitam/rejeitam isso?

    • 3 respostas
  • Marko Smith

    Quando devo usar um std::inplace_vector em vez de um std::vector?

    • 3 respostas
  • Marko Smith

    Um programa vazio que não faz nada em C++ precisa de um heap de 204 KB, mas não em C

    • 1 respostas
  • Marko Smith

    PowerBI atualmente quebrado com BigQuery: problema de driver Simba com atualização do Windows

    • 2 respostas
  • Marko Smith

    AdMob: MobileAds.initialize() - "java.lang.Integer não pode ser convertido em java.lang.String" para alguns dispositivos

    • 1 respostas
  • Marko Smith

    Estou tentando fazer o jogo pacman usando apenas o módulo Turtle Random e Math

    • 1 respostas
  • Martin Hope
    Aleksandr Dubinsky Por que a correspondência de padrões com o switch no InetAddress falha com 'não cobre todos os valores de entrada possíveis'? 2024-12-23 06:56:21 +0800 CST
  • Martin Hope
    Phillip Borge Por que esse código Java simples e pequeno roda 30x mais rápido em todas as JVMs Graal, mas não em nenhuma JVM Oracle? 2024-12-12 20:46:46 +0800 CST
  • Martin Hope
    Oodini Qual é o propósito de `enum class` com um tipo subjacente especificado, mas sem enumeradores? 2024-12-12 06:27:11 +0800 CST
  • Martin Hope
    sleeptightAnsiC `(expression, lvalue) = rvalue` é uma atribuição válida em C ou C++? Por que alguns compiladores aceitam/rejeitam isso? 2024-11-09 07:18:53 +0800 CST
  • Martin Hope
    The Mad Gamer Quando devo usar um std::inplace_vector em vez de um std::vector? 2024-10-29 23:01:00 +0800 CST
  • Martin Hope
    Chad Feller O ponto e vírgula agora é opcional em condicionais bash com [[ .. ]] na versão 5.2? 2024-10-21 05:50:33 +0800 CST
  • Martin Hope
    Wrench Por que um traço duplo (--) faz com que esta cláusula MariaDB seja avaliada como verdadeira? 2024-05-05 13:37:20 +0800 CST
  • Martin Hope
    Waket Zheng Por que `dict(id=1, **{'id': 2})` às vezes gera `KeyError: 'id'` em vez de um TypeError? 2024-05-04 14:19:19 +0800 CST
  • Martin Hope
    user924 AdMob: MobileAds.initialize() - "java.lang.Integer não pode ser convertido em java.lang.String" para alguns dispositivos 2024-03-20 03:12:31 +0800 CST
  • Martin Hope
    MarkB Por que o GCC gera código que executa condicionalmente uma implementação SIMD? 2024-02-17 06:17:14 +0800 CST

Hot tag

python javascript c++ c# java typescript sql reactjs html

Explore

  • Início
  • Perguntas
    • Recentes
    • Highest score
  • tag
  • help

Footer

AskOverflow.Dev

About Us

  • About Us
  • Contact Us

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve