質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

90.48%

  • Python 3.x

    6935questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。

Pythonで RuntimeWarning: divide by zero encountered in log -infのエラー

解決済

回答 4

投稿

  • 評価
  • クリップ 0
  • VIEW 461

gymgym

score 68

Pythonを用いてウィシャート分布を得たいと考えています

import pandas as pd
import numpy as np
import os
import math
import matplotlib.pyplot as plt
from scipy import arange, float64
from scipy.linalg import det
from scipy.special import digamma

def wishart(D, v_k, W):
    arr = float64(v_k + 1 - arange(1,D)) / 2
    return digamma(arr).sum() + D * np.log(2)+ np.log(det(W))

N=29
# 平均ベクトル
row_sum = np.sum(data, axis=0)
x_ave = row_sum/N

#分散
S = np.cov(data, rowvar=1, bias=1)

b0 = 0.1
m0 = 0
α = 0.1
W = np.eye(500)*α
v0 = 500

a1 = x_ave - m0
a2 = a1.T
b_k = b0 + N
m_k = 1/b_k*(b0*m0 + N*x_ave)
Wk_inverse = np.linalg.inv(W) + N*S + (b0*N/b0+N)* np.dot(a1,a2)


以上のコードで

 RuntimeWarning: divide by zero encountered in log  -infのエラー 


このようなエラーがでます

np.log(det(W))


この部分でlogに0が入っているみたいなのですが、0が入っているか確かめる方法はありますか。

また、0が入っているというエラーではなかったらどのようなことが原因でしょうか

よろしくお願い致します

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 4

+2

前以て行列の階数を調べておけばいいんじゃないでしょうか。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

checkベストアンサー

+1

det(W)が0になっているということですね。det(W)の結果が浮動小数点になると思うのでちょっと危なっかしいですが、if det(W) == 0:で捕まえられんじゃないでしょうか。

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

+1

上記のコードではwishartがコールさていませんが、関係ないものということでしょうか?

投稿

編集

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/05/16 18:37

    あら、ホントだ。

    キャンセル

  • 2018/05/17 09:56

    申し訳ございません。抜けてました。

    キャンセル

+1

... この部分でlogに0が入っているみたい ...

まずは、
print(W)
orint(det(W))
を追加して、値を確認してもみるとよいです。

W の内容が正しいものだとしたら、 det(W) が 0 のときは、 log をしないようにコードを修正すればよいです。
W のないようが想定しているものでなければ、W を計算してい部分をもなおしていく...
というように、確認・修正を繰返してみてください。

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2018/05/17 08:27

    orint

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 90.48%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる

関連した質問

  • 解決済

    cythonの型定義

    cythonファイルに型を定義したいです もともとのpyxファイルと 自分で型を定義したときのコードをのせます gamma1.pyx from math import exp

  • 解決済

    効率よく粒子の密度を求める方法

    与えられた多数の粒子の位置情報から粒子密度を求めたい。 粒子の位置情報を乱数で与え、そこからグリッドで分割された(一次元)空間における粒子の密度を推定するプログラムを作成しています

  • 受付中

    Python プログラム エラー

    学校の課題で 10から99までの数字をn回ランダムに出力したあと(条件:nは1000以上の数字)n個の数字の平均と分散、n個の数字の中での最大値、最小値を求めなさい。 という問題を

  • 解決済

    Pythonのtkinterのエラー

    2017-12-19 15:48:07.790 python[51777:1128073] *** Terminating app due to uncaught exceptio

  • 解決済

    高階関数の構造が理解できません。

    Python初心者です。書籍を購入して勉強していますが、内容を読んでも意味が分からない箇所があったので質問します。 書籍:Pythonプログラミング パーフェクトマスター (秀

  • 解決済

    関数について。 (コードの短縮)

     前提・実現したいこと 関数を使ったコードの短縮。 pythonの学習を始めたのですが書いたコードを関数を使って短くしたいと思ってもどうやればいいのか検討がつきません。どこで関数を

  • 解決済

    python エラー処理について

    エラー発生後 エラーが無くなるまで ほげほげ のコードを実行したいのですがどうすればいいでしょうか?宜しくお願い致します。 try: # tryの中に行いたい処理を書く

  • 解決済

    pythonで,ある条件のエラーを起こしたときに分岐する文を書くには

    タイトル通りです math.log()でlogの計算をしようとおもっているのですが 値が非常に小さいときに ValueError: math domain error のエラー

同じタグがついた質問を見る

  • Python 3.x

    6935questions

    Python 3はPythonプログラミング言語の最新バージョンであり、2008年12月3日にリリースされました。