Integrative sparse canonical correlation analysis

Description

This function provides a penalty-based integrative sparse canonical correlation analysis method to handle the multiple datasets with high dimensions generated under similar protocols, which consists of two built-in penalty items for selecting the important variables for users to choose, and two contrasted penalty functions for eliminating the diffierence (magnitude or sign) between estimators within each group.

Usage

iscca(x, y, L, mu1, mu2, mu3, mu4, eps = 1e-04, pen1 = "homogeneity",
  pen2 = "magnitude", scale.x = TRUE, scale.y = TRUE, maxstep = 50,
  submaxstep = 10, trace = FALSE, draw = FALSE)

Arguments

Value

An 'iscca' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• loading.x: the estimated canonical vector of variables x.

• loading.y: the estimated canonical vector of variables y.

• variable.x: the screening results of variables x.

• variable.y: the screening results of variables y.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

See Also

See Also as preview.cca, iscca.cv, meta.scca, scca.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.cca")
x <- simData.cca$x
y <- simData.cca$y
L <- length(x)
mu1 <- mu3 <- 0.4
mu2 <- mu4 <- 2.5

prev_cca <- preview.cca(x = x, y = y, L = L, scale.x = TRUE, scale.y = TRUE)
res_homo_m <- iscca(x = x, y = y, L = L, mu1 = mu1, mu2 = mu2, mu3 = mu3, mu4 = mu4,
                    eps = 5e-2, maxstep = 50, submaxstep = 10, trace = TRUE, draw = TRUE)


