User Tools

Site Tools


c:ms:2025:schedule:week06_t-test_and_anova_note

Differences

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

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
c:ms:2025:schedule:week06_t-test_and_anova_note [2025/04/10 13:41] – [t-test, different approach] hkimscilc:ms:2025:schedule:week06_t-test_and_anova_note [2025/04/14 08:35] (current) – [t-test, different approach] hkimscil
Line 462: Line 462:
 abline(v = mean.a, lty=2, lwd=3, col="blue" abline(v = mean.a, lty=2, lwd=3, col="blue"
 abline(v = mean.b, lty=2, lwd=3, col="darkgreen" abline(v = mean.b, lty=2, lwd=3, col="darkgreen"
 +
 +# or
 +hist(A, xlim=c(0,50), ylim=c(0,8), col="red")
 +hist(B, add=T, col=rgb(0, 1, 0, 0.5))
 +abline(v=mean(dat), col="black")
 +abline(v=mean(A), col="red")
 +abline(v=mean(B), col="darkgreen")
 </code> </code>
  
Line 538: Line 545:
 f.calculated.pvalue f.calculated.pvalue
 </code> </code>
 +
 <code> <code>
-df.between  
-df.within  
 x <- rf(500000, df1 = df.between, df2 = df.within) x <- rf(500000, df1 = df.between, df2 = df.within)
 hist(x, hist(x,
      breaks = "Scott",      breaks = "Scott",
      freq = FALSE,      freq = FALSE,
-     xlim = c(0, 3), +     xlim = c(0, 7), 
-     ylim = c(0, 1),+     ylim = c(0, 3),
      xlab = "",      xlab = "",
      main = "Histogram for a F-distribution       main = "Histogram for a F-distribution 
Line 552: Line 558:
      cex.main = 0.9      cex.main = 0.9
 ) )
-curve(df(x, df1 = df.one, df2 = df.two),  +curve(df(x, df1 = df.between, df2 = df.within),  
-      from = 0, to = 4, n = 5000, +      from = 0, to = 8, n = 5000, 
       col = "red", lwd = 2,        col = "red", lwd = 2, 
       add = T)       add = T)
 abline(v=f.calculated, col="blue", lwd=2, lty="dotted") abline(v=f.calculated, col="blue", lwd=2, lty="dotted")
  
 +f.calculated.pvalue
 +1 - f.calculated.pvalue
 </code> </code>
 +
 <code> <code>
 # 한편,  t test를 했었을 때 (A, B 그룹을 가지고 independent  # 한편,  t test를 했었을 때 (A, B 그룹을 가지고 independent 
Line 671: Line 680:
 > abline(v = mean.b, lty=2, lwd=3, col="darkgreen" > abline(v = mean.b, lty=2, lwd=3, col="darkgreen"
  
 +> # or 
 +> hist(A, xlim=c(0,50), ylim=c(0,8), col="red")
 +> hist(B, add=T, col=rgb(0, 1, 0, 0.5))
 +> abline(v=mean(dat), col="black")
 +> abline(v=mean(A), col="red")
 +> abline(v=mean(B), col="darkgreen")
 </code> </code>
 +
 {{c:ms:2025:schedule:pasted:20250409-075631.png?550}} {{c:ms:2025:schedule:pasted:20250409-075631.png?550}}
 +{{:c:ms:2025:schedule:pasted:20250414-083514.png?550}}
 <code> <code>
 > # mean.total 에서 그룹a의 평균까지의 차이를 구한 후 > # mean.total 에서 그룹a의 평균까지의 차이를 구한 후
c/ms/2025/schedule/week06_t-test_and_anova_note.1744260114.txt.gz · Last modified: 2025/04/10 13:41 by hkimscil

Donate Powered by PHP Valid HTML5 Valid CSS Driven by DokuWiki