#Typical call is %srsmean 20 c1 c10-c12. c10 and c11 contain a table
#displaying the estimators and c12 contains the 20 points sampled 
#from c1
macro
srsmean samsize popcol labcol descol samcol

mcolumn popcol labcol descol samcol
mconstant samsize sammean samsd samse popsize



sample samsize popcol samcol
let popsize=count(popcol)
let sammean=mean(samcol)
let samsd=stan(samcol)
let samse=(samsd/sqrt(samsize))*sqrt((popsize-samsize)/popsize)
set labcol;
format(a30).
Population size
Sample size
Sample mean
Sample mean standard error
end
set descol
popsize
samsize
sammean
samse
end
name samcol 'Sample' 

endmacro