[正經] 用腦時間

[ko7751 (雨茜)]

  double data[20]={22,34,52,62,30,40,64,84,56,59,
                   52,71,76,54,67,83,66,90,77,84};
  Univar DATA(20,data);
  DATA.print(0);
  Univar X,Y;
  DATA.split(X,10,Y);
  X.print(0);  Y.print(0);
  cout.precision(6);
  double t =t_stat(X,Y);
  cout<< "two sample t stat.=" << t << endl;
  cout<< t_2_tailed_prob(t,18) << endl;
  double ind[20]={1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0};
  Univar I(20,ind);
  I.print(0);
  double T1=DATA*I;
  cout<< " nonparametric randomization test stat. = " << T1 << endl;
  Univar T;
  for(i=0;i<1000000;i++)
  {
    I.shuffle();
    T.push_back(DATA*I);
  }
  T.freq_dist(-1);
  Univar x=time_axis(460,760,10);
  multi_dis G(x);
  G.push_freq_dist(T.freq_dist(-1));
  G.print_to_file("C:\\devcpp\\G.dat",5);

--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.136.161.80
推 ww770829:感謝~
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 14:56)
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 14:59)
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 15:01)
推 weiching0323:今天輪到你接手!!!!
推 weiching0323:幫你改標題!!!
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 15:10)
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 15:12)
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 15:23)
推 abidog:樓上人咧  竟然離開工作崗位!!