How to do modeling for seismological receiver function data sets

Hiii…friends, I am new to the pymc3. My aim is to do the modeling using seismological data sets for structure purposes. I actually use receiver function technique for deleneating the subsurface structure. Can it be done using pymc3.please suggest the proper documents for this.Thanks