def errvar_test():
import os
from pyemu import ErrVar
w_dir = os.path.join("..", "verification", "henry")
forecasts = ["pd_ten", "c_obs10_2"]
ev = ErrVar(jco=os.path.join(w_dir, "pest.jcb"), forecasts=forecasts)
print(ev.prior_forecast)
print(ev.get_errvar_dataframe())
def errvar_test():
import os
from pyemu import ErrVar
w_dir = os.path.join("..","..","verification","henry")
forecasts = ["pd_ten","c_obs10_2"]
ev = ErrVar(jco=os.path.join(w_dir,"pest.jcb"),forecasts=forecasts)
print(ev.prior_forecast)
print(ev.get_errvar_dataframe())
def errvar_test_nonpest():
import numpy as np
from pyemu import ErrVar, Matrix, Cov
#non-pest
pnames = ["p1","p2","p3"]
onames = ["o1","o2","o3","o4"]
npar = len(pnames)
nobs = len(onames)
j_arr = np.random.random((nobs,npar))
jco = Matrix(x=j_arr,row_names=onames,col_names=pnames)
parcov = Cov(x=np.eye(npar),names=pnames)
obscov = Cov(x=np.eye(nobs),names=onames)
forecasts = "o2"
omitted = "p3"
e = ErrVar(jco=jco,parcov=parcov,obscov=obscov,forecasts=forecasts,
omitted_parameters=omitted)
svs = [0,1,2,3,4,5]
print(e.get_errvar_dataframe(svs))
def errvar_test_nonpest():
import numpy as np
from pyemu import ErrVar, Matrix, Cov
#non-pest
pnames = ["p1","p2","p3"]
onames = ["o1","o2","o3","o4"]
npar = len(pnames)
nobs = len(onames)
j_arr = np.random.random((nobs,npar))
jco = Matrix(x=j_arr,row_names=onames,col_names=pnames)
parcov = Cov(x=np.eye(npar),names=pnames)
obscov = Cov(x=np.eye(nobs),names=onames)
forecasts = "o2"
omitted = "p3"
e = ErrVar(jco=jco,parcov=parcov,obscov=obscov,forecasts=forecasts,
omitted_parameters=omitted)
svs = [0,1,2,3,4,5]
print(e.get_errvar_dataframe(svs))