Differences

This shows you the differences between two versions of the page.

--- t-test_summing_up [2025/09/18 08:44] – [8] hkimscil
+++ t-test_summing_up [2025/09/18 08:45] (current) – [t-test summing up] hkimscil
@@ Line 1: / Line 1: @@
 ====== t-test summing up ======
 <code>
 rm(list=ls())
@@ Line 210: / Line 209: @@
 dat <- stack(comb)
 head(dat)
+tail(dat)
 m.tot <- mean(s.all)
@@ Line 216: / Line 216: @@
 ss.tot <- ss(s.all)
+bet.s1 <- (m.tot - m.s1)^2 * length(s1)
+bet.s2 <- (m.tot - m.s2)^2 * length(s2)
+ss.bet <- bet.s1 + bet.s2
 ss.s1 <- ss(s1)
 ss.s2 <- ss(s2)
+ss.wit <- ss.s1+ss.s2
-df.tot <- length(s.all)-1
+ss.tot
-df.s1 <- length(s1)-1
-df.s2 <- length(s2)-1
-ms.tot <- var(s.all)
-ms.tot
-ss.tot/df.tot
-var(s1)
-ss.s1 / df.s1
-var(s2)
-ss.s2 / df.s2
-ss.b.s1 <- length(s1) * ((m.tot - m.s1)^2)
-ss.b.s2 <- length(s2) * ((m.tot - m.s1)^2)
-ss.bet <- ss.b.s1+ss.b.s2
 ss.bet
-ss.wit <- ss.s1 + ss.s2
 ss.wit
+ss.bet+ss.wit
-ss.bet + ss.wit
+df.tot <- length(s.all) - 1
-ss.tot
+df.bet <- nlevels(dat$ind) - 1
+df.s1 <- length(s1)-1
+df.s2 <- length(s2)-1
+df.wit <- df.s1 + df.s2
-library(dplyr)
+df.tot
-# df.bet <- length(unique(dat)) - 1
-df.bet <- nlevels(dat$ind) - 1
-df.wit <- df.s1+df.s2
 df.bet
 df.wit
 df.bet+df.wit
-df.tot
+ss.tot/df.tot
+ms.tot <- ss.tot/df.tot
 ms.bet <- ss.bet / df.bet
 ms.wit <- ss.wit / df.wit
-ms.bet
-ms.wit
 f.cal <- ms.bet / ms.wit
 f.cal
 pf(f.cal, df.bet, df.wit, lower.tail = F)
 f.test <-  aov(dat$values~ dat$ind, data = dat)
@@ Line 273: / Line 259: @@
 t.cal.ts
-# this is anova after all.
+# the above is anova after all.
 m1 <- lm(dat$values~dat$ind, data = dat)
@@ Line 287: / Line 273: @@
 sum.m1$fstatistic[1]
 ms.bet/ms.wit
 </code>
 ====== t-test summing up output ======