Permalink
Browse files

3번 문제 완료!!

  • Loading branch information...
1 parent 7fb03ee commit 879a7e93a8e47b3bd6cf8787526669bd93669f04 @am0c committed May 24, 2012
Showing with 28 additions and 4 deletions.
  1. +1 −0 index.tex
  2. +8 −2 section-3.R
  3. +19 −2 section-3.tex
View
@@ -5,6 +5,7 @@
\usepackage{verbatim,fancyvrb}
\usepackage{parskip}
\usepackage{graphicx,caption,subcaption}
+\usepackage{amsmath}
\begin{document}
View
@@ -22,9 +22,15 @@ resist <- data.frame(
sink("section-rout-3.txt");
ret <- var.test(resist$val[resist$case == 0], resist$val[resist$case == 1],
- alternative = c("two.sided"), conf.level = 0.95)
+ alternative = c("two.sided"), conf.level = 0.95)
print(ret)
+if (ret$estimate[[1]] == 1) {
+ is_equal = TRUE
+} else {
+ is_equal = FALSE
+}
+
ret <- t.test(resist$val[resist$case == 0], resist$val[resist$case == 1],
- alternative = c("two.sided"), mu = 0, var.equal = FALSE)
+ alternative = c("two.sided"), mu = 0, var.equal = is_equal)
print(ret)
View
@@ -24,11 +24,28 @@ \section{통계적 가설검정}
두 표본의 자료가 30개를 넘기기 때문에 두 모집단이 정규분포라고 가정할 수 있다.
-$H_0: u_1 - u_2, = 0, H_1: u_1 - u_2 \not= 0$
+\begin{align*}
+ H_0: u_1 - u_2, = 0, H_1: u_1 - u_2 \not= 0
+\end{align*}
두 모집단이 정규분포라고 말할 수 있으므로 $F$-검정을 하여 두 모분산도 같은지 추정한다.
$p$-값이 유의수준인 0.05보다 훨씬 작으므로 $H_1$은 기각되지 않는다.
즉, 두 모집단의 표본평균은 유의수준 $5\%$ 내에서 같다고 할 수 없다.
이를 베렌스-피셔 문제라고 한다.
-이를 고려하여 $t$-검정을 하면
+따라서 Welch 또는 새터스웨이터(Satterthwate) 방법을 사용해야한다.
+R에서는 \texttt{t.test} 함수의 인자에 \texttt{var.equal = FALSE}라고 인자를 지정해주어야 한다.
+이렇게 $t$-검정을 한 결과를 보면 $p$-값이 $1.264e-10$으로 유의수준인 0.05보다 훨씬 작기 때문에,
+`전선 1'의 모평균이 `전선 2'의 모평균과 같지 않다는 대립가설을 기각하지 못한다.
+
+따라서 `전선 1'의 모평균과 `전선 2'의 모평균은 유의수준 0.05에서 다르다고 말할 수 있다.
+
+\subsection{스크립트}
+아래와 같이 R 스크립트를 작성했다.
+
+\VerbatimInput{section-3.R}
+
+
+
+
+

0 comments on commit 879a7e9

Please sign in to comment.