res_homo_s <- iscca(x = x, y = y, L = L, mu1 = mu1, mu2 = mu2, mu3 = mu3, mu4 = mu4,
                    eps = 5e-2, pen1 = "homogeneity", pen2 = "sign", scale.x = TRUE,
                    scale.y = TRUE, maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

mu1 <- mu3 <- 0.3
mu2 <- mu4 <- 2
res_hete_m <- iscca(x = x, y = y, L = L, mu1 = mu1, mu2 = mu2, mu3 = mu3, mu4 = mu4,
                    eps = 5e-2, pen1 = "heterogeneity", pen2 = "magnitude", scale.x = TRUE,
                    scale.y = TRUE, maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

res_hete_s <- iscca(x = x, y = y, L = L, mu1 = mu1, mu2 = mu2, mu3 = mu3, mu4 = mu4,
                    eps = 5e-2, pen1 = "heterogeneity", pen2 = "sign", scale.x = TRUE,
                    scale.y = TRUE, maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

Cross-validation for iscca

Description

Performs K-fold cross validation for the integrative sparse canonical correlation analysis over a grid of values for the regularization parameter mu1, mu2, mu3 and mu4.

Usage

iscca.cv(x, y, L, K = 5, mu1, mu2, mu3, mu4, eps = 1e-04,
  pen1 = "homogeneity", pen2 = "magnitude", scale.x = TRUE,
  scale.y = TRUE, maxstep = 50, submaxstep = 10)

Arguments

Value

An 'iscca.cv' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• mu1: the sparsity penalty parameter selected from the feasible set of parameter mu1 provided by users.

• mu2: the contrasted penalty parameter selected from the feasible set of parameter mu2 provided by users.

• mu3: the sparsity penalty parameter selected from the feasible set of parameter mu3 provided by users.

• mu4: the contrasted penalty parameter selected from the feasible set of parameter mu4 provided by users.

• fold: The fold assignments for cross-validation for each observation.

• loading.x: the estimated canonical vector of variables x with selected tuning parameters.

• loading.y: the estimated canonical vector of variables y with selected tuning parameters.

• variable.x: the screening results of variables x.

• variable.y: the screening results of variables y.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

See Also

See Also as iscca.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.cca")
x <- simData.cca$x
y <- simData.cca$y
L <- length(x)
mu1 <- c(0.2, 0.4)
mu3 <- 0.4
mu2 <- mu4 <- 2.5

res_homo_m <- iscca.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, mu3 = mu3,
                       mu4 = mu4, eps = 1e-2, pen1 = "homogeneity", pen2 = "magnitude",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

res_homo_s <- iscca.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, mu3 = mu3,
                       mu4 = mu4, eps = 1e-2, pen1 = "homogeneity", pen2 = "sign",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

mu1 <- mu3 <- c(0.1, 0.3)
mu2 <- mu4 <- 2
res_hete_m <- iscca.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, mu3 = mu3,
                       mu4 = mu4, eps = 1e-2, pen1 = "heterogeneity", pen2 = "magnitude",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

res_hete_s <- iscca.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, mu3 = mu3,
                       mu4 = mu4, eps = 1e-2, pen1 = "heterogeneity", pen2 = "sign",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

Plot the results of iscca

Description

Plot the convergence path graph in the integrative sparse canonical correlation analysis method or show the the first pair of canonical vectors.

Usage

iscca.plot(x, type)

Arguments

Details

See details in iscca.

Value

the convergence path graph or the scatter diagrams of the first pair of canonical vectors.

Examples

library(iSFun)
data("simData.cca")
x <- simData.cca$x
y <- simData.cca$y
L <- length(x)
mu1 <- mu3 <- 0.4
mu2 <- mu4 <- 2.5

res_homo_m <- iscca(x = x, y = y, L = L, mu1 = mu1, mu2 = mu2, mu3 = mu3,
                    mu4 = mu4, eps = 5e-2, maxstep = 100, trace = FALSE, draw = FALSE)
iscca.plot(x = res_homo_m, type = "path")
iscca.plot(x = res_homo_m, type = "loading")

Integrative sparse principal component analysis

Description

This function provides a penalty-based integrative sparse principal component analysis method to obtain the direction of first principal component of the multiple datasets with high dimensions generated under similar protocols, which consists of two built-in penalty items for selecting the important variables for users to choose, and two contrasted penalty functions for eliminating the diffierence (magnitude or sign) between estimators within each group.

Usage

ispca(x, L, mu1, mu2, eps = 1e-04, pen1 = "homogeneity",
  pen2 = "magnitude", scale.x = TRUE, maxstep = 50,
  submaxstep = 10, trace = FALSE, draw = FALSE)

Arguments

Value

An 'ispca' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• eigenvalue: the estimated first eigenvalue.

• eigenvector: the estimated first eigenvector.

• component: the estimated first component.

• variable: the screening results of variables.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

References

• Fang K, Fan X, Zhang Q, et al. Integrative sparse principal component analysis[J]. Journal of Multivariate Analysis, 2018, 166: 1-16.

See Also

See Also as preview.pca, ispca.cv, meta.spca, spca.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.pca")
x <- simData.pca$x
L <- length(x)

prev_pca <- preview.pca(x = x, L = L, scale.x = TRUE)
res_homo_m <- ispca(x = x, L = L, mu1 = 0.5, mu2 = 0.002, trace = TRUE, draw = TRUE)


res_homo_s <- ispca(x = x, L = L, mu1 = 0.5, mu2 = 0.002,
                    pen1 = "homogeneity", pen2 = "sign", scale.x = TRUE,
                    maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

res_hete_m <- ispca(x = x, L = L, mu1 = 0.1, mu2 = 0.05,
                    pen1 = "heterogeneity", pen2 = "magnitude", scale.x = TRUE,
                    maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

res_hete_s <- ispca(x = x, L = L, mu1 = 0.1, mu2 = 0.05,
                    pen1 = "heterogeneity", pen2 = "sign", scale.x = TRUE,
                    maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

Cross-validation for ispca

Description

Performs K-fold cross validation for the integrative sparse principal component analysis over a grid of values for the regularization parameter mu1 and mu2.

Usage

ispca.cv(x, L, K = 5, mu1, mu2, eps = 1e-04, pen1 = "homogeneity",
  pen2 = "magnitude", scale.x = TRUE, maxstep = 50,
  submaxstep = 10)

Arguments

Value

An 'ispca.cv' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• mu1: the sparsity penalty parameter selected from the feasible set of parameter mu1 provided by users.

• mu2: the contrasted penalty parameter selected from the feasible set of parameter mu2 provided by users.

• fold: The fold assignments for cross-validation for each observation.

• eigenvalue: the estimated first eigenvalue with selected tuning parameters mu1 and mu2.

• eigenvector: the estimated first eigenvector with selected tuning parameters mu1 and mu2.

• component: the estimated first component with selected tuning parameters mu1 and mu2.

• variable: the screening results of variables.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

References

• Fang K, Fan X, Zhang Q, et al. Integrative sparse principal component analysis[J]. Journal of Multivariate Analysis, 2018, 166: 1-16.

See Also

See Also as ispca.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.pca")
x <- simData.pca$x
L <- length(x)
mu1 <- c(0.3, 0.5)
mu2 <- 0.002

res_homo_m <- ispca.cv(x = x, L = L, K = 5, mu1 = mu1, mu2 = mu2, pen1 = "homogeneity",
                       pen2 = "magnitude", scale.x = TRUE, maxstep = 50, submaxstep = 10)

res_homo_s <- ispca.cv(x = x, L = L, K = 5, mu1 = mu1, mu2 = mu2, pen1 = "homogeneity",
                       pen2 = "sign", scale.x = TRUE, maxstep = 50, submaxstep = 10)

mu1 <- c(0.1, 0.15)
mu2 <- 0.05
res_hete_m <- ispca.cv(x = x, L = L, K = 5, mu1 = mu1, mu2 = mu2, pen1 = "heterogeneity",
                       pen2 = "magnitude", scale.x = TRUE, maxstep = 50, submaxstep = 10)

res_hete_s <- ispca.cv(x = x, L = L, K = 5, mu1 = mu1, mu2 = mu2, pen1 = "heterogeneity",
                       pen2 = "sign", scale.x = TRUE, maxstep = 50, submaxstep = 10)

Plot the results of ispca

Description

Plot the convergence path graph or estimated value of the first eigenvector u in the integrative sparse principal component analysis method.

Usage

ispca.plot(x, type)

Arguments

Details

See details in ispca.

Value

the convergence path graph or the scatter diagrams of the first eigenvector u.

Examples

library(iSFun)
data("simData.pca")
x <- simData.pca$x
L <- length(x)

res_homo_m <- ispca(x = x, L = L, mu1 = 0.5, mu2 = 0.002, trace = FALSE, draw = FALSE)
ispca.plot(x = res_homo_m, type = "path")
ispca.plot(x = res_homo_m, type = "loading")

Integrative sparse partial least squares

Description

This function provides a penalty-based integrative sparse partial least squares method to handle the multiple datasets with high dimensions generated under similar protocols, which consists of two built-in penalty items for selecting the important variables for users to choose, and two contrasted penalty functions for eliminating the diffierence (magnitude or sign) between estimators within each group.

Usage

ispls(x, y, L, mu1, mu2, eps = 1e-04, kappa = 0.05,
  pen1 = "homogeneity", pen2 = "magnitude", scale.x = TRUE,
  scale.y = TRUE, maxstep = 50, submaxstep = 10, trace = FALSE,
  draw = FALSE)

Arguments

Value

An 'ispls' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• betahat: the estimated regression coefficients.

• loading: the estimated first direction vector.

• variable: the screening results of variables x.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

References

• Liang W, Ma S, Zhang Q, et al. Integrative sparse partial least squares[J]. Statistics in Medicine, 2021, 40(9): 2239-2256.

See Also

See Also as preview.pls, ispls.cv, meta.spls, spls.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.pls")
x <- simData.pls$x
y <- simData.pls$y
L <- length(x)

prev_pls <- preview.pls(x, y, L, scale.x = TRUE, scale.y = TRUE)
res_homo_m <- ispls(x = x, y = y, L = L, mu1 = 0.05, mu2 = 0.25,
                    eps = 5e-2, trace = TRUE, draw = TRUE)


res_homo_s <- ispls(x = x, y = y, L = L, mu1 = 0.05, mu2 = 0.25,
                    eps = 5e-2, kappa = 0.05, pen1 = "homogeneity",
                    pen2 = "sign", scale.x = TRUE, scale.y = TRUE,
                    maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

res_hete_m <- ispls(x = x, y = y, L = L, mu1 = 0.05, mu2 = 0.25,
                    eps = 5e-2, kappa = 0.05, pen1 = "heterogeneity",
                    pen2 = "magnitude", scale.x = TRUE, scale.y = TRUE,
                    maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

res_hete_s <- ispls(x = x, y = y, L = L, mu1 = 0.05, mu2 = 0.25,
                    eps = 5e-2, kappa = 0.05, pen1 = "heterogeneity",
                    pen2 = "sign", scale.x = TRUE, scale.y = TRUE,
                    maxstep = 50, submaxstep = 10, trace = FALSE, draw = FALSE)

Cross-validation for ispls

Description

Performs K-fold cross validation for the integrative sparse partial least squares over a grid of values for the regularization parameter mu1 and mu2.

Usage

ispls.cv(x, y, L, K, mu1, mu2, eps = 1e-04, kappa = 0.05,
  pen1 = "homogeneity", pen2 = "magnitude", scale.x = TRUE,
  scale.y = TRUE, maxstep = 50, submaxstep = 10)

Arguments

Value

An 'ispls.cv' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• mu1: the sparsity penalty parameter selected from the feasible set of parameter mu1 provided by users.

• mu2: the contrasted penalty parameter selected from the feasible set of parameter mu2 provided by users.

• fold: The fold assignments for cross-validation for each observation.

• betahat: the estimated regression coefficients with selected tuning parameters mu1 and mu2.

• loading: the estimated first direction vector with selected tuning parameters mu1 and mu2.

• variable: the screening results of variables x.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

References

• Liang W, Ma S, Zhang Q, et al. Integrative sparse partial least squares[J]. Statistics in Medicine, 2021, 40(9): 2239-2256.

See Also

See Also as ispls.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.pls")
x <- simData.pls$x
y <- simData.pls$y
L <- length(x)
mu1 <- c(0.04, 0.05)
mu2 <- 0.25

res_homo_m <- ispls.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, eps = 1e-2,
                       kappa = 0.05, pen1 = "homogeneity", pen2 = "magnitude",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

res_homo_s <- ispls.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, eps = 1e-2,
                       kappa = 0.05, pen1 = "homogeneity", pen2 = "sign",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

res_hete_m <- ispls.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, eps = 1e-2,
                       kappa = 0.05, pen1 = "heterogeneity", pen2 = "magnitude",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

res_hete_s <- ispls.cv(x = x, y = y, L = L, K = 5, mu1 = mu1, mu2 = mu2, eps = 1e-2,
                       kappa = 0.05, pen1 = "heterogeneity", pen2 = "sign",
                       scale.x = TRUE, scale.y = TRUE, maxstep = 50, submaxstep = 10)

Plot the results of ispls

Description

Plot the convergence path graph of the first direction vector w in the integrative sparse partial least squares model or show the regression coefficients.

Usage

ispls.plot(x, type)

Arguments

Details

See details in ispls.

Value

show the convergence path graph of the first direction vector w or the regression coefficients.

Examples

library(iSFun)
data("simData.pls")
x <- simData.pls$x
y <- simData.pls$y
L <- length(x)

res_homo_m <- ispls(x = x, y = y, L = L, mu1 = 0.05, mu2 = 0.25,
                    eps = 5e-2, trace = FALSE, draw = FALSE)
ispls.plot(x = res_homo_m, type = "path")
ispls.plot(x = res_homo_m, type = "loading")
ispls.plot(x = res_homo_m, type = "heatmap")

Meta-analytic sparse canonical correlation analysis method in integrative study

Description

This function provides penalty-based sparse canonical correlation meta-analytic method to handle the multiple datasets with high dimensions generated under similar protocols, which is based on the principle of maximizing the summary statistics S.

Usage

meta.scca(x, y, L, mu1, mu2, eps = 1e-04, scale.x = TRUE,
  scale.y = TRUE, maxstep = 50, trace = FALSE)

Arguments

Value

A 'meta.scca' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• loading.x: the estimated canonical vector of variables x.

• loading.y: the estimated canonical vector of variables y.

• variable.x: the screening results of variables x.

• variable.y: the screening results of variables y.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

References

• Cichonska A, Rousu J, Marttinen P, et al. metaCCA: summary statistics-based multivariate meta-analysis of genome-wide association studies using canonical correlation analysis[J]. Bioinformatics, 2016, 32(13): 1981-1989.

See Also

See Also as iscca, scca.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.cca")
x <- simData.cca$x
y <- simData.cca$y
L <- length(x)
mu1 <- 0.08
mu2 <- 0.08

res <- meta.scca(x = x, y = y, L = L, mu1 = mu1, mu2 = mu2, trace = TRUE)

Meta-analytic sparse principal component analysis method in integrative study

Description

This function provides penalty-based sparse principal component meta-analytic method to handle the multiple datasets with high dimensions generated under similar protocols, which is based on the principle of maximizing the summary statistics S.

Usage

meta.spca(x, L, mu1, eps = 1e-04, scale.x = TRUE, maxstep = 50,
  trace = FALSE)

Arguments

Value

A 'meta.spca' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• eigenvalue: the estimated first eigenvalue.

• eigenvector: the estimated first eigenvector.

• component: the estimated first component.

• variable: the screening results of variables.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

References

• Kim S H, Kang D, Huo Z, et al. Meta-analytic principal component analysis in integrative omics application[J]. Bioinformatics, 2018, 34(8): 1321-1328.

See Also

See Also as ispca, spca.

Examples

library(iSFun)
data("simData.pca")
x <- simData.pca$x
L <- length(x)

res <- meta.spca(x = x, L = L, mu1 = 0.5, trace = TRUE)

Meta-analytic sparse partial least squares method in integrative study

Description

This function provides penalty-based sparse canonical correlation meta-analytic method to handle the multiple datasets with high dimensions generated under similar protocols, which is based on the principle of maximizing the summary statistics.

Usage

meta.spls(x, y, L, mu1, eps = 1e-04, kappa = 0.05, scale.x = TRUE,
  scale.y = TRUE, maxstep = 50, trace = FALSE)

Arguments

Value

A 'meta.spls' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• betahat: the estimated regression coefficients.

• loading: the estimated first direction vector.

• variable: the screening results of variables x.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

See Also

See Also as ispls, spls.

Examples

library(iSFun)
data("simData.pls")
x <- simData.pls$x
y <- simData.pls$y
L <- length(x)

res <- meta.spls(x = x, y = y, L = L, mu1 = 0.03, trace = TRUE)

Statistical description before using function iscca

Description

The function describes the basic statistical information of the data, including sample mean, sample variance of X and Y, and the first pair of canonical vectors.

Usage

preview.cca(x, y, L, scale.x = TRUE, scale.y = TRUE)

Arguments

Value

An 'preview.cca' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• loading.x: the estimated canonical vector of variables x.

• loading.y: the estimated canonical vector of variables y.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

See Also

See Also as iscca.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.cca")
x <- simData.cca$x
y <- simData.cca$y
L <- length(x)

prev_cca <- preview.cca(x = x, y = y, L = L, scale.x = TRUE, scale.y = TRUE)

Statistical description before using function ispca

Description

The function describes the basic statistical information of the data, including sample mean, sample co-variance of X and Y, the first eigenvector, eigenvalue and principal component, etc.

Usage

preview.pca(x, L, scale.x = TRUE)

Arguments

Value

An 'preview.pca' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• eigenvalue: the estimated first eigenvalue.

• eigenvector: the estimated first eigenvector.

• component: the estimated first component.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

See Also

See Also as ispca.

Examples

# Load a list with 3 data sets
library(iSFun)
data("simData.pca")
x <- simData.pca$x
L <- length(x)

prev.pca <- preview.pca(x = x, L = L, scale.x = TRUE)

Statistical description before using function ispls

Description

The function describes the basic statistical information of the data, including sample mean, sample variance of X and Y, the first direction of partial least squares method, etc.

Usage

preview.pls(x, y, L, scale.x = TRUE, scale.y = TRUE)

Arguments

Value

A 'preview.pls' object that contains the list of the following items.

• x: list of data matrices, L datasets of explanatory variables with centered columns. If scale.x is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• y: list of data matrices, L datasets of dependent variables with centered columns. If scale.y is TRUE, the columns of L datasets are standardized to have mean 0 and standard deviation 1.

• loading: the estimated first direction vector.

• meanx: list of numeric vectors, column mean of the original datasets x.

• normx: list of numeric vectors, column standard deviation of the original datasets x.

• meany: list of numeric vectors, column mean of the original datasets y.

• normy: list of numeric vectors, column standard deviation of the original datasets y.

See Also

See Also as ispls.

Examples

library(iSFun)
data("simData.pls")
x <- simData.pls$x
y <- simData.pls$y
L <- length(x)

prev_pls <- preview.pls(x = x, y = y, L = L, scale.x = TRUE, scale.y = TRUE)

Sparse canonical correlation analysis

Description

This function provides penalty-based sparse canonical correlation analysis to get the first pair of canonical vectors.

Usage

scca(x, y, mu1, mu2, eps = 1e-04, scale.x = TRUE, scale.y = TRUE,
  maxstep = 50, trace = FALSE)

Arguments

Value

An 'scca' object that contains the list of the following items.

• x: data matrix of explanatory variables with centered columns. If scale.x is TRUE, the columns of data matrix are standardized to have mean 0 and standard deviation 1.

• y: data matrix of dependent variables with centered columns. If scale.y is TRUE, the columns of data matrix are standardized to have mean 0 and standard deviation 1.

• loading.x: the estimated canonical vector of variables x.

• loading.y: the estimated canonical vector of variables y.

• variable.x: the screening results of variables x.

• variable.y: the screening results of variables y.

• meanx: column mean of the original dataset x.

• normx: column standard deviation of the original dataset x.

• meany: column mean of the original dataset y.

• normy: column standard deviation of the original dataset y.

See Also

See Also as iscca, meta.scca.

Examples

library(iSFun)
data("simData.cca")
x.scca <- do.call(rbind, simData.cca$x)
y.scca <- do.call(rbind, simData.cca$y)
res_scca <- scca(x = x.scca, y = y.scca, mu1 = 0.1, mu2 = 0.1, eps = 1e-3,
                 scale.x = TRUE, scale.y = TRUE, maxstep = 50, trace = FALSE)

Example data for method iscca

Description

Example data for users to apply the method iscca, iscca.cv, meta.scca or scca.

Format

list

Example data for method ispca

Description

Example data for users to apply the method ispca, ispca.cv, meta.spca or spca.

Format

list

Example data for method ispls

Description

Example data for users to apply the method ispls, ispls.cv, meta.spls or spls.

Format

list

Sparse principal component analysis

Description

This function provides penalty-based integrative sparse principal component analysis to obtain the direction of first principal component of a given dataset with high dimensions.

Usage

spca(x, mu1, eps = 1e-04, scale.x = TRUE, maxstep = 50,
  trace = FALSE)

Arguments

Value

An 'spca' object that contains the list of the following items.

• x: data matrix of explanatory variables with centered columns. If scale.x is TRUE, the columns of data matrix are standardized to have mean 0 and standard deviation 1.

• eigenvalue: the estimated first eigenvalue.

• eigenvector: the estimated first eigenvector.

• component: the estimated first principal component.

• variable: the screening results of variables.

• meanx: column mean of the original dataset x.

• normx: column standard deviation of the original dataset x.

See Also

See Also as ispca, meta.spca.

Examples

library(iSFun)
data("simData.pca")
x.spca <- do.call(rbind, simData.pca$x)
res_spca <- spca(x = x.spca, mu1 = 0.08, eps = 1e-3, scale.x = TRUE,
                 maxstep = 50, trace = FALSE)

Sparse partial least squares

Description

This function provides penalty-based sparse partial least squares analysis for single dataset with high dimensions., which aims to have the direction of the first loading.

Usage

spls(x, y, mu1, eps = 1e-04, kappa = 0.05, scale.x = TRUE,
  scale.y = TRUE, maxstep = 50, trace = FALSE)

Arguments

Value

An 'spls' object that contains the list of the following items.

• x: data matrix of explanatory variables with centered columns. If scale.x is TRUE, the columns of data matrix are standardized to have mean 0 and standard deviation 1.

• y: data matrix of dependent variables with centered columns. If scale.y is TRUE, the columns of data matrix are standardized to have mean 0 and standard deviation 1.

• betahat: the estimated regression coefficients.

• loading: the estimated first direction vector.

• variable: the screening results of variables.

• meanx: column mean of the original dataset x.

• normx: column standard deviation of the original dataset x.

• meany: column mean of the original dataset y.

• normy: column standard deviation of the original dataset y.

See Also

See Also as ispls, meta.spls.

Examples

library(iSFun)
data("simData.pls")
x.spls <- do.call(rbind, simData.pls$x)
y.spls <- do.call(rbind, simData.pls$y)
res_spls <- spls(x = x.spls, y = y.spls, mu1 = 0.05, eps = 1e-3, kappa = 0.05,
                 scale.x = TRUE, scale.y = TRUE, maxstep = 50, trace = FALSE)