code: def test_bernoulli_sampling_thin(posterior):
      fit = posterior.sample(num_thin=2)
      assert fit["theta"].shape[-1] == 500 * 4
  
file_path: pystan/tests/test_basic_bernoulli.py
imports: import numpy
  import pytest
  import stan
label_desc: Use image node_id 3 for calling a global function with example usage: test_bernoulli_sampling_thin(posterior) without return types