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

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

ただいまの
回答率

87.38%

scipy.optimize.fmin_bfgsでValue Errorが起きる

受付中

回答 1

投稿

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

score 10

前提・実現したいこと

CourseraのML(ex2)の課題をPythonで書いているのですが、コスト関数の最小値を求めるところでscipy.optimize.fmin_bfgsのモジュールを使いましたが、エラーが起きます。色々試してみましたが、手詰まりで助けていただきたいです...

発生している問題・エラーメッセージ

ValueError: shapes (3,) and (100,1) not aligned: 3 (dim 0) != 100 (dim 0)

行列掛け算の行と列がマッチしてないといったエラーだと思いますが, そこは確認しており, test_thetaを使ってコストを計算するとちゃんと計算できます。

該当のソースコード

import numpy as np 
import pandas as pd 
import matplotlib.pyplot as plt 

#Load data
data = pd.read_csv('ex2data1.txt', header=None)
X = np.array([data[0],data[1]]).T
y = np.array([data[2]]).T
m, n = np.shape(X)

#Add intercept term to X
X = np.hstack((np.ones((m, 1)), X))
#Initialize theta
initial_theta = np.zeros((n + 1, 1))

#Define sigmoid function
def sigmoid(z):
    g = 1 / (1 + np.exp(-z))
    return g

#Define cost function for logistic regression
def computeCost(X, y, theta):
    h = sigmoid(X @ theta)
    J = np.sum(-y * np.log(h) - (1-y) * np.log(1-h)) / m
    return J

import scipy.optimize as opt
result = opt.fmin_bfgs(computeCost, x0=initial_theta, args=(X, y))

試したこと

変数test_theta=[-24, 0.2, 0.2]としてコスト関数を計算させると、

test_theta = [-24, 0.2, 0.2]
cost = computeCost(X, y, test_theta)
print(cost)
# >> 223.5416509878586


と出るので、行列の形はミスってはいないと思うのですが...

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

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

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

    クリップを取り消します

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

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

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

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

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

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

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

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

    質問の評価を下げる

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

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

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

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

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

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

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

    詳細な説明はこちら

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

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

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

回答 1

0

エラーは何行目で起きているんでしょうか?
スタックトレースの全文をペーストしたほうが回答しやすそうです。

投稿

  • 回答の評価を上げる

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

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

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

  • 回答の評価を下げる

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

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

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

  • 2019/09/29 21:19

    ```
    Traceback (most recent call last):
    File "/Users/taro/Desktop/coursera/pythonByMyself/ex2/ex2.py", line 44, in <module>
    result = opt.fmin_bfgs(computeCost, x0=initial_theta, args=(X, y))
    File "/anaconda3/lib/python3.7/site-packages/scipy/optimize/optimize.py", line 916, in fmin_bfgs
    res = _minimize_bfgs(f, x0, args, fprime, callback=callback, **opts)
    File "/anaconda3/lib/python3.7/site-packages/scipy/optimize/optimize.py", line 970, in _minimize_bfgs
    gfk = myfprime(x0)
    File "/anaconda3/lib/python3.7/site-packages/scipy/optimize/optimize.py", line 300, in function_wrapper
    return function(*(wrapper_args + args))
    File "/anaconda3/lib/python3.7/site-packages/scipy/optimize/optimize.py", line 730, in approx_fprime
    return _approx_fprime_helper(xk, f, epsilon, args=args)
    File "/anaconda3/lib/python3.7/site-packages/scipy/optimize/optimize.py", line 664, in _approx_fprime_helper
    f0 = f(*((xk,) + args))
    File "/anaconda3/lib/python3.7/site-packages/scipy/optimize/optimize.py", line 300, in function_wrapper
    return function(*(wrapper_args + args))
    File "/Users/taro/Desktop/coursera/pythonByMyself/ex2/ex2.py", line 34, in computeCost
    h = sigmoid(X @ theta)
    ValueError: matmul: Input operand 1 has a mismatch in its core dimension 0, with gufunc signature (n?,k),(k,m?)->(n?,m?) (size 100 is different from 3)
    ```

    こんな感じです。

    キャンセル

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

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

関連した質問

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