Apache Maths (何でも) を使用してクラスをコーディングしましたが、奇妙な問題が発生しました (または、Java コードに問題があるのでしょうか?)。
ここにコードがあります
package distributions;
import java.util.HashMap;
import java.util.Map;
import java.util.Random;
import org.apache.commons.math3.distribution.AbstractIntegerDistribution;
import org.apache.commons.math3.util.FastMath;
public class CustomDistribution extends AbstractIntegerDistribution {
long[] v;
double[] p;
double[] cp;
Map<Long, Double> mp = new HashMap<Long, Double>();
Map<Long, Double> mcp = new HashMap<Long, Double>();
public CustomDistribution (long[] values, double[] probabilities) {
v = values;
p = probabilities;
int len = values.length;
for (int i = 0; i < len; i++) {
mp.put(v[i], p[i]);
}
cp = new double[len];
cp[0] = p[0];
for (int i = 1; i < len; i++) {
cp[i] = cp[i-1] + p[i];
}
for (int i = 0; i < len; i++) {
mcp.put(v[i], cp[i]);
}
}
public Map<Long, Double> getMCP () {return mcp;}
@Override
public double cumulativeProbability(int v) {
return mcp.get(v);
}
// Not implemented, can be done in Clojure code
@Override
public double getNumericalMean() {
return 0;
}
// Not implemented, can be done in Clojure code
@Override
public double getNumericalVariance() {
return 0;
}
@Override
public int getSupportLowerBound() {
return 0;
}
@Override
public int getSupportUpperBound() {
return 0;
}
@Override
public boolean isSupportConnected() {
return false;
}
@Override
public double probability(int v) {
return mp.get(v);
}
// Uses a naive search implementation, should be ok due to data size
public int sample () {
double r = FastMath.random();
int len = p.length;
Boolean flag = false;
int i = 0;
int result = -1;
while (i < len && flag == false) {
if (cp[i] < r) {
i = i + 1;
}
else {
result = (int) v[i];
flag = true;
}
}
return result;
}
}
Clojureラッパー:
(defn create-distribution
"Creates a distribution from empirical data"
[data]
(let [values (long-array (data :values))
probabilities (double-array (data :probabilities))]
(CustomDistribution. values probabilities)))
(create-distribution {:values [1 2 3 4] :probabilities [0.3 0.2 0.2 0.3]})
基本的に、HashMaps と Arrays の 2 つの形式で内部にデータを含むディストリビューションです。
Clojureで、私は試しました:
(.probability dist 4)
(.probability dist (int 4))
どちらも NullPointerException、NativeMethodAccessorImpl..blabla を返します .cumulativeProbabilityMethod については同じです (一方、.sample は正常に動作するため、HashMap の問題である可能性があります)
おそらくコンストラクターで mcp と mp (HashMaps) がコンストラクターで適切に計算されていないことがわかりましたが、試してみると:
(.get (.getMCP dist) 4)
それは私に正しい結果を返します。MCP ハッシュマップは実際、私が期待した結果です。
では、なぜ私のコードが機能しないのですか? 私の2行は「機能的に」同じです。タイプミスをした場合、表示されません。多分それは不変性のせいですか?
ありがとう !