|
| |
LAB 2 - 9/15/97
Tomato Experiment: two sample t-test.
> fertA<-c(29.9,11.4,25.3,16.5,21.1)
> fertA
[1] 29.9 11.4 25.3 16.5 21.1
> fertB_c(26.6,23.7,28.5,14.2,17.9,24.3)
> fertB
[1] 26.6 23.7 28.5 14.2 17.9 24.3
> mean(fertA)
[1] 20.84
> mean(fertB)
[1] 22.53333
> var(fertA)
[1] 52.498
> var(fertB)
[1] 29.50667
> s2<-(var(fertA)*4+var(fertB)*5)/(5+6-1-1)
> s2
[1] 39.72504
> obst_(mean(fertB)-mean(fertA))/(sqrt(s2)*sqrt((1/5)+(1/6)))
> obst
[1] 0.4436848
> 1-pt(.44,9)
[1] 0.335157
> 2*(1-pt(.44,9))
[1] 0.6703141
> gamma(11+1)/(gamma(5+1)*gamma(6+1))
[1] 462
> sample(11,5)
[1] 2 5 11 4 6
> sample(11,5)
[1] 7 11 6 1 5
> tomato_c(fertA,fertB)
> fert<-c(1,1,1,1,1,2,2,2,2,2,2)
> cbind(fert,tomato)
fert tomato
[1,] 1 29.9
[2,] 1 11.4
[3,] 1 25.3
[4,] 1 16.5
[5,] 1 21.1
[6,] 2 26.6
[7,] 2 23.7
[8,] 2 28.5
[9,] 2 14.2
[10,] 2 17.9
[11,] 2 24.3
> motif()
> plot(fert,tomato)
> text(1.2,15,"A")
> text(1.8,15,"B")
Pulp Experiment.
Question: Do two sample t-test for any two operators.
Also computer 95% confidence interval for difference in mean between the
same two operators.
You have to turn in output by the end of the lab.
Output will be graded.
Vadim Kutsyy
kutsyy@umich.edu
|