Displaying 1 result from an estimated 1 matches for "t_statisit".
Did you mean:
t_statisitic
2009 Nov 05
4
The equivalence of t.test and the hypothesis testing of one way ANOVA
...,N1),rep('y',N2))
> Y = c(x,y)
> fr = data.frame(Y=Y,A=as.factor(A))
> afit=aov(Y ~ A,fr)
>
> X=model.matrix(afit)
> B=afit$coefficients
> V=solve(t(X) %*% X)
>
> mse=tail(summary(afit)[[1]]$'Mean Sq',1)
> df=tail(summary(afit)[[1]]$'Df',1)
> t_statisitic=(B/(mse * sqrt(diag(V))))[[2]]
> 2*(1-pt(abs(t_statisitic),df))#the p-value from aov
[1] 0.1090802
>