# 곤충 스프레이 예제
# IV: 6종류의 살충제, DV: 각 살충제 이후 곤충의 수
>>> data(InsectSprays)     # 예제 데이터를 불러오기
>>> attach(InsectSprays)   # 예제 데이터의 변수를 편하게 사용하기 위해 데이터프레임을 부착
>>> str(InsectSprays)      # 예제 데이터의 구조를 확인
'data.frame': 72 obs. of 2 variables:
 $ count: num 10 7 20 14 14 12 10 23 17 20 ...
 $ spray: Factor w/ 6 levels "A","B","C","D",..: 1 1 1 1 1 1 1 1 1 1 ...

>>> mean(InsectSprays$count)
[1] 9.5
>>> var(InsectSprays$count)
[1] 51.88732
>>> table(InsectSprays$spray)
A B C D E F
12 12 12 12 12 12
>>> mean(InsectSprays[InsectSprays$spray=="A", 1])
[1] 14.5

# 살충제 종류별 평균의 차이가 있는지를 보기 위해 일원 분산분석을 이용
>>> attach(InsectSprays)
# 살충제 자료는 등분산을 갖고 있음을 var.equal=TRUE를 통해 표현
>>> oneway.test(count~spray, var.equal=TRUE)
    One-way analysis of means

data: count and spray
F = 34.702, num df = 5, denom df = 66, p-value < 2.2e-16