View Javadoc
1   /*
2    * Licensed to the Apache Software Foundation (ASF) under one or more
3    * contributor license agreements.  See the NOTICE file distributed with
4    * this work for additional information regarding copyright ownership.
5    * The ASF licenses this file to You under the Apache License, Version 2.0
6    * (the "License"); you may not use this file except in compliance with
7    * the License.  You may obtain a copy of the License at
8    *
9    *      http://www.apache.org/licenses/LICENSE-2.0
10   *
11   * Unless required by applicable law or agreed to in writing, software
12   * distributed under the License is distributed on an "AS IS" BASIS,
13   * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14   * See the License for the specific language governing permissions and
15   * limitations under the License.
16   */
17  package org.apache.commons.statistics.examples.distribution;
18  
19  import java.util.ArrayList;
20  import java.util.List;
21  import org.apache.commons.statistics.distribution.DiscreteDistribution;
22  import org.apache.commons.statistics.distribution.ZipfDistribution;
23  import picocli.CommandLine.ArgGroup;
24  import picocli.CommandLine.Command;
25  import picocli.CommandLine.Option;
26  
27  /**
28   * Command for the {@link ZipfDistribution}.
29   */
30  @Command(name = "zipf",
31           description = "Zipf distribution.",
32           subcommands = {
33               ZipfCommand.Check.class,
34               ZipfCommand.PMF.class,
35               ZipfCommand.LPMF.class,
36               ZipfCommand.CDF.class,
37               ZipfCommand.SF.class,
38               ZipfCommand.ICDF.class,
39               ZipfCommand.ISF.class,
40           })
41  class ZipfCommand extends AbstractDistributionCommand {
42  
43      /** Base command for the distribution that defines the parameters. */
44      private abstract static class BaseCommand extends AbstractDiscreteDistributionCommand {
45          /** Distribution parameters. */
46          @ArgGroup(validate = false, heading = "Distribution parameters:%n", order = 1)
47          private Params params = new Params();
48  
49          /** Parameters class. */
50          static class Params {
51              /** The distribution number of elements. */
52              @Option(names = {"-n", "--number-of-elements"},
53                      arity = "1..*",
54                      split = ",",
55                      description = {"number of elements (default: ${DEFAULT-VALUE})."})
56              private int[] n = {10};
57  
58              /** The distribution exponent. */
59              @Option(names = {"-e", "--exponent"},
60                      arity = "1..*",
61                      split = ",",
62                      description = {"exponent (default: ${DEFAULT-VALUE})."})
63              private double[] e = {1, 2, 3, 4};
64          }
65  
66          /** Extend the options to set the default values for this distribution. */
67          static final class Options extends DiscreteDistributionOptions {
68              /** Set defaults. */
69              private Options() {
70                  min = 1;
71                  max = 10;
72              }
73          }
74  
75          @Override
76          protected List<Distribution<DiscreteDistribution>> getDistributions() {
77              int[] n = params.n;
78              double[] e = params.e;
79              final int max = DistributionUtils.validateLengths(n.length, e.length);
80  
81              n = DistributionUtils.expandToLength(n, max);
82              e = DistributionUtils.expandToLength(e, max);
83  
84              // Create distributions
85              final ArrayList<Distribution<DiscreteDistribution>> list = new ArrayList<>();
86              for (int i = 0; i < max; i++) {
87                  final DiscreteDistribution d = ZipfDistribution.of(n[i], e[i]);
88                  list.add(new Distribution<>(d, "n=" + n[i] + ",p=" + e[i]));
89              }
90              return list;
91          }
92      }
93  
94      /** Base command for the distribution that defines the parameters. */
95      private abstract static class ProbabilityCommand extends BaseCommand {
96          /** The distribution options. */
97          @ArgGroup(validate = false, heading = "Evaluation options:%n", order = 2)
98          private Options distributionOptions = new Options();
99  
100         @Override
101         protected DistributionOptions getDistributionOptions() {
102             return distributionOptions;
103         }
104     }
105 
106     /** Base command for the distribution that defines the parameters for inverse probability functions. */
107     private abstract static class InverseProbabilityCommand extends BaseCommand {
108         /** The distribution options. */
109         @ArgGroup(validate = false, heading = "Evaluation options:%n", order = 2)
110         private InverseDiscreteDistributionOptions distributionOptions = new InverseDiscreteDistributionOptions();
111 
112         @Override
113         protected DistributionOptions getDistributionOptions() {
114             return distributionOptions;
115         }
116     }
117 
118     /** Verification checks command. */
119     @Command(name = "check",
120              hidden = true,
121              description = "Zipf distribution verification checks.")
122     static class Check extends ProbabilityCommand {}
123 
124     /** PMF command. */
125     @Command(name = "pmf",
126              aliases = {"pdf"},
127              description = "Zipf distribution PMF.")
128     static class PMF extends ProbabilityCommand {}
129 
130     /** LPMF command. */
131     @Command(name = "lpmf",
132              aliases = {"lpdf"},
133              description = "Zipf distribution natural logarithm of the PMF.")
134     static class LPMF extends ProbabilityCommand {}
135 
136     /** CDF command. */
137     @Command(name = "cdf",
138              description = "Zipf distribution CDF.")
139     static class CDF extends ProbabilityCommand {}
140 
141     /** SF command. */
142     @Command(name = "sf",
143              description = "Zipf distribution survival probability.")
144     static class SF extends ProbabilityCommand {}
145 
146     /** ICDF command. */
147     @Command(name = "icdf",
148              description = "Zipf distribution inverse CDF.")
149     static class ICDF extends InverseProbabilityCommand {}
150 
151     /** ISF command. */
152     @Command(name = "isf",
153              description = "Zipf distribution inverse SF.")
154     static class ISF extends InverseProbabilityCommand {}
155 }