maxes_1seg: Get peak heights and permuted peak heights for one segment

Description Usage Arguments Value

Description

Finds peak heights and permuatation peak heights for one segment

Usage

1
2
3
maxes_1seg(dat.file, pheno.file, s0, zmin, seed, n.perm,
  save.perm.stats = FALSE, segment.bounds = NULL, out.file = NULL,
  max1.file = NULL, z0 = zmin * 0.3, bandwidth = 50, maxit = 50)

Arguments

dat.file

Data file. Should have first column as position and a header giving sample names. One row per base pair.

pheno.file

Phenotype file. First column should correspond to the header of dat.file.

s0

Variance inflation constant.

zmin

Lower bound on significance thresholds.

seed

Set a seed for running permutations (Required)

n.perm

Number of permutations

save.perm.stats

Do you want to save all the permutation statistics?

segment.bounds

Length 2 vector giving start and stop base-pair for the segment. If NULL, all base-pairs in dat.file are used.

out.file

Output file (optional). If NULL will append "_mx.RData" to the leading portion of the dat.file name.

max1.file

File with data from running with n.perm=0

z0

Reference level for merging

bandwidth

Smoother bandwidth.

maxit

Maximum iterations for rlm.

Value

Saves an RData file with maxes list and other info. Returns the length of that list.


jean997/fret documentation built on May 18, 2019, 11:43 p.m.