作者weiching0323 (許惟妮)
看板FJU-STAT95B
標題[專業] 6/8無母數課堂筆記15:40-16:30
時間Tue Jun 8 16:01:23 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 $$$$$$$$$$$$$$$$$$$$$
Matrix Y;
Y.read_from_file("c:\\devcpp\\anacova.dat",2);
Y.print(0);
Matrix X2;
Y.split_column(Y,1,X2);
Y.print(0); X2.print(0);
Matrix D = gen_dummy(One_way_class(3,5));
D.print(0);
Matrix SST, SSR, SSE;
Matrix X = J_Mat(Y.row(),1);
X.append_column(X2);
X.print(0);
Part_SS_X_full(Y,X,SST,SSR,SSE);
SSR.print(2);
SSE.print(2);
SST.print(2);
X.append_column(D);
X.print(0);
Part_SS_X_full(Y,X,SST,SSR,SSE);
cout << "full Model SS = " << endl;
SSR.print(2);
SSE.print(2);
SST.print(2);
double F = (13.28/2)/(27.99/11);
cout << "F = " << F << endl;
double P_value = F_r_tailed_prob(F,2,11);
cout.precision(4);
cout << " P_value = " << P_value << endl;
Partial_F_test(Y,X,2,3);
Partial_F_test(Y,X,1,1);
Matrix B = (X.T()*X).inv()*X.T()*Y;
B.print(4);
Partial_t_test(Y,X);
Univar SSER;
for(i=0; i < 10000; i++)
{
X.randomize_row(1,1);
Part_SS_X_full(Y, X, SST, SSR, SSE);
SSER.push_back(SSE);
}
SSER.freq_dist();
cout.precision(4);
cout << "nonparametric randomization P_value = ";
cout << SSER.emp_CDF(27.99) << endl;
Univar x = time_axis(45,206,1);
multi_dis G(x);
G.push_freq_dist(SSER.freq_dist());
G.print_to_file("c:\\devcpp\\G2.dat",4);
// $$$$$$$ End of your program $$$$$$$$$$$$$$$$$$$$$$$$$
cin.get();
return 0;
}
=====================================複製結束==================================
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.136.161.81
※ 編輯: weiching0323 來自: 140.136.161.81 (06/08 16:15)