#include "stdio.h"
#include "math.h"
/*形象与函数类型 参数意义
int n 自变量x的个数
int k 观测数据的点数
double x[k][n+1] 其中前n列存放自变量因子xi(i=0，1，…，n-1)的k次观测值；最后一列存放因变量y的k次观测值
double f1 欲选入因子时显著性检验的F-分布值
double f2 欲剔除因子时显著性检验的F-分布值
double eps 防止系数相关矩阵退化的判据
double xx[n+1] 前n个分量返回n个自变量因子的算术平均值 ；最后一个分量返回因变量y的算术平均值
double b[n+1] 返回回归方程中各因子的回归系数及常数项b0，b1，…，bn
double v[n+1] 前n个分量返回各因子的偏回归平方和Vi（i=0，1，…，n-1）；最后一个分量返回残差平方和q
double s[n+1] 前n个分量返回各因子回归系数的标准偏差si(i=0，1，…，n-1)；最后一个分量返回估计的标准偏差s
double dt[2] dt(0)返回复相关系数，dt(1)返回F-
double ye[k] 返回对应于k个观测值的因变量条件期望值的k个估计值ei(i=0,1,…，k-1)
double yr[k] 返回因变量的k个观测值的残差δi(i=0，1，…，k-1)
double r[n+1][n+1] 返回最终的规格化的系数相关矩阵R
void sqt3( ) 过程
*/
void sqt3(int n,int k,double x[],double f1,double f2,double eps,double xx[],double b[],double v[],double s[],double dt[],double ye[],double yr[],double r[])
{
int i,j,ii,m,imi,imx,l,it ;
double z,phi,sd,vmi,vmx,q,fmi,fmx ;
m=n+1 ;
q=0.0 ;
for(j=0;j<=n;j++)
{
z=0.0 ;
for(i=0;i<=k-1;i++)
z=z+x[i*m+j]/k ;
xx[j]=z ;
}
for(i=0;i<=n;i++)
for(j=0;j<=i;j++)
{
z=0.0 ;
for(ii=0;ii<=k-1;ii++)
z=z+(x[ii*m+i]-xx[i])*(x[ii*m+j]-xx[j]);
r[i*m+j]=z ;
}
for(i=0;i<=n;i++)
ye[i]=sqrt(r[i*m+i]);
for(i=0;i<=n;i++)
for(j=0;j<=i;j++)
{
r[i*m+j]=r[i*m+j]/(ye[i]*ye[j]);
r[j*m+i]=r[i*m+j];
}
phi=k-1.0 ;
sd=ye[n]/sqrt(k-1.0);
it=1 ;
while(it==1)
{
it=0 ;
vmi=1.0e+35 ;
vmx=0.0 ;
imi=-1 ;
imx=-1 ;
for(i=0;i<=n;i++)
