Agriculture Reference
In-Depth Information
q1obs
1
2
3
q2obs
1
total
50.00000 70.00000
50.00000
SE
20.78024 24.32731
20.78024
2
total
40.00000 70.00000 100.00000
SE
18.68397 24.32731
28.60388
3
total
70.00000 90.00000
80.00000
SE
24.32731 27.28636
25.86679
4
total
40.00000 80.00000 110.00000
SE
18.68397 25.86679
29.83287
5
total
70.00000 40.00000
40.00000
SE
24.32731 18.68397
18.68397
Two more useful functions are subset , which transforms a design object by
extracting the part that satisfies a given rule, and svyby , which repeats the
estimation process for every code of a vector. These two actions are very useful
in domain estimation (see Chap. 11 ).
> esrs0 < - subset(dsrs,q1obs ¼¼ 1)
> domesrs0 < - svytotal(~yobs,esrs0)
> domesrs0
total SE
yobs 25647 4134.4
> domesrs1 < - svyby(~yobs,~q1obs,dsrs,svytotal)
> domesrs1
q1obs
yobs
se
1
1 25647.31 4134.426
2
2 30463.50 4211.684
3
3 35555.94 4483.966
> domesrs2 < - svyby(~yobs,~q1obs+q2obs,dsrs,svytotal)
> domesrs2
q1obs q2obs
yobs
se
1.1
1
1 3259.115 1358.853
2.1
2
1 2992.660 1073.066
3.1
3
1 2846.868 1211.068
1.2
1
2 3263.294 1525.462
2.2
2
2 5703.026 1984.628
3.2
3
2 8123.017 2329.283
1.3
1
3 6586.044 2292.237
2.3
2
3 8525.725 2586.441
3.3
3
3 7835.151 2534.085
1.4
1
4 4164.286 1945.902
2.4
2
4 8412.343 2722.331
Search WWH ::




Custom Search