Actually you are sampling the pm.free_RVs, which in this case is w_stick_breaking__ and tau_log__ and they dont have a random property/function. But nonetheless for transformed distribution we should be able to also add the random method - the way to go is to sample from the basic_RVs, and use the transformed forward value function to get the sample in the transformed space.