datamining hw2

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

ttest= function(X1,X2,delta,type){
n1 = length(X1)
n2 = length(X2)
barX1 = mean(X1)
barX2 = mean(X2)
s1 = sd(X1)
s2 = sd(X2)

T = ((barX1 - barX2) - delta)/sqrt(((s1)^2/n1)+((s2)^2/n2))

if (type == "two.sided") {
  cat("H0 : mu1 - mu2 = " , delta , "\n" , "H1 : mu1 - mu2 !=" , delta , "\n")
    cat("검정통계량 T = " , T , "\n" , "유의확률 = ", 1-(pnorm(abs(T)) - pnorm(-abs(T))) , "\n")
    if (abs(T) >= qnorm(0.975)){
        cat("결론 : H0를 기각할 수 있습니다. \n")}
    else {cat("결론 : H0를 기각할 수 없습니다. \n")}   
}

if (type == "less") {
    cat("H0 : mu1 - mu2 = " , delta , "\n" , "H1 : mu1 - mu2 < " , delta , "\n")
    cat("검정통계량 T = " , T , "\n" , "유의확률 = " , pnorm(T), "\n")
    if (T <= -qnorm(0.95)) {
      cat("결론 : H0를 기각할 수 있습니다. \n")}
    else {cat ("결론 : H0를 기각할 수 없습니다. \n")}

}

 if (type == "greater"){
    cat(" H0 : mu1 - mu2 =" , delta , "\n" , "H1 : mu1 - mu2 >" , delta , "\n")
    cat("검정통계량 T =" , T , "유의확률 =" , 1-pnorm(T) , "\n")
    if (T >= qnorm(0.95)) {
      cat("결론 : H0를 기각할 수 있습니다. \n")
    }
    else {cat ("결론 : H0를 기각할 수 없습니다. \n")}
  }
}
##함수 예##
X1 = rnorm(30,0,1)
X2 = rnorm(30,1,2)
de.2 = 0

ttest(X1,X2,de.2,"two.sided")

## H0 : mu1 - mu2 =  0 
##  H1 : mu1 - mu2 != 0 
## 검정통계량 T =  -2.886 
##  유의확률 =  0.003906 
## 결론 : H0를 기각할 수 있습니다.