スポンサーサイト

上記の広告は１ヶ月以上更新のないブログに表示されています。新しい記事を書く事で広告が消せます。
Pagetop

PHPでExcelの統計関数を実装する

PHPで統計関数を使いたかったのだが、探しても探しても存在しないので、自前で実装しました。
わかりやすく、Excelの統計関数に[s_]を付ければ、引数も同様に使える様になっています。
[s_]を付けた理由は、var関数がPHPでは予約語になってしまっているからです。


↓↓同じ境遇の人の為に公開（ご自由にどうぞ）↓↓
↓↓（使ってもらえる際には、コメントもらえると嬉しいです）↓↓

【関数名一覧】＝＝＝＝＝
関数名説明
s_average単純平均を取得する
s_harmean調和平均を取得する
s_varp分散を取得する
(引数を母集団全体であると仮定して、母集団の分散を返す)
s_stdevp標準偏差を取得する
s_devp偏差値を取得する
s_correl相関係数を取得する
s_intercept回帰直線のｙ切片を求める
（[y = a * x + b]の[b]を求める）
s_slope回帰直線の傾きを求める
（[y = a * x + b]の[a]を求める）
s_forecast回帰直線上でｘに対応するｙを求める
s_trend（＊空実装）【現在未対応】
直線ｙ＝aｘ＋ｂのｙの値を求める
s_var（＊空実装）【現在未対応】
分散を取得する（不偏分散）
(引数を母集団の標本であると見なして、母集団に対する分散を返す)
s_stdev（＊空実装）【現在未対応】
標準偏差を取得する（不偏分散）


【ソース】＝＝＝＝＝＝＝

/**
 * 単純平均を取得する
 */
function s_average($list) {
    if(count($list) < 2 || !is_array($list)){
        return false;
    }
    $count = count($list);
    for ($i = 0; $i < $count; $i++) {
        $sum += $list[$i];
    }
    return $sum / $count;
}

/**
 * 調和平均を取得する
 */
function s_harmean($list) {
    if(count($list) < 2 || !is_array($list)){
        return false;
    }
    $count = count($list);
    for ($i = 0; $i < $count; $i++) {
        $sum += 1 / $list[$i];
    }
    return $count / $sum;
}

/**
 * 分散を取得する
 * (引数を母集団全体であると仮定して、母集団の分散を返す)
 */
function s_varp($list) {
    if(count($list) < 2 || !is_array($list)){
        return false;
    }
    $avg = s_average($list);
    $dec = 0;
    $count = count($list);
    for($i = 0; $i < $count; $i++){
        $dec += pow($avg - $list[$i], 2);
    }
    return $dec / $count;
}

/**
 * 標準偏差を取得する
 */
function s_stdevp($list) {
    if(count($list) < 2 || !is_array($list)){
        return false;
    }
    return sqrt(s_varp($list));
}

/**
 * 偏差値を取得する
 */
function s_devp($list, $val) {
    if(count($list) < 2 || !is_array($list)){
        return false;
    }
    return 10 * ($val - s_average($list)) / s_stdevp($list) + 50;
}

/**
 * 相関係数を取得する
 */
function s_correl($list_y, $list_x){
    if(count($list_x) < 2 || count($list_y) < 2 || count($list_x) != count($list_y)){
        return false;
    }
    $avg_x = s_average($list_x);
    $stdevp_x = s_stdevp($list_x);
    $avg_y = s_average($list_y);
    $stdevp_y = s_stdevp($list_y);
   
    $count = count($list_x);
    for($i = 0; $i < $count; $i++){
        $x = $list_x[$i] - $avg_x;
        $y = $list_y[$i] - $avg_y;
        $devsum[$i] = $x * $y;
    }
   
    $avg_devsum = s_average($devsum);
    return $avg_devsum / ($stdevp_x * $stdevp_y);
}

/**
 * 回帰直線のｙ切片を求める
 * （[y = a * x + b]の[b]を求める）
 */
function s_intercept($list_y, $list_x){
    if(count($list_x) < 2 || count($list_y) < 2 || count($list_x) != count($list_y)){
        return false;
    }
    $count = count($list_x);
    for ($i = 0; $i < $count; $i++) {
        $x = $list_x[$i];
        $y = $list_y[$i];
        $x_sum += $x;
        $y_sum += $y;
        $xx_sum += $x * $x;
        $xy_sum += $x * $y;
    }
    $a = ($count * $xy_sum - $x_sum * $y_sum)/($count * $xx_sum - $x_sum * $x_sum);
    $b = ($y_sum - $a * $x_sum)/$count;
   
    return $b;
}

/**
 * 回帰直線の傾きを求める
 * （[y = a * x + b]の[a]を求める）
 */
