『Expert advisor』は、おもしろい！: 2012/02/19

2012/02/19

GSLでヒストグラム■■■

前回の内容に、ヒストグラムを追加しました。ファイルはココにまとめてあります。

DLL追加コード

#include ".\include\gsl\gsl_histogram.h"//ヒストグラム
//+-------------------------------------------------------------------+
//|ヒストグラム (Histgram)
//|input[inp_N] 入力用バッファ
//|out[N] 出力用バッファ
//|inp_N 入力データ数
//|N 階級数
//|min 最小値(含む）
//|max 最大値(含まない）
//+-------------------------------------------------------------------+
__declspec(dllexport) void __stdcall Histgram(double *input,double *out,int inp_N,int N,double min,double max){
int i;
//メモリの確保
gsl_histogram * h = gsl_histogram_alloc(N);
if(h==NULL)exit(0);
//階級のセット
gsl_histogram_set_ranges_uniform(h,min,max);
//実行
for(i=0;i<inp_N;i++)gsl_histogram_increment(h,input[i]);
//出力
for(i=0;i<N;i++)out[i]=gsl_histogram_get(h,i);
//メモリの解放
gsl_histogram_free(h);
}

#include ".\include\gsl\gsl_histogram.h"//ヒストグラム

//+-------------------------------------------------------------------+
//|ヒストグラム (Histgram)
//|input[inp_N]  入力用バッファ
//|out[N]      出力用バッファ
//|inp_N         入力データ数
//|N             階級数
//|min           最小値(含む）
//|max           最大値(含まない）
//+-------------------------------------------------------------------+
__declspec(dllexport) void __stdcall Histgram(double *input,double *out,int inp_N,int N,double min,double max){
	int i;
	//メモリの確保
	gsl_histogram * h = gsl_histogram_alloc(N);
	if(h==NULL)exit(0);
	//階級のセット
	gsl_histogram_set_ranges_uniform(h,min,max);
	//実行
	for(i=0;i<inp_N;i++)gsl_histogram_increment(h,input[i]);
	
	//出力
	for(i=0;i<N;i++)out[i]=gsl_histogram_get(h,i);
	

	//メモリの解放
	gsl_histogram_free(h);
}

MQLコード

//+------------------------------------------------------------------+
//|                                               ヒストグラム.mq4 　|
//|                                       Copyright ｩ 2012, bighope  |
//|                      http://expertadviser-bighope.blogspot.com/  |
//+------------------------------------------------------------------+
#property copyright "Copyright ｩ 2012, bighope"
#property link      "http://expertadviser-bighope.blogspot.com/"
#property indicator_separate_window
#property indicator_minimum 0
#include <GSL.mqh>
#property indicator_buffers 3
#property indicator_color1 Snow
#property indicator_color2 Yellow
#property indicator_color3 Red
extern string  SymbA   =   "USDJPY";
extern int     TframeA =   PERIOD_H1;
extern string  SymbB   =   "EURUSD";
extern int     TframeB =   PERIOD_H1;
extern string  SymbC   =   "AUDUSD";
extern int     TframeC =   PERIOD_H1;
extern int     Range   =   2000;//データ数
extern int     N       =    100;//階級数
extern double  min     =    0.0;//階級最小値(pips)
extern double  max     =    2000.0;//階級最大値(pips)

double Ap[];
double Bp[];
double Cp[];
//+------------------------------------------------------------------+
//| Custom indicator initialization function                         |
//+------------------------------------------------------------------+
int init()
  {
//---- indicators
       SetIndexStyle(0,DRAW_LINE);
       SetIndexStyle(1,DRAW_LINE);
       SetIndexStyle(2,DRAW_LINE);

       SetIndexBuffer(0,Ap);
       SetIndexBuffer(1,Bp);
       SetIndexBuffer(2,Cp);

//----
   return(0);
  }
//+------------------------------------------------------------------+
//| Custom indicator deinitialization function                       |
//+------------------------------------------------------------------+
int deinit()
  {
//----

//----
   return(0);
  }
//+------------------------------------------------------------------+
//|GetVal
//|symb        通貨ペア名
//|timeframe   タイムフレーム
//|shift       Barの位置
//+------------------------------------------------------------------+
double GetVal(string symb,int timeframe,int shift){
   double Hi,Lw;

   Hi = iHigh(symb,timeframe,shift);
   Lw = iLow(symb,timeframe,shift);

return(Hi-Lw);
}
//+------------------------------------------------------------------+
//|GetData
//|symb        通貨ペア名
//|timeframe   タイムフレーム
//|out[Ns]     データの取得配列
//|Ns          データの取得数
//+------------------------------------------------------------------+
void GetData(string symb,int timeframe,double &out[],int Ns){
   for(int i=0;i<Ns;i++)out[i]=GetVal(symb,timeframe,i);
}
//+------------------------------------------------------------------+
//|GetHist
//|symb        通貨ペア名
//|timeframe   タイムフレーム
//|out[Ns]     データの取得配列
//|Ns          データの取得数
//+------------------------------------------------------------------+
void GetHist(string symb,int timeframe,int Ns,double &data[]){
   int i;
   double count,minpt,maxpt,pt;
   double input[];
   double output[];
   ArrayResize(input,Ns);
   ArrayResize(output,N);
   ArrayInitialize(data,EMPTY_VALUE);

   GetData(symb,timeframe,input,Ns);
   pt = MarketInfo(symb,MODE_POINT);
   minpt = min*pt;
   maxpt = max*pt;
   Histgram(input,output,Ns,N,minpt,maxpt);
   for(i=0;i<N;i++)count += output[i];
   if(count==0)count=1.0;
   for(i=0;i<N;i++)data[i]=output[i]/count;
}

//+------------------------------------------------------------------+
//| Custom indicator iteration function                              |
//+------------------------------------------------------------------+
int start()
  {
   int    counted_bars=IndicatorCounted();
//----
   if(counted_bars<1){
      GetHist(SymbA,TframeA,Range,Ap);
      GetHist(SymbB,TframeB,Range,Bp);
      GetHist(SymbC,TframeC,Range,Cp);
   }

//----
   return(0);
  }
//+------------------------------------------------------------------+

作動・確認

デフォルトで作動させると以下の様になります。
※縦軸が、発生確率で横軸が階級です。（チャート右が最少数）(白：USDJPY 1H 黄：EURUSD 1H 赤：AUDUSD 1H)

※お化粧が全くなされていないので見苦しいかもしれませんがご了承ください。＾＾；

USDJPYが特異な特徴と示す？

いろいろ試してみるとUSDJPYだけがなんだか特異な形状を示します。なぜ？
例えば、(白：EURJPY 1H 黄：EURUSD 1H 赤：USDCHF 1H)だと

の様に似たような形になります。
他のペアも試してみましたが、試したペアの中でUSDJPYだけが極端に形状が異なっています。

また、(白：USDJPY 4H 黄：EURUSD 1H 赤：GBPJPY 1H)とした場合は、

となり、似たような形が出現します。

しかし、(白：USDJPY 1D 黄：EURUSD 1D 赤：AUDUSD 1D)となると、

特異な形状が解消されたように見えます。（ちなみに、max=5000で表示)

まとめ

全てのペア及びタイムフレームを確認したわけでありませんが、ドル円だけがなぜ特異な形となるのか疑問です。他ペアに比べて突発的事情（イベントやニュース）に敏感に反応（突発的なボラ）し、その他のボラは、低い事を示しているのでしょうか？疑問が残るところです。