MATLABでの分布出力

Question

ベイズ推論における1次元ガウス分布の精度が未知の時のパラメータの分布と予測分布の出力方法が分かりません。
また、多次元ガウスの場合の2つの分布の出力方法もどのように応用するか教えていただきたいです。
x = allsubj_eyedata(:,1);
x = rmmissing(x);
mu = mean(x);
%視点のX軸とY軸の2次元データを使う
N=length(x);

disp(N)

%prior hyper-parameters
a = 1;
b = 1;

%poster hyper-parameters
a_hat = N/2 + a;
b_hat = sum((x-mu).^2)/2 + b;

%pposterior distribution of lambda
%gampdfのbの定義とtextのbの定義が違うので1/b_hatに修正する必要有
x_axis = 0:0.01:15;
y_axis = gampdf(x_axis,a_hat,1/b_hat);

% plot distribution of x
figure
plot(x_axis,y_axis)

%predictive distribution
lambda = a_hat/b_hat;
nu = 2*a_hat;
x_axis = -10:0.01:10;
I =length(x_axis);
y_axis = zeros(1,I);
for i =1:I
    y_axis(1,i) = student_t_dist(x_axis(1,i),mu,lambda,nu);
end

figure
plot(plot(x_axis,y_axis))

function [y] = student_t_dist(x,mu,lambda,nu)
%y = gamma((nu+1)/2)/gamma(nu/2)*(lambda?(pi*nu))^(1/2)*(1+lambda/nu*(x-mu)^2)^(-(nu+1)/2);
G = (nu/(2*exp(1)*(nu+1)))^(1/2)*((nu+1)/nu)^(nu/2)*(nu+1)^(1/2);
y = G*(lambda/(pi*nu))
end

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

関数または変数 'student_t_dist' が未定義です。

### 
該当するコマンド
``` MATLAB
y_axis(1,i) = student_t_dist(x_axis(1,i),mu,lambda,nu);

Answer 1

student_t_distを別ファイルで定義する、または以下のように直接定義するとよいでしょう。

x = allsubj_eyedata(:,1);%<--------------- ??????????????????
x = rmmissing(x);
mu = mean(x);

student_t_dist = @(x,mu,lambda,nu) (nu/(2*exp(1)*(nu+1)))^(1/2)*((nu+1)/nu)^(nu/2)*(nu+1)^(1/2)*(lambda/(pi*nu));

%視点のX軸とY軸の2次元データを使う
N=length(x);

disp(N)

%prior hyper-parameters
a = 1;
b = 1;

%poster hyper-parameters
a_hat = N/2 + a;
b_hat = sum((x-mu).^2)/2 + b;

%pposterior distribution of lambda
%gampdfのbの定義とtextのbの定義が違うので1/b_hatに修正する必要有
x_axis = 0:0.01:15;
y_axis = gampdf(x_axis,a_hat,1/b_hat);

% plot distribution of x
figure
plot(x_axis,y_axis)

%predictive distribution
lambda = a_hat/b_hat;
nu = 2*a_hat;
x_axis = -10:0.01:10;
I =length(x_axis);
y_axis = zeros(1,I);
for i =1:I
    y_axis(1,i) = student_t_dist(x_axis(1,i),mu,lambda,nu);
end

figure
plot(x_axis,y_axis) %<--------------- plot(plot(x_axis,y_axis))を修正