st: R: Is there a better way of allocating individuals into different categories based on the probability of an event?

["Carlo Lazzaro" <carlo.lazzaro@tiscalinet.it>]

Dear Gwinyai,
the first trivial step to make your code running is setting the number of
observations:
-----------------------------------------------------
set obs  100 //or whatever you like or need
gen smoker1=int(1+(50000-1)*runiform())
gen smoker=0 if smoker1 < 42000
replace smoker=1 if smoker1 >= 42000
label define smoker_label 0"Non-smoker" 1"Smoker"
label values smoker smoker_label
label variable smoker "Smoking status"
drop smoker1
-----------------------------------------------------

Joseph has already commented on more substantive details.

Kindest regards,
Carlo
-----Messaggio originale-----
Da: owner-statalist@hsphsun2.harvard.edu
[mailto:owner-statalist@hsphsun2.harvard.edu] Per conto di Gwinyai Masukume
Inviato: sabato 23 marzo 2013 04:16
A: statalist@hsphsun2.harvard.edu
Oggetto: st: Is there a better way of allocating individuals into different
categories based on the probability of an event?

Dear Stata-listers,

Is there a better way of executing the code below which assigns individuals
into the Non-smoker and Smoker categories based on the probability of
smoking in pregnancy?

****0=Non-smoker, 1=Smoker***
**** Assuming 84% of women do not smoke during pregnancy
**** Johnston V, Thomas DP, McDonnell J, Andrews RM. Maternal smoking and
smoking in the household during pregnancy and postpartum: findings from an
Indigenous cohort in the Northern Territory. Med J Aust. 2011 May
16;194(10):556-9.
gen smoker1=int(1+(50000-1)*runiform())
gen smoker=0 if smoker1 < 42000
replace smoker=1 if smoker1 >= 42000
label define smoker_label 0"Non-smoker" 1"Smoker"
label values smoker smoker_label
label variable smoker "Smoking status"
drop smoker1

Kind regards,
Gwinyai
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/

*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/faqs/resources/statalist-faq/
*   http://www.ats.ucla.edu/stat/stata/