以下を実行し、
import numpy as np
import matplotlib.pyplot as plt
from scipy import stats, integrate
from scipy.optimize import minimize_scalar
%precision 3
%matplotlib inline
x_range=[0,2]
y_range=[0,1]
def f_xy(x,y):
if 0<=y<=1 and 0<=x-y<=1:
return 4*y*(x-y)
else:
return 0
XY=[x_range,y_range,f_xy]
from functools import partial
# 周辺密度関数を定義
def f_X(x):
return integrate.quad(partial(f_xy, x), -np.inf, np.inf)[0]
X = [x_range, f_X]
def f_Y(y):
return integrate.quad(partial(f_xy, y), -np.inf, np.inf)[0]
Y = [y_range, f_Y]
xs=np.linspace(*x_range,100)
ys=np.linspace(*y_range,100)
fig=plt.figure(figsize=(12,4))
ax1=fig.add_subplot(121)
ax2=fig.add_subplot(122)
ax1.plot(xs,[f_X(x)for x in xs],color='gray')
ax2.plot(ys,[f_Y(y)for y in ys],color='gray')
ax1.set_title('Xの周辺密度関数',fontname='MS Gothic')
ax2.set_title('Yの周辺密度関数',fontname='MS Gothic')
plt.show()
x_range
やy_range
が定義されていないようですが、質問時の転記ミスでしょうか? そうでなければ以前に実行されていた内容が残って影響している可能性があるので、いったんjupyter notebookをリセットか終了・再起動してから、それらを定義して実行してみてください。いずれにしろ質問にそれらを追記してください。