For a Bayesian linear regression problem like the following basic_model = pm.Model() with basic_model: # Priors for unknown model parameters alpha = pm.Normal("alpha", mu=0, sigma=10) beta = pm.Normal("beta", mu=0, sigma=10, shape=2) sigma = pm.HalfNormal("sigma", sigma=1) # …

Aesara tensors have a .dot method that lets you easily do matrix multiplications: mu = alpha + X.dot(beta) So will need to make sure your design matrix (X) is an aesara tensor. Using pm.Data is probably the cleanest way.

Write in matrix multiplication form

Questions version agnostic

jackshi June 19, 2022, 2:47am 3

Really helpful, thanks a lot!!

Topic		Replies	Views
Matrix multiplication in pymc3 Questions development	8	6389	April 24, 2018
Create model matrix Questions	7	1518	June 23, 2021
Matrix \mu in multivariate normal Questions	2	594	May 10, 2019
Different beta functions in a model Questions	4	486	January 17, 2021
Matrix Multiplication With Multiple Dimensions in PYMC Model v5 modeling	4	1114	July 7, 2022

Write in matrix multiplication form

Related topics