作者ko7751 (雨茜)
看板FJU-STAT95B
標題[專業] 無母5/11 黃小娘生日耶 13:40~16:30 part2
時間Tue May 11 16:23:58 2010
#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[8]={12,14,15,10,12,15,17,16};
Univar DATA(8,data);
DATA.print(0);
Univar X,Y;
DATA.split(X,5,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,6) << endl;
double ind[8]={1,1,1,1,1,0,0,0};
Univar I(8,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();
Univar x=time_axis(63,77,1);
multi_dis G(x);
G.push_freq_dist(T.freq_dist());
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\\G2.dat",5);
// $$$$$$$ End of your program $$$$$$$$$$$$$$$$$$$$$$$$$
cin.get();
return 0;
}
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.136.161.80
1F：推 weiching0323:辛苦你了 我今天翹班一天 05/11 18:45
2F：→ ko7751:那怎麼沒扣你薪水 發給我啊!!! 05/11 23:14
3F：推 abidog: 這一篇文章值 1 元 05/12 17:42