code: def _kernel_rvs_beta1(x, bw):
      # Beta kernel for density, pdf, estimation
      # Kiriliouk, Segers, Tsukuhara 2020 arxiv, using bandwith 1/nobs sample
      return stats.beta.rvs(x / bw, (1 - x) / bw + 1)
  
file_path: statsmodels/statsmodels/distributions/copula/other_copulas.py
imports: import numpy
  from scipy import stats
  from statsmodels.tools.rng_qrng import check_random_state
  from statsmodels.distributions.copula.copulas import Copula
label_desc: Use image node_id 3 for calling a global function with example usage: _kernel_rvs_beta1(x, bw) and returns: stats