[專業] 無母5/11 黃小娘生日耶 13:40~16:30

[ko7751 (雨茜)]

#include <C:\devcpp\all.h>
double mean,stdev;

double Normal_den(double x)
{  return nor_den(x, mean,stdev);
}

int main()
{
  long int i, j;
  srand((unsigned)time(NULL));
  cout.setf(ios::fixed);  cout.precision(10);


// $$$$$$$$  Start your program here $$$$$$$$$$$$$$$$$$$$$

  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);
  cout.precision(6);
  double P=T.emp_CDF(T1);
  cout<< " nonparametric randomization test stat. = " << 2.0*P << endl;
  double L=P-1.96*sqrt(P*(1.0-P)/1000000),
         R=P+1.96*sqrt(P*(1.0-P)/1000000);
  cout<< "C.I. for P_value = ("<<2.0*L<<","<<2.0*R<<")"<<endl;

  mean=T.mean();
  stdev=T.stdev();
  G.push_function(Normal_den);
  G.print_to_file("C:\\devcpp\\G.dat",5);

// $$$$$$$   End of your program $$$$$$$$$$$$$$$$$$$$$$$$$


  cin.get();
  return 0;
}
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.136.161.80
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 16:01)
※ 編輯: ko7751          來自: 140.136.161.80       (05/11 16:11)
推 abidog:生日快樂  我的男朋友
推 wen771020:生日快樂 小胖的女朋友
推 weiching0323: 快樂 小胖的女朋友
→ ko7751:   生日快樂 小胖的女朋友 淑齡的男朋友 = =
推 timmac: 生日快 樂林淑齡的男朋友
→ ko7751:   生日快樂 小胖的女朋友 淑齡的男朋友  樓上再吃醋!!!!
推 abidog:   生日快樂 小胖的女朋友 現在是我搶輸的意思嗎
推 weiching0323:樓上是妾
→ ko7751:黃小娘很搶手的意思??!