Я столкнулся с ошибкой при построении CDF для треугольного распределения; Я могу построить гистограмму, используя Создание треугольного распределения в MATLAB , но как сделать Я строю CDF? n = 10000; %Random number a = 0.26; %Min b = 0.46;...

Я не могу понять, почему мой простой код, приведенный ниже, не работает. Я получаю NaN для значений P_j . % Plotting the CDF for normally distributed variable mu_j = 0.008; simga_j = -0.002; x= 0.008 -0.002*randn(100000,1); X_j = min(x) :...

У меня есть заданная функция (назовем ее f (x)), и я использовал метод Монте-Карло, чтобы нормализовать ее. Я вычислил функцию плотности вероятности и после ее интегрирования получил кумулятивную функцию распределения. f = function(x) ......

Я пытаюсь заставить Rcpp-версию pmvnorm работать по крайней мере так же быстро, как mvtnorm :: pmvnorm в R. Я нашел https://github.com/zhanxw/libMvtnorm и создал скелетный пакет Rcpp с соответствующие исходные файлы. Я добавил следующие...

Итак, я перевожу Matlab в код Python и наткнулся на эту строку Pv = chi2cdf(Devk,Md,'upper'); Где Devk — трехмерная матрица, а Md — df=10. Я искал эквивалент в Python. Я нашел это: from scipy.stats import chi2 Pv = chi2.cdf(Devk,Md)...

Я пытаюсь построить распределение, используя scipy.stats.pearson3, но я получаю противоречивые результаты при просмотре кумулятивной функции распределения (или функции выживания). В качестве примера: scipy.stats.pearson3.cdf(-1, skew=0, loc=0,...

Я использую пакет ks для оценки плотности ядра. Вот простой пример: n <- 70 x <- rnorm(n) library(ks) f_kde <- kde(x) На самом деле меня интересуют соответствующие вероятности превышения моих входных данных, которые могут быть...