scoreMDS: Scores MDS

scoreMDSR Documentation

Scores MDS

Description

Scores MDS

Usage

scoreMDS(
  df,
  id,
  label_scheme_sub,
  anal_type,
  scale_log2r,
  center_features,
  scale_features,
  dist_co = log2(1),
  adjEucDist = FALSE,
  choice = "cmdscale",
  method = "euclidean",
  p = 2L,
  k = 3L,
  col_group,
  folds,
  out_file,
  ...
)

Arguments

df

The name of a primary data file. By default, it will be determined automatically after matching the types of data and analysis with an id among c("pep_seq", "pep_seq_mod", "prot_acc", "gene"). A primary file contains normalized peptide or protein data and is among c("Peptide.txt", "Peptide_pVal.txt", "Peptide_impNA_pVal.txt", "Protein.txt", "Protein_pVal.txt", "protein_impNA_pVal.txt"). For analyses require the fields of significance p-values, the df will be one of c("Peptide_pVal.txt", "Peptide_impNA_pVal.txt", "Protein_pVal.txt", "protein_impNA_pVal.txt").

id

Character string; one of pep_seq, pep_seq_mod, prot_acc and gene.

label_scheme_sub

A data frame. Subset entries from label_scheme for selected samples.

anal_type

Character string; the type of analysis that are preset for method dispatch in function factories. The value will be determined automatically. Exemplary values include anal_type = c("PCA", "Corrplot", "EucDist", "GSPA", "Heatmap", "Histogram", "MDS", "Model", "NMF", "Purge", "Trend", "LDA", ...).

scale_log2r

Logical; if TRUE, adjusts log2FC to the same scale of standard deviation across all samples. The default is TRUE. At scale_log2r = NA, the raw log2FC without normalization will be used.

center_features

Logical; if TRUE, adjusts log2FC to center zero by features (proteins or peptides). The default is TRUE. Note the difference to data alignment with method_align in standPrn or standPep where log2FC are aligned by observations (samples).

scale_features

Logical; if TRUE, adjusts log2FC to the same scale of variance by features (protein or peptide entries). The default is TRUE. Note the difference to data scaling with scale_log2r where log2FC are scaled by observations (samples).

dist_co

Numeric; The cut-off in the absolute distance measured by d = abs(x_i - x_j). Data pairs, x_i and x_j, with corresponding d smaller than dist_co will be excluded from distance calculations by dist. The default is no distance cut-off at dist_co = log2(1).

adjEucDist

Logical; if TRUE, adjusts the inter-plex Euclidean distance by 1/sqrt(2) at method = "euclidean". The option adjEucDist = TRUE may be suitable when reference samples from each TMT plex undergo approximately the same sample handling process as the samples of interest. For instance, reference samples were split at the levels of protein lysates. Typically, adjEucDist = FALSE if reference samples were split near the end of a sample handling process, for instance, at the stages immediately before or after TMT labeling. Also see online README, section MDS for a brief reasoning.

choice

Character string; the MDS method in c("cmdscale", "isoMDS"). The default is "cmdscale".

method

Character string; the distance measure in one of c("euclidean", "maximum", "manhattan", "canberra", "binary") for dist. The default method is "euclidean".

p

Numeric; The power of the Minkowski distance in dist. The default is 2.

k

Numeric; The desired dimension for the solution passed to cmdscale. The default is 3.

col_group

Character string to a column key in expt_smry.xlsx. Samples corresponding to non-empty entries under col_group will be used for sample grouping in the indicated analysis. At the NULL default, the column key Group will be used. No data annotation by groups will be performed if the fields under the indicated group column is empty.

folds

Not currently used. Integer; the degree of folding data into subsets. The default is one without data folding.

out_file

A file path object to an output file.

...

filter_: Variable argument statements for the row filtration against data in a primary file linked to df. See also normPSM for the format of filter_ statements.

Additional parameters for ggsave:
width, the width of plot;
height, the height of plot
...


qzhang503/proteoQ documentation built on Dec. 14, 2024, 12:27 p.m.