#
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
# 
#-----

tot_num_crypt_foci = c(87,53,72,90,78,85,83)

print( mean( tot_num_crypt_foci ) )
print( var( tot_num_crypt_foci ) ) # the mean and the variance of the crypt data are not equal ... not a Poisson process 

x = rpois( 7, 78.3 ) # use the sample mean from the crypt foci data 
print( mean(x) ); print( var(x) ) # the mean and the variance are "equal"