I appreciate your help! I use pm.Uniform.dist(lower=0, upper=0) to specify a mass point in my distribution. I applied the suggested function and found out this error message is because there are too many 0 in my observation. I solved this problem by adding 0.0000001 to all the observations with a value of 0. Thank you so much for your help again, and the suggested function is very useful!