2007年4月9日月曜日

データの不偏分散,不偏標準偏差

> #統計R解析の場合
> data <- runif(10000) #一様乱数を発生させる関数( 0.0 - 1.0 ) > d_var <- var(data) # 不偏分散 > d_var
[1] 0.08373053
> d_sd <- sd(data) # 不偏標準偏差 > d_sd
[1] 0.2893623
>

%MATLABの場合
echo on
N = 10000;
data = rand(N,1); % 一様分布の擬似乱数
d_var = var(data); % 不偏分散
disp(d_var);
0.0842
d_std = std(data); % 不偏標準偏差
disp(d_std);
0.2901

0 件のコメント: