Point estimate covariates branch works kind of

apatil · Apr 23, 2009 · 8a82adb · 8a82adb
1 parent d9f2f94
commit 8a82adb
Showing 1 changed file with 4 additions and 2 deletions.
diff --git a/cov_test/test_marg.py b/cov_test/test_marg.py
@@ -4,8 +4,8 @@
 import pylab as pl
 from marginal_model import make_model
 from util import *
-n_data = 9
-n_pred = 1
+n_data = 90
+n_pred = 10
 
 # make_model(d,lon,lat,t,covariate_values)
 names = ['rain','temp','ndvi']
@@ -51,6 +51,8 @@
 
 prior_var=1e20
 N = pm.MAP(make_model(transform_bin_data(pos[:n_data],neg[:n_data]),lon[:n_data],lat[:n_data],t[:n_data],cv_data,2,prior_var))
+N.fit()
+print N.C.value.params
 
 # C1 = N.T.parents['C'].value(N.logp_mesh, N.logp_mesh)+N.V.value*np.eye(N.logp_mesh.shape[0]) + N.u.T*N.u*prior_var
 # C2 = N.marg_T.value.I