twoStageiLCA: Two-staged Independent Linked Component Analysis

In CHuanSite/SJD: Structured Joint Decomposition (SJD)

Two-staged Independent Linked Component Analysis

Description

Two-staged Independent Linked Component Analysis, a generalization based on the Two-staged Independent Linked Component Analysis

Usage

twoStageiLCA(
  dataset,
  group,
  comp_num,
  weighting = NULL,
  backup = 0,
  plotting = FALSE,
  proj_dataset = NULL,
  proj_group = NULL,
  enable_normalization = TRUE,
  column_sum_normalization = FALSE,
  screen_prob = NULL
)

Arguments

dataset	A list of dataset to be analyzed
group	A list of grouping of the datasets, indicating the relationship between datasets
comp_num	A vector indicates the dimension of each compoent
weighting	Weighting of each dataset, initialized to be NULL
backup	A positive scalar to determine how many ICs to over select
plotting	A boolean value to determine whether to plot the scree plot or not, default to be False
proj_dataset	The dataset(s) to be projected on.
proj_group	A listed of boolean combinations indicating which groupings should be used for each projected dataset.The length of proj_group should match the length of proj_dataset, and the length of each concatenated boolean combination should match the length of the parameter group.
enable_normalization	An argument to decide whether to use normalizaiton or not, default is TRUE
column_sum_normalization	An argument to decide whether to use column sum normalization or not, default it FALSE
screen_prob	A vector of probabilies for genes to be chosen

Value

A list contains the component and the score of each dataset on every component after 2siLCA algorithm

Examples

dataset = list(matrix(runif(5000, 1, 2), nrow = 100, ncol = 50),
matrix(runif(5000, 1, 2), nrow = 100, ncol = 50),
matrix(runif(5000, 1, 2), nrow = 100, ncol = 50),
matrix(runif(5000, 1, 2), nrow = 100, ncol = 50))
group = list(c(1,2,3,4), c(1,2), c(3,4), c(1,3), c(2,4), c(1), c(2), c(3), c(4))
comp_num = c(2,2,2,2,2,2,2,2,2)
proj_dataset = matrix(runif(5000, 1, 2), nrow = 100, ncol = 50)
proj_group = c(TRUE, TRUE, FALSE, FALSE, FALSE, TRUE, FALSE, FALSE, FALSE)
res_twoStageiLCA = twoStageiLCA(
dataset,
group,
comp_num,
proj_dataset = proj_dataset,
proj_group = proj_group)

CHuanSite/SJD documentation built on Nov. 29, 2024, 5:52 a.m.