function s_slope($list_y, $list_x){
    if(count($list_x) < 2 || count($list_y) < 2 || count($list_x) != count($list_y)){
        return false;
    }
    $count = count($list_x);
    for ($i = 0; $i < $count; $i++) {
        $x = $list_x[$i];
        $y = $list_y[$i];
        $x_sum += $x;
        $y_sum += $y;
        $xx_sum += $x * $x;
        $xy_sum += $x * $y;
    }
    $a = ($count * $xy_sum - $x_sum * $y_sum)/($count * $xx_sum - $x_sum * $x_sum);
    $b = ($y_sum - $a * $x_sum)/$count;
   
    return $a;
}

/**
 * 回帰直線上でｘに対応するｙを求める
 */
function s_forecast($target_x, $list_y, $list_x){
    if(count($list_x) < 2 || count($list_y) < 2 || count($list_x) != count($list_y)){
        return false;
    }
    $a = s_slope($list_y, $list_x);
    $b = s_intercept($list_y, $list_x);
   
    return $a * $target_x + $b;
}

/**
 * 回帰直線上の上昇率を求める
 */
function s_rate($list_y, $list_x){
    if(count($list_x) < 2 || count($list_y) < 2 || count($list_x) != count($list_y)){
        return false;
    }
    $a = s_slope($list_y, $list_x);
    $b = s_intercept($list_y, $list_x);
   
    return ($a * $list_x[count($list_x) - 1] + $b) / ($a * $list_x[count($list_x) - 2] + $b);
}

/**
 * 【現在未対応】
 * 直線ｙ＝aｘ＋ｂのｙの値を求める
 */
function s_trend($list_y, $list_x, $target_x, $flg=true){
}

/**
 * 【現在未対応】
 * 分散を取得する（不偏分散）
 * (引数を母集団の標本であると見なして、母集団に対する分散を返す)
 */
function s_var(){
}

/**
 * 【現在未対応】
 * 標準偏差を取得する（不偏分散）
 */
function s_stdev($list) {
    //return sqrt(var($list));
}


参考：
http://noraphp.blog73.fc2.com/blog-entry-7.html
http://q.hatena.ne.jp/1084333525
http://kogolab.jp/elearn/icecream/chap2/sec2.html
http://www.geocities.co.jp/Technopolis/2890/tool/Excel/regression.html


テーマ : プログラミング
ジャンル : コンピュータ

tag : PHP プログラミング 統計 Excel

Pagetop

勝手にメモ

引用させて頂いた皆様へ
まずは、素晴らしい情報を有難う御座います。
自分用にメモ、トラックバック貼らせて頂きましたm(_ _)m
ご迷惑でしたら、削除いたしますので、ご連絡頂ければ幸いです。

ファイルまたはディレクトリのみへの chmod コマンド

特定のユーザー（アカウント）に「パスワード無し」でsudoコマンド操作を許したい

テーマ : CentOS
ジャンル : コンピュータ

Pagetop

VMWareでネットワークに接続出来ない場合の対応

ちょっと、VMwareで「？」な現象が発生したので、メモ

VMWareで突然ネットワークに接続出来なくなった。
但し、ネットワークの設定は変更していない。
元々接続はできており、あるミドルウェアを入れたら突然通信不能になった。

簡単な確認で、デフォルトゲートウェイにpingを打ってみたら、通らない。
試しに、外部からVMwareのゲストOSにpingを打ってみたが、やはり通らない。
試しに、service network restart をやってみたが、やはり通らない。

なんでだ？？
と思い、色々調べてみた所・・・・原因不明ｗｗ

とりあえず、一旦ゲストOSを止め、ネットワークの設定で、
ラジオボタンのところの、BrigedをCustomに変更したら解決した。
何でこの設定で直ったんだ？この変更は、特段設定の内容に変更はないはずだが？

とにかく、よく分からないがこれで解決した。
（考察するに、VMWare内のデフォルトネットワークの設定のどこかが壊れたのかな？）

【筆者が勧める関連書籍】
↓VMWareを使用している人は是非一読してみて下さい。私的にはかなりヒットでした。

テーマ : VMware
ジャンル : コンピュータ

Pagetop

【連載】DRBD+HeartBeatでクラスタを構築する��

さて、今回は、
■＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝■
１、CentOSのインストール＆設定
２、DRBDのインストール＆設定
３、PostgreSQLのインストール＆設定
４、Apacheのインストール＆設定
５、HeartBeatのインストール＆設定
■＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝■
の、
「２、DRBDのインストール＆設定」
をご紹介します。

#OSインストール時に作成した自動マウントを解除する
作業対象：ha_active、ha_stanby
vi /etc/fstab
　　/var/lib/pgsqlの行を削除
　　/var/lib/php/sessionの行を削除
　　/var/www/htmlの行を削除

