Has anyone here done work on modelling order statistics within PyMC3? Suppose my data comes in a set of rankings (or real numbers) for n items X_1 < X_2 < … < X_n over multiple observations and im interested in determining the distribution of the ranking itself such as: given item i, what is the …

OK, I cooked up something that kind of work: import numpy as np import pymc3 as pm import theano.tensor as tt # data K = 5 # number of items being ranked J = 20 # number of raters yreal = np.argsort(np.random.randn(K, 1), axis=0) print(yreal) y = np.argsort(yreal + np.random.randn(K, J), axis=0) p…

Order statistics in PyMC3

Questions

dycontri June 21, 2021, 6:32pm 8

I’m having trouble understanding how the model, which doesn’t have any observed RVs, is able to fit to the data

Can you help me understand this, @junpenglao ?

Topic		Replies	Views
Error using Ordered() in Mixture of Uniform Distributions Questions	14	1899	February 7, 2019
Sampling uniformly in a triangular support Questions	15	2723	May 3, 2018
Mixture models and breaking class symmetry Questions	11	3932	October 12, 2018
Why does `transform=pm.distributions.transforms.ordered` lead to worse convergence? Questions	3	1579	August 27, 2021
Sampling ordered normal distributions with array-valued bounds Questions	0	347	December 7, 2020

Order statistics in PyMC3

Related topics