RSSp_estimate: Estimate heritability using RSSp
In CreRecombinase/RSSp: Polygenic RSS

The three most important arguments

RSSp_estimate(
  quh,
  D,
  sample_size = NULL,
  trait_id = gen_trait_uuid(),
  pve_bounds = c(.Machine$double.eps, 1 - .Machine$double.eps),
  nterms = 1,
  p = sum(D),
  useGradient = T
)

`quh`	transformed summary statistics
`D`	vector of eigenvalues of the reference LD matrix
`sample_size`	sample size of GWAS
`trait_id`	identifier for the trait
`pve_bounds`	boundaries for pve estimation
`nterms`	number of terms in the model, by default there is one term that corresponds to heritability
`p`	number of variants in the dataset. By default (and in almost all scenarios) this is the sum of the eigenvalues.
`useGradient`	a boolean indicating whether or not to optimize using an analytic form for the gradient.