Advertisement
Guest User

Untitled

a guest
May 24th, 2015
206
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.03 KB | None | 0 0
  1. Ticker = ["ALTR", "BCP", "BPI", "EDP", "GALP", "IPR", "JMT", "EGL", "NOS", "PTI", "PTC", "RENE", "SEM", "SON"];
  2.  
  3. d1 = date(2010, 1, 1);
  4.  
  5. d2 = date(2014, 12, 31);
  6.  
  7. K = 0(#Ticker,3);
  8.  
  9. PvKS = 0(#Ticker,3);
  10.  
  11. for (t in 1:#Ticker) {
  12.  
  13. V0 = loaddaily(d1,d2,Ticker(t)+".csv","Volume");
  14.  
  15. I0 = loaddaily(d1,d2,Ticker(t)+".csv","Close");
  16.  
  17. V = select (V0,(I0>0));
  18.  
  19. I = select (I0,(I0>0));
  20.  
  21. L = [0;100*dif(ln(I))];
  22.  
  23. E = elaw(V);
  24.  
  25. Q1 = E.invpl(1/3);
  26.  
  27. Q2 = E.invpl(2/3);
  28.  
  29. L1 = select (L,(V<=Q1));
  30.  
  31. E1 = elaw(L1);
  32.  
  33. L2 = select (L,(V>Q1 & V<=Q2));
  34.  
  35. E2 = elaw(L2);
  36.  
  37. L3 = select (L,(V>Q2));
  38.  
  39. E3 = elaw(L3);
  40.  
  41. x = [min(L):max(L):(max(L)-min(L))/100];
  42.  
  43. K(t,1) = ((#L1*#L2/(#L1+#L2))^0.5)*(max(abs(E1.pl(x)-E2.pl(x)))); // К - статистика критерия Колмогорова-Смирнова
  44.  
  45. K(t,2) = ((#L1*#L3/(#L1+#L3))^0.5)*(max(abs(E1.pl(x)-E3.pl(x))));
  46.  
  47. K(t,3) = ((#L3*#L2/(#L3+#L2))^0.5)*(max(abs(E3.pl(x)-E2.pl(x))));
  48.  
  49. PvKS(t,1) = pvKolm(K(t,1)); // Р-значения критерия Колмогорова-Смирнова
  50.  
  51. PvKS(t,2) = pvKolm(K(t,2));
  52.  
  53. PvKS(t,3) = pvKolm(K(t,3));}
  54.  
  55. savetable(["Тикер", "Малый и средний", "Малый и большой", "Средний и большой"; 'Ticker,K],"Табл.9.Значения статистики критерия Колмогорова-Смирнова.csv");
  56.  
  57. savetable(["Тикер", "Малый и средний", "Малый и большой", "Средний и большой"; 'Ticker,round(PvKS,4)],"Табл.10.P-значения критерия Колмогорова-Смирнова.csv");
  58.  
  59. //Задаем интервал от нуля до единицы сшагом 0,01
  60.  
  61. P=0:1:0.01;
  62.  
  63. //Рассчитываем значение критерия Колмогорова для проверки гипотезы о равномерности распределения P-значений критерия Колмогорова-Смирнова
  64.  
  65. d= maxabs(elaw(PvKS).pl(P) - ulaw(0,1).pl(P));
  66.  
  67. u=d*(((#PvKS)^0.5));
  68.  
  69. //Задаем условие, что если критерий Колмогорова < критического значения статистики для уровня значимости 0,05, то Р-значения распределены равномерно
  70.  
  71. //А если условие не выполняется, то Р-значения не распределены равномерно
  72.  
  73. //if (u<crKolm(0.05)) itog = "P-значения распределены равномерно";
  74.  
  75. //else itog = "Р-значения распределены не равномерно";
  76.  
  77. if (u<crKolm(0.01)) itog = "P-значения распределены равномерно";
  78.  
  79. else itog = "Р-значения распределены не равномерно";
  80.  
  81. //Построение гистограммы частот P-значений критерия Колмогорова-Смирнова
  82.  
  83. wintitle(itog); hist(PvKS.intfrgel(0:1:0.1),blue);axes();
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement