87
|
1 #!/usr/bin/env python
|
|
2
|
|
3 """
|
|
4 generate a Gaussian dataset
|
|
5 """
|
|
6
|
|
7 import random
|
|
8 import sys
|
|
9 from .cli import DatasetGenerationParser
|
|
10
|
|
11
|
|
12 class GaussianSampler(object):
|
|
13
|
|
14 def __init__(self, D, sigma=1., center=None):
|
|
15 assert D > -1
|
|
16 self.D = D
|
|
17 if center is None:
|
|
18 self.center = [0]*self.D
|
|
19 else:
|
|
20 assert len(center) == self.D
|
|
21 self.center = center
|
|
22 self.sigma = sigma
|
|
23
|
|
24 def __call__(self):
|
|
25 """yield one point of a Gaussian distribution"""
|
|
26
|
|
27 return [random.gauss(x, self.sigma)
|
|
28 for x in self.center]
|
|
29
|
|
30
|
|
31 def main(args=sys.argv[1:]):
|
|
32 """CLI"""
|
|
33
|
|
34 # parse command line
|
|
35 parser = DatasetGenerationParser(description=__doc__)
|
|
36 options = parser.parse_args(args)
|
|
37
|
|
38 # point generator
|
|
39 generator = GaussianSampler(D=options.dimensions)
|
|
40 points = [generator() for _ in xrange(options.number)]
|
|
41
|
|
42 # output
|
|
43 parser.writer().writerows(points)
|
|
44
|
|
45 if __name__ == '__main__':
|
|
46 main()
|