# Scaling covariance function in Gaussian Processes
**URL:** https://discourse.pymc.io/t/scaling-covariance-function-in-gaussian-processes/4013
**Category:** Questions
**Created:** [October 24, 2019, 9:26am UTC](https://discourse.pymc.io/t/scaling-covariance-function-in-gaussian-processes/4013 "2019-10-24T09:26:30Z")
**Posts on this page:** 4
**Page:** 1
### Author:  [@Mcb](https://discourse.pymc.io/u/Mcb)
#### Post date: [October 24, 2019, 9:26am UTC](https://discourse.pymc.io/t/scaling-covariance-function-in-gaussian-processes/4013/1 "2019-10-24T09:26:30Z")
Hi,
This one is directed to @bwengals since it’s a direct question on one of his [notebooks](https://bwengals.github.io/gaussian-processes-models.html).
In particular there is a scaling RV σsq\_gp which is multiplied by the overall covariance function.
Can you expand on the choice for this variable and why it was used?
Below there is additional context for the model.
```
# define positive normal distribution with a lower bound
BoundedNormal = pm.Bound(pm.Normal, lower=0)
# define gp model in pymc3
with pm.Model() as conjugate_gp:
σsq_gp = pm.HalfCauchy("σsq_gp", beta=2) # scales the covariance function
ℓ = BoundedNormal("ℓ", mu=0.5, sd=1) # lengthscale of the covariance function (weakly informative)
k = σsq_gp * pm.gp.cov.Matern52(1, ℓ) # the covariance function
....
```
Thanks!
---
### Author:  [@bwengals](https://discourse.pymc.io/u/bwengals)
#### Post date: [October 24, 2019, 4:17pm UTC](https://discourse.pymc.io/t/scaling-covariance-function-in-gaussian-processes/4013/2 "2019-10-24T16:17:33Z")
That parameters role is to scale the covariance function. It controls the magnitude of the GP. You can try setting up a GP prior with different values of that parameter and see what the function samples look like to see how it affects things.
I’d want to refer you to the pymc3 documentation and examples, the code on that blog is extremely out of date.
---
### Author:  [@Mcb](https://discourse.pymc.io/u/Mcb)
#### Post date: [October 28, 2019, 9:45am UTC](https://discourse.pymc.io/t/scaling-covariance-function-in-gaussian-processes/4013/3 "2019-10-28T09:45:37Z")
Thanks for the input.
I’m somewhat familiar with the examples on the pymc3 docs. Great stuff by the way.
I was wondering if the scaling factor was following a formal rule to control the GP.
---
### Author:  [@bwengals](https://discourse.pymc.io/u/bwengals)
#### Post date: [October 28, 2019, 2:48pm UTC](https://discourse.pymc.io/t/scaling-covariance-function-in-gaussian-processes/4013/4 "2019-10-28T14:48:33Z")
Yes it is. It comes from the rules for constructing kernels. Multiplying a kernel (or covariance function) by some function a(x) is a valid kernel: a(x) k(x, x') a(x'), so to design a changepoint a(x) you can use a sigmoid function for a(x).