code: def emulate(self, X):
      """Emulates model output with new input data.
  
      Constructs orthonormal polynomials with new input matrix, X,
      and multiplies it with solution array, hdmr.x
  
      Compares emulated results with observed vector, Y.
  
      Returns
      ========
      Y_em : numpy.array
          Emulated output
      """
      # Calculate HDMR Basis Matrix
      b_m = _basis_matrix(X, self["hdmr"])
  
      # Get component functions
      Y_em = _comp_func(b_m, self["hdmr"], emulator=True)
  
      return np.sum(Y_em, axis=1)
  
imports: import math
  import warnings
  from typing import Dict, Tuple
  from types import MethodType
  from itertools import combinations, product
  from collections import defaultdict, namedtuple
  import numpy
  from pandas import DataFrame
  from numpy.linalg import det, pinv, matrix_rank
  from scipy.linalg import svd, LinAlgError, solve
  from scipy import stats, special
  from .None import common_args
  from ..util import read_param_file, ResultDict
label_desc: Use image node_id 18 for calling a global function with example usage: emulate(X) and returns: np