#現在のマウントを解除する
作業対象：ha_active、ha_stanby
umount /var/lib/pgsql
umount /var/lib/php/session
umount /var/www/html

#名前解決用にhostsに登録
作業対象：ha_active
vi /etc/hosts
　　# Do not remove the following line, or various programs
　　# that require network functionality will fail.
　　127.0.0.1                ha_active localhost.localdomain localhost
　　192.168.0.12         ha_stanby
作業対象：ha_stanby
vi /etc/hosts
　　# Do not remove the following line, or various programs
　　# that require network functionality will fail.
　　127.0.0.1                ha_stanby localhost.localdomain localhost
　　192.168.0.11         ha_active

#############
# drbd.conf  - DRBDデバイスの設定ファイル
# drbd    - DRBDを起動/停止するためのスクリプト
# drbdadm   - DRBDの管理ツール
# drbddisk   - デバイスをプライマリに切り替えてファイルシステムをマウントする
# drbdmeta   - DRBDのメタデータ管理ツール
# drbdsetup  - DRBDのセットアップツール
#############

作業対象：ha_active、ha_stanby
yum -y install drbd83 kmod-drbd83
chkconfig drbd on
vi /etc/drbd.conf
#
# please have a a look at the example configuration file in
# /usr/share/doc/drbd83/drbd.conf
#
global {
  usage-count yes;
}
common {
  protocol C;
}
resource r0 {
  on ha_active {
    device /dev/drbd0;
    disk /dev/mapper/VolGroup01-LogVol00;
    address 172.16.0.11:7789;
    meta-disk internal;
  }
  on ha_stanby {
    device /dev/drbd0;
    disk /dev/mapper/VolGroup01-LogVol00;
    address 172.16.0.12:7789;
    meta-disk internal;
  }
}
resource r1 {
  on clsact.itoparts-net.st.wakwak.ne.jp {
    device /dev/drbd1;
    disk /dev/mapper/VolGroup01-LogVol01;
    address 172.16.0.11:7790;
    meta-disk internal;
  }
  on clsstn.itoparts-net.st.wakwak.ne.jp {
    device /dev/drbd1;
    disk /dev/mapper/VolGroup01-LogVol01;
    address 172.16.0.12:7790;
    meta-disk internal;
  }
}
resource r2 {
  on clsact.itoparts-net.st.wakwak.ne.jp {
    device /dev/drbd2;
    disk /dev/mapper/VolGroup01-LogVol02;
    address 172.16.0.11:7791;
    meta-disk internal;
  }
  on clsstn.itoparts-net.st.wakwak.ne.jp {
    device /dev/drbd2;
    disk /dev/mapper/VolGroup01-LogVol02;
    address 172.16.0.12:7791;
    meta-disk internal;
  }
}

#OSインスト時に領域を作成してしまっていた為、DRBD用のデバイスを破壊する
作業対象：ha_active、ha_stanby
dd if=/dev/zero bs=1M count=1 of=/dev/mapper/VolGroup01-LogVol00; sync
dd if=/dev/zero bs=1M count=1 of=/dev/mapper/VolGroup01-LogVol01; sync
dd if=/dev/zero bs=1M count=1 of=/dev/mapper/VolGroup01-LogVol02; sync

#DRBD管理用のメタデータを同一パーティション上に作成する（別パーティションにする設定もある）
作業対象：ha_active、ha_stanby
drbdadm create-md r0
　　yes
drbdadm create-md r1
drbdadm create-md r2

#DRBDの実行
作業対象：ha_active、ha_stanby
/etc/init.d/drbd start

# DRBDディスクをフォーマット＆同期する
作業対象：ha_active
drbdsetup /dev/drbd0 primary --overwrite-data-of-peer
drbdsetup /dev/drbd1 primary --overwrite-data-of-peer
drbdsetup /dev/drbd2 primary --overwrite-data-of-peer
mkfs.ext3 /dev/drbd0
mkfs.ext3 /dev/drbd1
mkfs.ext3 /dev/drbd2

##### 同期が完了するまでに3~10時間かかります #####
##### ↓↓↓↓↓待てない人向け↓↓↓↓↓↓ #####
　　#同期速度の一時的な速度変更
　　#設定(Gigabit Eathの場合)
　　drbdsetup /dev/drbd0 syncer -r 3M
　　drbdsetup /dev/drbd1 syncer -r 1M
　　drbdsetup /dev/drbd2 syncer -r 2M
　　#元に戻す
　　drbdadm adjust r0
　　drbdadm adjust r1
　　drbdadm adjust r2
##### ↑↑↑↑↑↑↑待てない人向け↑↑↑↑↑↑↑ #####

