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

Question

前提・実現したいこと

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

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

Answer 1

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

Comments

[TaroNoguchi]

```
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)
```

こんな感じです。