Matlab进行假设检验程序

Matlab进行假设检验程序：（以下均是m文件的程序）

统计原理可以参考浙江大学第四版教程的假设检验这一章1 %单个总体均值的检验%

function p1=T1_test(x,mu,alpha);

x=input('输入x的值：');

mu=input('输入mu的值：');%mu须检验的值%

alpha=input('输入alpha的值：');%alpha为显著水平%

n1=length(x);

x1=mean(x);

s1=std(x);

t1=abs((x1-mu)/(s1/n1^0.5));

p1=2*(1-tcdf(t1,n1-1));

disp('单个总体均值的检验的p值为：'),disp(p1);

2 %独立样本t检验%

%这里均是用p值法进行双边检验%

%条件：两总体方差相等，且总体样本x1,x2独立% function p=T_test2(x1,x2,alpha);

x1=input('输入x1的值：');

x2=input('输入x2的值：');

alpha=input('输入alpha的值：');

n1=length(x1);

n2=length(x2);

s1=std(x1);s2=std(x2);

sw=((n1-1)*s1^2+(n2-1)*s2^2)/(n1+n2-2);

t=mean(x1-x2)/(sw^0.5*(1/n1+1/n2)^0.5);

p=2*(1-tcdf(t,n1+n2-2));

disp('配对样本的t检验的p值为：'),disp(p);

3 %配对样本T检验%

function p=Paired_Samples_Test(x1,x2,alpha);

x1=input('输入x1的值：');

x2=input('输入x2的值：');

alpha=input('输入alpha的值：');

n=length(x1);

d=x1-x2;

d1=mean(d);

t=abs(d1/(std(d)/n^0.5));

p=2*(1-tcdf(t,n-1));

disp('配对样本的t检验的p值为：'),disp(p);