code: def test_pdf_continuous(self):
      # Test for only continuous data
      dens = nparam.KDEMultivariate(
          data=[self.growth, self.Italy_gdp], var_type="cc", bw="cv_ls"
      )
      # take the first data points from the training set
      sm_result = np.squeeze(dens.pdf()[0:5])
      R_result = [1.6202284, 0.7914245, 1.6084174, 2.4987204, 1.3705258]
  
      ## CODE TO REPRODUCE THE RESULTS IN R
      ## library(np)
      ## data(oecdpanel)
      ## data (Italy)
      ## bw <-npudensbw(formula = ~oecdpanel$growth[1:50] + Italy$gdp[1:50],
      ## bwmethod ='cv.ls')
      ## fhat <- fitted(npudens(bws=bw))
      ## fhat[1:5]
      npt.assert_allclose(sm_result, R_result, atol=1e-3)
  
imports: import numpy
  import numpy.testing
  import pytest
  from numpy.testing import assert_allclose, assert_equal
  import statsmodels.api
label_desc: Use image node_id 4 for calling the TestKDEMultivariate obj's underlying member method code with example usage: obj.test_pdf_continuous() without return types