相関係数

統計

概要

 閲覧数:1909  投稿日:2014-05-14  更新日:2014-05-14  

相関係数


コード

 function corrcoef($var_x, $var_y)

    {

        $mean_x = NULL;
        $mean_y = NULL;
        $sum_sq_x = NULL;
        $sum_sq_y = NULL;
        $sum_coproduct = NULL;

        if (!count($var_x) || !is_array($var_x))    die('配列じゃない・・・');

        if (!count($var_y) || !is_array($var_y))    die('配列じゃない・・・');

        

        $c = (count($var_x) > count($var_y)) ? count($var_x) : count($var_y);

        

        $c = 0;

        while (true) {

//Notice: Undefined offsetとなるので変更
//          if (is_null($var_x[$c]) || is_null($var_y[$c]))    break;
            if (!isset($var_x[$c]) || !isset($var_y[$c]))    break;
            $_var_x[] = $var_x[$c];

            $_var_y[] = $var_y[$c];

            $c++;

        }

        
				if ( isset( $_var_x[$c] ) ){
        	$mean_x = (float) $_var_x[$c];
				}

				if ( isset( $_var_y[$c] ) ){
        	$mean_y = (float) $_var_y[$c];
				}        

        $n = 0;

        

        while ($n < $c) {

            if (is_null($_var_x[$n]) || is_null($_var_y[$n]))    break;

            $n = (float) $n;

            $sweep = $n / ($n + 1.0);

            

            $float_x = (float) $_var_x[$n];

            $float_y = (float) $_var_y[$n];

            

            $delta_x = $float_x - $mean_x;

            $delta_y = $float_y - $mean_y;

            

            $sum_sq_x += $delta_x * $delta_x * $sweep;

            $sum_sq_y += $delta_y * $delta_y * $sweep;

            $sum_coproduct += $delta_x * $delta_y * $sweep;

            $mean_x += $delta_x / ($n + 1.0);

            $mean_y += $delta_y / ($n + 1.0);

            $n += 1;

        }

        $n = (float) $n;

        $pop_sd_x = sqrt($sum_sq_x / $n);

        $pop_sd_y = sqrt($sum_sq_y / $n);

        $cov_x_y = $sum_coproduct / $n;

        $result = $cov_x_y / ($pop_sd_x * $pop_sd_y);

        

        return $result;

    }



for ($i=0; $i<=100; $i++) {

    $ary1[] = rand(0, 100);

    $ary2[] = rand(0, 100);

}

print corrcoef($ary1, $ary2);



結果

0.12244969945723



タグ


sqrt 





PHPでExcelの統計関数を実装



週間人気ページランキング / 7-4 → 7-10
順位 ページタイトル抜粋 アクセス数
1 指定文字列より前を取得 | 文字列(テキスト処理) 9
2 URLから、トップページ(index.html)判定 | ルーティング 8
3 文字列に含まれるすべての半角空白と全角空白を削除 | テキスト処理 3
3 PHPで最後の「指定区切り文字」より後ろを取得 | 配列(型) 3
3 配列の空要素を削除した後、添字を振り直す … array_diff() | 配列(型) 3
3 後ろから3文字削除 / 「18:00:00」→「18:00」 | 文字列(テキスト処理) 3
3 配列キー内に「指定文字列が含まれるキー」と「その値」を抽出して、新しい配列を返す | 配列(型) 3
3 PHPコード 3
4 8桁の生年月日数字を分解 | 文字列(テキスト処理) 2
4 range( ) を使用して日付期間を配列生成 | 日付および時刻関連 2
5 平均値、分散、標準偏差を求める | 統計 1
5 ログイン 1
5 再帰処理で平坦な連想配列をツリー化 | 多次元配列(型) 1
5 array_map()第3引数 / 「2つの配列の要素数同士を順番に掛け合わせた配列」を返す | 配列(型) 1
5 対象文字列を、文字数に応じて半分に分割後、それぞれを取得 | 文字列(テキスト処理) 1
5 「半角スペースもしくは全角スペースで区切られた文字列」先頭に指定文字列を付与 | 文字列(テキスト処理) 1
5 配列の空要素を削除 / array_diff()で空配列と比較し配列の差分を返す | 配列(型) 1
5 マークダウンで「> 引用」を、HTMLエスケープ後、
へHTML変換 | エスケープ処理
1
5 foreachの中でswitch文 | 配列(型) 1
5 後ろから3文字取得 | 文字列(テキスト処理) 1
2025/7/11 1:02 更新