/usr/lib/python2.7/dist-packages/pymc/examples/gp/cov.py is in python-pymc 2.2+ds-1.1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 | from pymc.gp import *
from pymc.gp.cov_funs import matern
from numpy import *
C = Covariance(eval_fun = matern.euclidean, diff_degree = 1.4, amp = .4, scale = 1.)
# C = Covariance(eval_fun = matern.euclidean, diff_degree = 1.4, amp = .4, scale = 1., rank_limit=100)
# C = FullRankCovariance(eval_fun = matern.euclidean, diff_degree = 1.4, amp = .4, scale = 1.)
# C = NearlyFullRankCovariance(eval_fun = matern.euclidean, diff_degree = 1.4, amp = .4, scale = 1.)
#### - Plot - ####
if __name__ == '__main__':
from pylab import *
x=arange(-1.,1.,.01)
clf()
# Plot the covariance function
subplot(1,2,1)
contourf(x,x,C(x,x).view(ndarray),origin='lower',extent=(-1.,1.,-1.,1.),cmap=cm.bone)
xlabel('x')
ylabel('y')
title('C(x,y)')
axis('tight')
colorbar()
# Plot a slice of the covariance function
subplot(1,2,2)
plot(x,C(x,0).view(ndarray).ravel(),'k-')
xlabel('x')
ylabel('C(x,0)')
title('A slice of C')
# show()
|