***～会飞的蝎子～学习宁静致远 » 日志 » HW8Q9 UseR
HW8Q9 UseR
蝎子 发表于 2008-03-08 09:34:35
mid=c(0.75,0.65,0.55,0.45,0.35,0.25,0.15,0.05,-0.05,-0.15,-0.25,-0.35,-0.45,-0.55,-0.65,-0.75)
> mid
[1] 0.75 0.65 0.55 0.45 0.35 0.25 0.15 0.05 -0.05 -0.15 -0.25 -0.35 -0.45 -0.55 -0.65 -0.75
> length(mid)
[1] 16
> freq=c(2,3,7,13,13,19,22,24,26,19,17,12,11,5,4,3)
> length(freq)
[1] 16
> sum(mid*freq)/200
[1] 0.009
> sum((mid-0.009)^2*freq)/199
[1] 0.1035367
> sqrt(sum((mid-0.009)^2*freq)/199)
[1] 0.3217712
exp=(pnorm((mid-0.05-0.009)/0.3218)-pnorm((mid+0.05-0.009)/0.3218))*200
> exp
[1] -1.780026 -3.450749 -6.078439 -9.728911 -14.149182 -18.697975 -22.451990 -24.496991 -24.286703 -21.878738
[11] -17.909090 -13.320535 -9.002554 -5.528465 -3.084867 -1.564082
> exp=(pnorm((mid+0.05-0.009)/0.3218)-pnorm((mid-0.05-0.009)/0.3218))*200
> exp
[1] 1.780026 3.450749 6.078439 9.728911 14.149182 18.697975 22.451990 24.496991 24.286703 21.878738 17.909090 13.320535
[13] 9.002554 5.528465 3.084867 1.564082
> sum(exp)
[1] 197.4093
> data.frame(mid,freq,exp)
mid freq exp
1 0.75 2 1.780026
2 0.65 3 3.450749
3 0.55 7 6.078439
4 0.45 13 9.728911
5 0.35 13 14.149182
6 0.25 19 18.697975
7 0.15 22 22.451990
8 0.05 24 24.496991
9 -0.05 26 24.286703
10 -0.15 19 21.878738
11 -0.25 17 17.909090
12 -0.35 12 13.320535
13 -0.45 11 9.002554
14 -0.55 5 5.528465
15 -0.65 4 3.084867
16 -0.75 3 1.564082
> sum((freq-exp)^2/exp)
[1] 4.203124
> mid
[1] 0.75 0.65 0.55 0.45 0.35 0.25 0.15 0.05 -0.05 -0.15 -0.25 -0.35 -0.45 -0.55 -0.65 -0.75
> length(mid)
[1] 16
> freq=c(2,3,7,13,13,19,22,24,26,19,17,12,11,5,4,3)
> length(freq)
[1] 16
> sum(mid*freq)/200
[1] 0.009
> sum((mid-0.009)^2*freq)/199
[1] 0.1035367
> sqrt(sum((mid-0.009)^2*freq)/199)
[1] 0.3217712
exp=(pnorm((mid-0.05-0.009)/0.3218)-pnorm((mid+0.05-0.009)/0.3218))*200
> exp
[1] -1.780026 -3.450749 -6.078439 -9.728911 -14.149182 -18.697975 -22.451990 -24.496991 -24.286703 -21.878738
[11] -17.909090 -13.320535 -9.002554 -5.528465 -3.084867 -1.564082
> exp=(pnorm((mid+0.05-0.009)/0.3218)-pnorm((mid-0.05-0.009)/0.3218))*200
> exp
[1] 1.780026 3.450749 6.078439 9.728911 14.149182 18.697975 22.451990 24.496991 24.286703 21.878738 17.909090 13.320535
[13] 9.002554 5.528465 3.084867 1.564082
> sum(exp)
[1] 197.4093
> data.frame(mid,freq,exp)
mid freq exp
1 0.75 2 1.780026
2 0.65 3 3.450749
3 0.55 7 6.078439
4 0.45 13 9.728911
5 0.35 13 14.149182
6 0.25 19 18.697975
7 0.15 22 22.451990
8 0.05 24 24.496991
9 -0.05 26 24.286703
10 -0.15 19 21.878738
11 -0.25 17 17.909090
12 -0.35 12 13.320535
13 -0.45 11 9.002554
14 -0.55 5 5.528465
15 -0.65 4 3.084867
16 -0.75 3 1.564082
> sum((freq-exp)^2/exp)
[1] 4.203124
曾经的这一天...
- » 2007年: TA生活
- » 2006年: next week is called dead week:)
- » 2005年: 爱情和女人的节日……
相关日志:
收藏:
QQ书签
del.icio.us
订阅:
Google
抓虾
