run cTWAS analysis using summary statistics with "no LD" version

ctwas_sumstats_noLD(
  z_snp,
  weights,
  region_info,
  snp_map,
  z_gene = NULL,
  thin = 0.1,
  niter_prefit = 3,
  niter = 30,
  init_group_prior = NULL,
  init_group_prior_var = NULL,
  maxSNP = Inf,
  min_nonSNP_PIP = 0.5,
  min_p_single_effect = 0.8,
  use_null_weight = TRUE,
  outputdir = NULL,
  outname = "ctwas_noLD",
  ncore = 1,
  logfile = NULL,
  verbose = FALSE,
  ...
)

Arguments

z_snp: A data frame with four columns: "id", "A1", "A2", "z". giving the z scores for snps. "A1" is effect allele. "A2" is the other allele.
weights: a list of pre-processed prediction weights.
region_info: a data frame of region definitions.
snp_map: a list of data frames with SNP-to-region map for the reference.
z_gene: A data frame with columns: "id", "z", giving the z-scores for genes.
thin: The proportion of SNPs to be used for estimating parameters and screening regions.
niter_prefit: the number of iterations of the E-M algorithm to perform during the initial parameter estimation step
niter: the number of iterations of the E-M algorithm to perform during the complete parameter estimation step
init_group_prior: a vector of initial values of prior inclusion probabilities for SNPs and genes.
init_group_prior_var: a vector of initial values of prior variances for SNPs and gene effects.
maxSNP: Inf or integer. Maximum number of SNPs in a region. Default is Inf, no limit. This can be useful if there are many SNPs in a region and you don't have enough memory to run the program.
min_nonSNP_PIP: Regions with non-SNP PIP >= min_nonSNP_PIP will be selected to run finemapping using full SNPs.
min_p_single_effect: Regions with probability >= min_p_single_effect of having at most one causal effect will be selected for the final EM step.
use_null_weight: If TRUE, allow for a probability of no effect in susie.
outputdir: The directory to store output. If specified, save outputs to the directory.
outname: The output name.
ncore: The number of cores used to parallelize susie over regions
logfile: The log filename. If NULL, print log info on screen.
verbose: If TRUE, print detailed messages
...: Additional arguments of susie_rss.

Value

a list, including z_gene, estimated parameters, region_data, cross-boundary genes, screening region results, and fine-mapping results.