#同期状態の確認
cat /proc/drbd
##### ↓↓↓↓↓同期確認がメンドクサイ人向け↓↓↓↓↓↓ #####
#!/bin/bash
while :
do
        clear;
        cat /proc/drbd;
        sleep 0.5;
done
##### ↑↑↑↑↑↑↑同期確認がメンドクサイ人向け↑↑↑↑↑↑↑ #####


#動作確認
作業対象：ha_active
mount /dev/drbd0 /var/lib/pgsql/
mount /dev/drbd1 /var/lib/php/session/
mount /dev/drbd2 /var/www/html/
df

#サンプルファイルの作成
作業対象：ha_active
echo "test" > /var/lib/pgsql/test
echo "test" > /var/lib/php/session/test
echo "test" > /var/www/html/test

#切り替え
作業対象：ha_active
umount /var/lib/pgsql/
umount /var/lib/php/session/
umount /var/www/html/
df
drbdadm secondary all

# Secondaryになっているか確認
作業対象：ha_active
cat /proc/drbd

#セカンダリで実行
作業対象：ha_stanby
drbdadm primary all

#プライマリになっているか確認
作業対象：ha_stanby
cat /proc/drbd

#切り替え後のプライマリでマウント
作業対象：ha_stanby
mount /dev/drbd0 /var/lib/pgsql/
mount /dev/drbd1 /var/lib/php/session/
mount /dev/drbd2 /var/www/html/
df

# 作成したファイルがあるかどうか確認
作業対象：ha_stanby
ls /var/lib/pgsql/test
ls /var/lib/php/session/test
ls /var/www/html/test

#もとに戻す

 




【筆者が勧める関連書籍】



テーマ : CentOS
ジャンル : コンピュータ

tag : CentOS クラスタリング HEARTBEAT DRBD

Pagetop

【連載】DRBD+HeartBeatでクラスタを構築する��

さて、今回は、
■＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝■
１、CentOSのインストール＆設定
２、DRBDのインストール＆設定
３、PostgreSQLのインストール＆設定
４、Apacheのインストール＆設定
５、HeartBeatのインストール＆設定
■＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝■
の、
「１、CentOSのインストール＆設定」
をご紹介します。

＊ディスクの設定は各々の環境に合わせて読み替えていただきたい。

インストール方式
　　GUIインストール
言語
　　Japanese(日本語)
キーボード
　　日本語
パーティション
　　カスタムレイアウト
　　（今回は、DRBDで共有するパーティションを先に切っておきます。）
　　（理由としては、後からパーティションを追加とかするのがメンドクサイからですｗ）
構成FS容量マウントポイント備考
/dev/sda    
　　/dev/sda1ext3101MB/bootブート領域
　　/dev/sda2――  
　　/dev/sda3――  
LVMFS容量マウントポイント備考
VolGroup00   VG
　　LogVol00―4096MBswapスワップ領域（メモリ量の2倍）
　　LogVol01ext361440MB/その他全て
VolGroup01   VG
　　LogVol00ext361440MB/var/lib/pgsqlPostgreSQLのデータ領域
　　LogVol01ext36144MB/bar/lib/php/sessionPHPのセッションファイル格納先
　　LogVol02ext315360MB/var/www/htmlApacheのコンテンツ領域

ブートローダ
　　/dev/sda

ホスト・ネットワーク

ホスト名備考Ethernet自動起動IPアドレス備考
ha_activeクラスタの稼動側ifcfg-eth0on192.168.0.11

ルータと接続
通常の通信用
HeartBeat（ICMP）通信用

  ifcfg-eth1on172.16.0.11

ha_stanbyと直接接続
DRBD通信用
HeartBeat（UDP）通信用

ha_stanbyクラスタの待機側ifcfg-eth0on192.168.0.12ルータと接続
  ifcfg-eth1on172.16.0.12ha_activeと直接接続
DRBD通信用
HEARTBEAT（UDP）通信用


ゲートウェイ・DNS

項目名IPアドレス
デフォルトゲートウェイ192.168.0.1
DNS１111.111.111.111
DNS２222.222.222.222

時間
　　アジア/東京
　　UTC使用する

パスワード
　　tekitou

パッケージ

言語
　　日本語のサポート
ベースシステム
　　ベース
アプリケーション
　　エディタ
　　テキストベースのインターネット
開発
　　開発ツール

＊上記以外はインストールしない

 

続きは次回



【筆者が勧める関連書籍】



テーマ : CentOS
ジャンル : コンピュータ

tag : CentOS クラスタリング HEARTBEAT DRBD

Pagetop
ホーム  次のページ »
Pagetop