methods of hrdu (pdu) population estimates
DESCRIPTION
Methods of HRDU (PDU) population estimates. Bruno Sopko , Ph.D. 11 th June 2014. Content. Capture-recapture method Multiplier method (overview) RDS method. Two population study. Two population study. Small sample (Chapman estimator):. Two population – non coded clients. - PowerPoint PPT PresentationTRANSCRIPT
![Page 1: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/1.jpg)
Methods of HRDU (PDU) population estimates
Bruno Sopko, Ph.D. 11th June 2014
![Page 2: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/2.jpg)
Content Capture-recapture method Multiplier method (overview) RDS method
![Page 3: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/3.jpg)
Two population study
z
X Y
xy
![Page 4: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/4.jpg)
Two population study
Small sample (Chapman estimator):
![Page 5: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/5.jpg)
Two population – non coded clients
![Page 6: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/6.jpg)
Three and more samples R – statistical program R add-ons (Rcapture) Log-linear estimation (Poisson, Chao, AIC) Data preparation
![Page 7: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/7.jpg)
Robust design
![Page 8: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/8.jpg)
R-project http://www.r-project.org/– http://cran.rstudio.com/– http://cran.rstudio.com/manuals.html
Required packages:– Rcapture– foreign
![Page 9: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/9.jpg)
![Page 10: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/10.jpg)
Rcapture installation
![Page 11: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/11.jpg)
Data preparationSubst 2010 Needle 2010 Subst 2011 Needle 2011 Subst 2012 Frequency
[1,] 1 1 1 1 1 1520
[2,] 1 1 1 1 0 0
[3,] 1 1 1 0 1 0
[4,] 1 1 1 0 0 0
[5,] 1 1 0 1 1 0
[6,] 1 1 0 1 0 0
[7,] 1 1 0 0 1 0
[8,] 1 1 0 0 0 0
[9,] 1 0 1 1 1 0
[10,] 1 0 1 1 0 0
[11,] 1 0 1 0 1 0
[12,] 1 0 1 0 0 0
[13,] 1 0 0 1 1 0
[14,] 1 0 0 1 0 0
[15,] 1 0 0 0 1 0
[16,] 1 0 0 0 0 0
[17,] 0 1 1 1 1 284
[18,] 0 1 1 1 0 0
[19,] 0 1 1 0 1 0
[20,] 0 1 1 0 0 0
[21,] 0 1 0 1 1 10
[22,] 0 1 0 1 0 0
[23,] 0 1 0 0 1 205
[24,] 0 1 0 0 0 667
[25,] 0 0 1 1 1 0
[26,] 0 0 1 1 0 0
[27,] 0 0 1 0 1 0
[28,] 0 0 1 0 0 0
[29,] 0 0 0 1 1 0
[30,] 0 0 0 1 0 0
[31,] 0 0 0 0 1 0
Total 2686
57 % use all services
These data were wronglyconverted
![Page 12: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/12.jpg)
Results Closed population model for every period: M0
Model fit: deviance df AICfitted model 3974.074 11 4011.492
Capture probabilities: estimate stderrperiod 1 0.9466 0.0028period 2 1.0000 0.0000
Survival probabilities: estimate stderrperiod 1 -> 2 0.6632 0.0092
Abundances: estimate stderrperiod 1 2735.2 7.7period 2 1814.0 0.1
Number of new arrivals: estimate stderrperiod 1 -> 2 0 0
Total number of units who ever inhabited the survey area: estimate stderrall periods 2735.2 7.7
Total number of captured units: 2686
![Page 13: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/13.jpg)
Multiplier method Network size reliability Statistical weight Standard error of weighted mean – SPSS x “rest
of the world”
![Page 14: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/14.jpg)
Response driven sampling MethodRespondent-driven sampling (RDS), combines “snowball sampling" (getting individuals to refer those they know, these individuals in turn refer those they know and so on) with a mathematical model that weights the sample to compensate for the fact that the sample was collected in a non-random way.
![Page 15: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/15.jpg)
RDS – essential information Personal Network Size (Degree) - Number of
people the respondent knows within the target population.
Respondent's Serial Number - Serial number of the coupon the respondent was recruited with.
Respondent's Recruiting Serial Numbers - Serial numbers from the coupons the respondent is given to recruit others.
![Page 16: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/16.jpg)
RDS - program http://
wiki.stat.ucla.edu/hpmrg/index.php/RDS_Analyst_Install
![Page 17: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/17.jpg)
RDS Analyst
![Page 18: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/18.jpg)
RDS computing – data formatIdentifikacioni broj Mreza Kupon1 Kupon2 Kupon3 Pol Godina rodjenja Godine
010000000000 20 011000000000 012000000000 013000000000 m 1978 35
030000000000 10 031000000000 032000000000 033000000000 m 1989 24
020000000000 29 021000000000 022000000000 023000000000 ž 1990 23
050000000000 30 051000000000 052000000000 053000000000 m 1985 28
040000000000 19 041000000000 042000000000 043000000000 m 1991 22
060000000000 6 061000000000 062000000000 063000000000 m 1978 35
070000000000 30 071000000000 072000000000 073000000000 ž 1978 35
032000000000 7 032100000000 032200000000 032300000000 m 1980 33
011000000000 20 011100000000 011200000000 011300000000 m 1982 31
021000000000 20 021100000000 021200000000 021300000000 m 1984 29
012000000000 10 012100000000 012200000000 012300000000 m 1985 28
013000000000 50 013100000000 013200000000 013300000000 m 1978 35
022000000000 50 022100000000 022200000000 022300000000 m 1981 32
053000000000 20 053100000000 053200000000 053300000000 m 1978 35
![Page 19: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/19.jpg)
Video
![Page 20: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/20.jpg)
Results – recruitment tree
![Page 21: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/21.jpg)
![Page 22: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/22.jpg)
RDS – convergence (theory)
![Page 23: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/23.jpg)
RDS - result
Mean Median Mode 90% 3% 98% Prior 1935 1284 600 4422 766 2444 Posterior 360 359 359 384 325 401
Summary of Population Size Estimation
![Page 24: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/24.jpg)
RDS – convergence results
![Page 25: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/25.jpg)
Thank you for your attention
![Page 26: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/26.jpg)
Corrected CRM – preliminary bias observation
The codes supplied were from Syringe exchange program and from Substitution treatment program
These programs are mutually exclusive – therefore negatively biased
The data were collected over three years, therefore data from the same program have been positively biased
The data from Syringe exchange program for year 2012 were incomplete
![Page 27: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/27.jpg)
Corrected CRM - method Due to the previously described bias problems,
the classical robust design method (robustd.0 or robustd.t functions in R) could not be deployed.
The capture-recapture method with bias correction over all years have been used (closedp.mX function in R), the total number of estimated PDU has been broken into years values by the corresponding code counts.
![Page 28: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/28.jpg)
Corrected CRM - dataNumber of captured units: 6782 Frequency statistics: fi i = 1 3937i = 2 1360i = 3 1387i = 4 51i = 5 32i = 6 15fi: number of units captured i times
![Page 29: Methods of HRDU (PDU) population estimates](https://reader036.vdocument.in/reader036/viewer/2022062321/5681305c550346895d961dcf/html5/thumbnails/29.jpg)
Corrected CRM - results
Year Population STD Error Arrived Left Stayed STD Error Stayed all years STD Error
2010 24169 1744
2011 20790 1500 8799 12178 11911 3244
2012 19813 1429 8817 7894 11919 2929 6355 4673