x=read.table(“C:\Users\dell\Desktop\one.csv”,sep=’,’,header=F)#读取#
xbar=colMeans(x)#样本均值向量#
y=cov(x)#样本协方差#
m=(n-1)y#li离差#
n=cor(x)#相关矩阵#
u=c(7,5,4,8)
T2=19t(xbar-u)%%solve(n)%%(xbar-u)#统计量#
16/(19*4)*T2#F#
qf(1-0.05,4,16)#F分布临界值0.05#
x=read.table(“C:\Users\dell\Desktop\one.csv”,sep=’,’,header=F)#读取#
xbar=colMeans(x)#样本均值向量#
y=cov(x)#样本协方差#
m=(n-1)y#li离差#
n=cor(x)#相关矩阵#
u=c(7,5,4,8)
T2=19t(xbar-u)%%solve(n)%%(xbar-u)#统计量#
16/(19*4)*T2#F#
qf(1-0.05,4,16)#F分布临界值0.05#