I am trying to learn PyMC by starting off with a simple parameter identification using linear regression. I have a very simple mass-spring-damper model which I integrate numerically to get the displacement, velocity, acceleration and applied force as a function of time (I add noise to simulate measurement error).

My question: Is it possible to set up a PyMC model to identify the stiffness, damping and mass of the system form these noisy time series? (The displacement, velocity and acceleration time series form the regressor matrix, the force time series forms the observed vector, and the unknown parameters are the regression coefficients).

Any suggestions on how to go about solving this problem via PyMC would be greatly appreciated!

… unfortunately it is not up to date. But there might be other useful traces in the v4 examples.

Generally, I find probabilistic models very useful to estimate parameters underlying noisy data, in very simple models. I have a simple measurement documented here, but that was also quite a long time ago

Yeah the search engine indexing for any specific example notebooks seems like a big mess. It’s best to just look up the example gallery first, then look for what you want from there.