一連の観測値を、それらの平均と共分散とともにgnuplotに表示したいと思います。誰かが私がこのようなことをする方法を知っていますか?これをうまく行う別のプログラムはありますか?Octaveのようなプログラムはこれを行うことができますか?統計を取得する方法はありますが、それらの統計を表示する方法については何も見つかりませんでした。
明確にするために、データポイント(2dまたは3d)を表示し、データセットにプロットされたそれらの平均を表示し、共分散を表す楕円を表示したいと思います。
はい、Octaveでこれを行うことができます。
## create 2x1000 matrix of random data with normal distribution
data = randn (1000, 2);
## plot data (the o is for dots)
plot (data(:,1), data(:,2), "o");
## get mean from each
mu = mean (data);
## calculate covariance matrix
R = cov (data);
## calculate the ellipse points
A = chol (R, "lower");
theta = linspace (0, 2*pi, 1000);
x = mu' + 2.5 .* A * [cos(theta); sin(theta)];
## plot the ellipse
hold on;
plot(x(1,:), x(2,:), "r", "LineWidth", 2);
オクターブ4.0.2でエラーが発生し、carandraugの回答の「x=」行を次のように変更しました。
x = mu' + 2.5 .* A * [cos(theta); sin(theta)];