作者weiching0323 (许惟妮)
看板FJU-STAT95B
标题[专业] 5/4无母数课堂笔记15:40-16:30
时间Tue May 4 15:53:52 2010
================================从以下开始复制=================================
#include <C:\devcpp\all.h>
int main()
{
long int i, j;
srand((unsigned)time(NULL));
cout.setf(ios::fixed); cout.precision(10);
// $$$$$$$$ Start your program here $$$$$$$$$$$$$$$$$$$$$
double x[12] = {0,1,1,0,-2,6,7,7,4,-3,9,14};
Univar X(12,x);
double ind[12] = {1,1,1,1,1,0,0,0,0,0,0,0};
Univar IND(12,ind);
IND.print(0);
double C = 0.0;
Univar T1;
do
{
T1.push_back( X * IND );
C++;
}while(prev_permutation(IND.begin(),IND.end()));
cout << "number of permutation is " << C << endl;
cout << exp(ln_comb(12,5)) << endl;
T1.freq_dist();
cout << "P_value = " << endl;
cout.precision(4);
cout << 2.0*T1.emp_CDF(0) << endl;
Univar T2;
for(i=0; i < 10000; i++)
{
IND.shuffle();
T2.push_back( X * IND );
}
T2.freq_dist();
Univar t = time_axis(-4,43,1);
multi_dis G(t);
G.push_freq_dist(T1.freq_dist());
G.push_freq_dist(T2.freq_dist());
// G.print_to_file("c:\\devcpp\\G.dat",5);
cout << "Monte Carlo P_value = " << endl;
cout.precision(5);
cout << 2.0*T2.emp_CDF(0) << endl;
// $$$$$$$ End of your program $$$$$$$$$$$$$$$$$$$$$$$$$
==================================复制结束====================================
上色是新的
请继续follow^_<
又要开EXCLE 画图了
**今天画图不要用平滑线**
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 140.136.161.83
1F：推 ww770829:揪甘心 05/04 16:02
※ 编辑: weiching0323 来自: 140.136.161.83 (05/04 16:04)
2F：→ weiching0323:今天的问题好难喔ˊˋ 05/04 16:10
3F：嘘 seankevin:一楼交大哥快回答 05/04 16:17
4F：→ seankevin:等等补推 05/04 16:17
5F：推 seankevin:我补 05/04 16:20
6F：推 ww770829:让政大跟成大的来吧ˊˋ 05/04 16:22