,
Guanqun Meng [aut],
Wen Tang [aut]| Title: | Cell-Type-Specific Power Assessment |
|---|---|
| Description: | CYPRESS is a cell-type-specific power tool. This package aims to perform power analysis for the cell-type-specific data. It calculates FDR, FDC, and power, under various study design parameters, including but not limited to sample size, and effect size. It takes the input of a SummarizeExperimental(SE) object with observed mixture data (feature by sample matrix), and the cell-type mixture proportions (sample by cell-type matrix). It can solve the cell-type mixture proportions from the reference free panel from TOAST and conduct tests to identify cell-type-specific differential expression (csDE) genes. |
| Authors: | Shilin Yu [aut, cre] ,
Guanqun Meng [aut],
Wen Tang [aut] |
| Maintainer: | Shilin Yu <[email protected]> |
| License: | GPL-2 | GPL-3 |
| Version: | 1.1.1 |
| Built: | 2024-07-27 03:12:28 UTC |
| Source: | https://github.com/bioc/cypress |
The cypress package is specifically designed to perform comprehensive cell-type-specific power assessment for differential expression using RNA-sequencing experiments. It accepts real Bulk RNAseq data as input for parameter estimation and simulation. The tool provides flexibility by allowing users to customize sample sizes, number of cell types, and effect sizes (log-fold change). Additionally, it computes statistical power, true discovery rate (TDR), and false discovery cost (FDC) under different scenarios as results.
cypress is the first statistical tool to evaluate the power in cell-type-specific Differentially Expressed (csDE) genes detection experiments from a prospective way by letting researchers be flexible in tuning sample sizes, effect sizes, csDE genes percentage, total number of genetic features, type I error control, etc.
cypress offers 3 options for simulation and power evaluation: simFromData() ,simFromParam() and quickPower(). If users have their own bulk RNA-seq count data, they can use the simFromData() function; otherwise, they can use simFromParam(), which uses one of the three sets of simulation parameters estimated from existing studies, to perform power evaluation under user-defined simulation settings. If users prefer to quickly examine the power evaluation results and do not want to run the simulation, they can use the quickPower() function to view results from our existing simulations. The output of these 3 functions is a S4 object including a list of simulation results under various experimental settings, including as statistical power, TDR, and FDC.
Once users have obtained an S4 object with a list of results output from either simFromData(), simFromParam() or quickPower(), they can use the following functions to generate basic evaluation plots: plotPower(), plotTDR() and plotFDC().
Shilin Yu <[email protected]> Guanqun Meng <[email protected]> Wen Tang <[email protected]>
An S4 object that stores parameter estimates associated with the Autism Spectrum Disorder (ASD) dataset. This object contains a variety of numerical vectors and matrices representing different statistical parameters used in the simulation.
data('quickParaASD')data('quickParaASD')
Simulation parameters for simFromParam function.
health_lm_meanA numeric vector containing the log-mean parameter estimates for each cell type from healthy samples.
health_lm_mean_dA matrix containing the variance-covariance estimates of log-mean parameters across cell types from healthy samples.
lod_mA numeric vector containing the log-dispersion parameter estimates for each cell type from healthy samples.
lod_dA matrix containing the variance-covariance estimates of log-dispersion parameters across cell types from healthy samples.
health_alphaA numeric vector of the estimated alpha parameter used to simulate cell type proportions for healthy samples.
case_alphaA numeric vector of the estimated alpha parameter used to simulate cell type proportions for case samples.
One S4 object.
data('quickParaASD')data('quickParaASD')
Pre-calculated power evaluation results from Autism Spectrum Disorder (ASD) study. The results
can be used to create plots using plot functions (plotFDC, plotPower, plotTDR).
data('quickPowerASD')data('quickPowerASD')
A S4 object.
ct_TDR_bio_smryCell-type-specific target TDR.
TDR_bio_smryAverage target TDR.
ct_PWR_bio_smryCell-type-specific target power.
PWR_bio_smryAverage target power.
PWR_strata_ct_bio_smryCell type specific target power by gene expression stratification.
PWR_strata_bio_smryAverage target power by gene expression stratification.
ct_FDC_bio_smryCell type specific target FDC.
FDC_bio_smryAverage target FDC.
One S4 object.
data('quickPowerASD')data('quickPowerASD')
ASD_prop is an example of SummarizedExperiment (SE) object input for the simFromData function. It contains the following elements:
countsA gene expression value dataset from Autism Spectrum Disorder (ASD) study, in the form of raw read counts, 29674 genes by 48 samples, with 24 cases and 24 controls
colDataSample meta-data. The first column is the group status (i.e. case/ctrl), the second column is the subject ID. The remaining are the cell type proportions of all samples.
data(ASD_prop_se)data(ASD_prop_se)
SE object.
One SE object.
data(ASD_prop_se)data(ASD_prop_se)
The cypress_out and est_out classes are custom S4 classes in the cypress package. both classes are designed as a comprehensive container for various types of analysis results.
Class for cypress.
The cypress_out class is a S4 class in the cypress package. This class is customized to better present results and use for cypress plot functions.
ct_TDR_bio_smryCell type specific target TDR
TDR_bio_smryTarget TDR
ct_PWR_bio_smryCell type specific target power
PWR_bio_smryTarget power
PWR_strata_bio_smryTarget power by gene expression stratification
PWR_strata_ct_bio_smryCell type specific target power by gene expression stratification
ct_FDC_bio_smryCell type specific target FDC
FDC_bio_smryTarget FDC.
The est_out class is designed to output the parameter estimated results, providing a structured representation of results.
health_alphaControl group proportion simulation parameter.
case_alphaCase group proportion simulation parameter.
health_lmean_mMean of genetic distribution mean for each cell.
health_lmean_dVar/cov of genetic distribution mean among cell types
lod_mMean of genetic distribution dispersion for each cell.
lod_dVar/cov of genetic distribution dispersion among cell types.
sample_CT_propMatrix of sample Cell type proportions.
genenameGene Name.
samplenameSample Name.
CTnameCell type names.
dimensions_Z_hat_arydimensions for the Z hat array.
Shilin Yu <[email protected]>
data(quickParaGSE60424)data(quickParaGSE60424)
Accessor function for getting or replace slots. Show methods for cypress object.
getcypress(object, name) setcypress(object, name, value)getcypress(object, name) setcypress(object, name, value)
object |
object from cypress. |
name |
name of the slot in cypress object. |
value |
value of the slot in cypress object. |
Methods for cypress.
data(quickPowerIBD) getcypress(ibd_propPower, "ct_TDR_bio_smry")data(quickPowerIBD) getcypress(ibd_propPower, "ct_TDR_bio_smry")
An S4 object that stores simulation parameter estimated from the immune-related disease (IAD) study (GSE60424). This object contains a variety of numerical vectors and matrices representing different statistical parameters used in the simulation. The patients were drawn from healthy subjects from the immune-related diseases study.
data('quickParaGSE60424')data('quickParaGSE60424')
Simulation parameters for simFromParam function.
health_lm_meanA numeric vector containing the log-mean parameter estimates for each cell type from healthy samples.
health_lm_mean_dA matrix containing the variance-covariance estimates of log-mean parameters across cell types from healthy samples.
lod_mA numeric vector containing the log-dispersion parameter estimates for each cell type from healthy samples.
lod_dA matrix containing the variance-covariance estimates of log-dispersion parameters across cell types from healthy samples.
health_alphaA numeric vector of the estimated alpha parameter used to simulate cell type proportions for healthy samples.
case_alphaA numeric vector of the estimated alpha parameter used to simulate cell type proportions for case samples.
One S4 object.
data('quickParaGSE60424')data('quickParaGSE60424')
Pre-calculated power evaluation results from immune-related disease (IAD) study (GSE60424). The results
can be used to create plots using plot functions (plotFDC, plotPower, plotTDR).
data('quickPowerGSE60424')data('quickPowerGSE60424')
A S4 object.
ct_TDR_bio_smryCell-type-specific target TDR.
TDR_bio_smryAverage target TDR.
ct_PWR_bio_smryCell-type-specific target power.
PWR_bio_smryAverage target power.
PWR_strata_ct_bio_smryCell type specific target power by gene expression stratification.
PWR_strata_bio_smryAverage target power by gene expression stratification.
ct_FDC_bio_smryCell type specific target FDC.
FDC_bio_smryAverage target FDC.
One S4 object.
data('quickPowerGSE60424')data('quickPowerGSE60424')
An S4 object that stores simulation parameter estimated from the inflammatory bowel disease (IBD) study (GSE57945). This object contains a variety of numerical vectors and matrices representing different statistical parameters used in the simulation. The patients were drawn from healthy subjects from the immune-related diseases study.
data('quickParaIBD')data('quickParaIBD')
Simulation parameters for simFromParam function.
health_lm_meanA numeric vector containing the log-mean parameter estimates for each cell type from healthy samples.
health_lm_mean_dA matrix containing the variance-covariance estimates of log-mean parameters across cell types from healthy samples.
lod_mA numeric vector containing the log-dispersion parameter estimates for each cell type from healthy samples.
lod_dA matrix containing the variance-covariance estimates of log-dispersion parameters across cell types from healthy samples.
health_alphaA numeric vector of the estimated alpha parameter used to simulate cell type proportions for healthy samples.
case_alphaA numeric vector of the estimated alpha parameter used to simulate cell type proportions for case samples.
One S4 object.
data('quickParaIBD')data('quickParaIBD')
Pre-calculated power evaluation results from pediatric inflammatory bowel disease (IBD) study (GSE57945). The results
can be used to create plots using plot functions (plotFDC, plotPower, plotTDR).
data('quickPowerIBD')data('quickPowerIBD')
A S4 object.
ct_TDR_bio_smryCell-type-specific target TDR.
TDR_bio_smryAverage target TDR.
ct_PWR_bio_smryCell-type-specific target power.
PWR_bio_smryAverage target power.
PWR_strata_ct_bio_smryCell type specific target power by gene expression stratification.
PWR_strata_bio_smryAverage target power by gene expression stratification.
ct_FDC_bio_smryCell type specific target FDC.
FDC_bio_smryAverage target FDC.
One S4 object.
data('quickPowerIBD')data('quickPowerIBD')
Plot false discovery cost results. This function plots false discovery cost results in a 2x1 panel. The illustration of each plot from left to right:
1: False discovery cost(FDC) by effect size, each line represents cell type. Sample size to be fixed at 10 if sample_size=10.
2: False discovery cost(FDC) by top effect size, each line represents sample size. FDC was the average value across cell types.
simulation_results |
A list of results produced by power evaluation functions. |
sample_size |
A numerical value indicating which sample size to be fixed. For example, 10 means when plotting the relationship between FDC and effect size, we fixed the scenario of sample size at 10. Default is 10. |
This function does not return a value. It generates a two-panel plot visualizing the false discovery cost (FDC) results. The first panel shows the FDC by effect size for each cell type at a fixed sample size (default is 10). The second panel illustrates the FDC by the top effect sizes, with each line representing a different sample size, averaged across cell types.
Wen Tang <[email protected]> Shilin Yu <[email protected]>
data(quickPowerGSE60424) ### Plot power results plotFDC(GSE60424Power,sample_size=10)data(quickPowerGSE60424) ### Plot power results plotFDC(GSE60424Power,sample_size=10)
This function plots all statistical power measurements in a 2x3 panel. The illustration of each plot from left to right and from up to bottom is as follows:
1: Statistical power by effect size, each line represents sample size. Statistical power was the average value across cell types.
2: Statistical power by effect size, each line represents cell type. Sample size is fixed at 10 if sample_size=10.
3: Statistical power by sample size, each line represents cell type. Effect size is fixed at 1 if effect.size=1.
4: Statistical power by strata, each line represents cell type. Sample size is fixed at 10 and effect size is fixed at 1 if sample_size=10 and effect.size=1.
5: Statistical power by strata, each line represents sample size. Statistical power was the average value across cell types and effect size is fixed at 1 if effect.size=1.
6: Statistical power by strata, each line represents effect size. Statistical power was the average value across cell types and sample size is fixed at 10 if sample_size=10.
simulation_results |
A list of results produced by power evaluation functions. |
effect.size |
A numerical value indicating which effect size is to be fixed. For example, 1 means when plotting the relationship between power and strata, we fixed the scenario of log fold change at 1. The default is 1. |
sample_size |
A numerical value indicating which sample size to be fixed. For example, 10 means when plotting the relationship between power and strata, we fixed the scenario of sample size at 10. The default is 10. |
This function generates a 2x3 panel plot visualizing various statistical power measurements, but does not return a programmable value. Each panel displays power metrics under different conditions such as effect size, sample size, and stratification, with lines representing either sample size, cell type, or effect size.
Wen Tang <[email protected]> Shilin Yu <[email protected]>
data(quickPowerGSE60424) ### Plot power results plotPower(GSE60424Power,effect.size=1,sample_size=10)data(quickPowerGSE60424) ### Plot power results plotPower(GSE60424Power,effect.size=1,sample_size=10)
This function plots all true discovery rate measurements in a 2x2 panel. The illustration of each plot is as follows:
1: True discovery rate(TDR) by top-rank genes, each line represents cell type. Sample size to be fixed at 10 and effect size to be fixed at 1 if sample_size=10 and effect.size=1.
2: True discovery rate(TDR) by top rank genes, each line represents effect size. TDR was the average value across cell types and sample size is fixed at 10 if sample_size=10.
3: True discovery rate(TDR) by top rank genes, each line represents sample size. TDR was the average value across cell types and effect size is fixed at 1 if effect.size=1.
4: True discovery rate(TDR) by effect size, each line represents sample size. TDR was calculated under the scenario of top rank gene equals 350.
simulation_results |
A list of results produced by power evaluation functions. |
effect.size |
A numerical value indicating which effect size is to be fixed. For example, 1 means when plotting the relationship between TDR and top rank genes for cell types or sample size, we fixed the scenario of log2 fold change at 1. The default is 1. |
sample_size |
A numerical value indicating which sample size to be fixed. For example, 10 means when plotting the relationship between TDR and top rank genes for cell types or effect size, we fixed the scenario of sample size at 10. The default is 10. |
This function creates a 2x2 panel plot showcasing various true discovery rate (TDR) measurements but does not return any values for further programmatic use. Each panel displays TDR analyses based on top rank genes, with lines representing different cell types, effect sizes, or sample sizes under specific conditions.
Wen Tang <[email protected]> Shilin Yu <[email protected]>
data(quickPowerGSE60424) ### Plot power results plotTDR(GSE60424Power,effect.size=1,sample_size=10)data(quickPowerGSE60424) ### Plot power results plotTDR(GSE60424Power,effect.size=1,sample_size=10)
This function quickly outputs pre-calculated power evaluation results from three datasets: (IAD, IBD, and ASD). The obtained results can be used to create plots from plot functions.
quickPower(data = "IAD")quickPower(data = "IAD")
data |
A character string specifying the dataset to be retrieved. Options are 'IAD', 'IBD', and 'ASD'. |
IAD: Whole transcriptome signatures of 6 immune cell subsets. The patients were drawn from subjects with a range of immune-related diseases.
IBD: Inflammatory Bowel Disease
ASD: Autism Spectrum Disorder.
ct_TDR_bio_smry |
Cell-type-specific target TDR. |
TDR_bio_smry |
Average target TDR. |
ct_PWR_bio_smry |
Cell-type-specific target power. |
PWR_bio_smry |
Average target power. |
PWR_strata_ct_bio_smry |
Cell type specific target power by gene expression stratification. |
PWR_strata_bio_smry |
Average target power by gene expression stratification. |
ct_FDC_bio_smry |
Cell type specific target FDC. |
FDC_bio_smry |
Average target FDC. |
Shilin Yu <[email protected]> Guanqun Meng <[email protected]>
# library(cypress) Quick_power <- quickPower(data = "IAD")# library(cypress) Quick_power <- quickPower(data = "IAD")
This function conducts simulations with various user-defined study design parameters, Users will need to provide SE object bulk data for parameter estimation purposes.
simFromData(INPUTdata = NULL, CT_index = NULL, CT_unk = FALSE, n_sim = 3, n_gene = 30000, DE_pct = 0.05, ss_group_set = c(10, 20, 50, 100), lfc_set = c(0, 0.5, 1, 1.5, 2), lfc_target = 0.5, fdr_thred = 0.1, DEmethod = "TOAST",BPPARAM=bpparam())simFromData(INPUTdata = NULL, CT_index = NULL, CT_unk = FALSE, n_sim = 3, n_gene = 30000, DE_pct = 0.05, ss_group_set = c(10, 20, 50, 100), lfc_set = c(0, 0.5, 1, 1.5, 2), lfc_target = 0.5, fdr_thred = 0.1, DEmethod = "TOAST",BPPARAM=bpparam())
INPUTdata |
The input SE (SummarizedExperiment) object should contain a count matrix, study design, and an optional cell type proportion matrix. The study design should have a column named ‘disease’, where the control by 1 and the case is indicated by 2. If provided, the cell type proportion matrix should sum to 1 for each sample. The cell type proportion matrix is optional, the CT_unk should be True if the user did not provide this matrix |
CT_index |
Column index for cell types proportion matrix in Coldata, the
input can also be a single number (>3) when the |
CT_unk |
Logical flag indicating whether unknown cell types are present.
Defaults to |
n_sim |
The total number of iterations users wish to conduct. Default to 3. In simulation results, it is set to 20. |
n_gene |
Total number of genetic features users with to conduct. Default to 30000. Must be greater than or equal to 1000. |
DE_pct |
Percentage of DEG on each cell type. Default to 0.05. |
ss_group_set |
Sample sizes per group users wish to simulate. The length should be less than or equal to 5. Default to 10,20,50. |
lfc_set |
Effect sizes users wish to simulate. The length should be less than or equal to 5. Default to 0,0.5,1,1.5. |
lfc_target |
Target effect size, should be greater than or equal to 0. The absolute LFC lower than this values will be treated as None-DEGs. Default to 0.5 |
fdr_thred |
Adjusted p value threshold. The parameter value should be within the range (0, 1). Default to 0.1 |
DEmethod |
Differential expression (DE) methods available include 'TOAST', 'DESeq2', and 'CeDAR'. The default method is 'TOAST' |
BPPARAM |
An instance of |
One SummarizedExperiment object containing the following elements:
countsA gene expression value dataset
colDataSample meta-data. The first column is the group status (i.e. case/ctrl) named as 'disease', and the second column is the subject ID. The remaining are the cell type proportions of all samples. The user could also input the Column index for cell types proportion matrix in Coldata. Example: CT_index= 3:8
ct_TDR_bio_smry |
Cell-type-specific target TDR. |
TDR_bio_smry |
Average target TDR. |
ct_PWR_bio_smry |
Cell-type-specific target power. |
PWR_bio_smry |
Average target power. |
PWR_strata_ct_bio_smry |
Cell type specific target power by gene expression stratification. |
PWR_strata_bio_smry |
Average target power by gene expression stratification. |
ct_FDC_bio_smry |
Cell type specific target FDC. |
FDC_bio_smry |
Average target FDC. |
Shilin Yu <[email protected]> Guanqun Meng <[email protected]>
data(ASD_prop_se) result <- simFromData(INPUTdata = ASD_prop, CT_index = (seq_len(6) + 2), CT_unk = FALSE, n_sim = 2,n_gene = 1000, DE_pct = 0.05, ss_group_set = c(8,10), lfc_set = c(1, 1.5))data(ASD_prop_se) result <- simFromData(INPUTdata = ASD_prop, CT_index = (seq_len(6) + 2), CT_unk = FALSE, n_sim = 2,n_gene = 1000, DE_pct = 0.05, ss_group_set = c(8,10), lfc_set = c(1, 1.5))
This function conducts simulations with various user-defined study design parameters, including but not limited to sample size, and log fold change.)
simFromParam(n_sim = 3, n_gene = 30000, DE_pct = 0.05, ss_group_set = c(10, 20, 50, 100), lfc_set = c(0, 0.5, 1, 1.5, 2), sim_param = "IAD", lfc_target = 0.5, fdr_thred = 0.1, DEmethod = "TOAST", BPPARAM=bpparam())simFromParam(n_sim = 3, n_gene = 30000, DE_pct = 0.05, ss_group_set = c(10, 20, 50, 100), lfc_set = c(0, 0.5, 1, 1.5, 2), sim_param = "IAD", lfc_target = 0.5, fdr_thred = 0.1, DEmethod = "TOAST", BPPARAM=bpparam())
n_sim |
The total number of iterations users wish to conduct. Default to 3. In simulation results, it is set to 20. |
n_gene |
Total number of genetic features users with to conduct. Default to 30000. Must be greater than or equal to 1000. |
DE_pct |
Percentage of DEG on each cell type. Default to 0.05. |
ss_group_set |
Sample sizes per group users wish to simulate. The length should be less than or equal to 5. Default to 10,20,50. |
lfc_set |
effect sizes users wish to simulate. The length should be less than or equal to 5. Default to 0,0.5,1,1.5. |
sim_param |
Users specify which embedded simulation parameters they wish to use. By default set to 'IAD', which is a cell line specific data. Other options include 'IBD' data, and 'ASD' data |
lfc_target |
Target effect size, should be greater than or equal to 0. The absolute LFC lower than this value will be treated as None-DEGs. Default to 0.5 |
fdr_thred |
Adjusted p value threshold. The parameter value should be within the range (0, 1). Default to 0.1 |
DEmethod |
Differential expression (DE) methods available include 'TOAST', 'DESeq2', and 'CeDAR'. The default method is 'TOAST' |
BPPARAM |
An instance of |
GSE60424: Immune-related disease (IAD) study. Whole transcriptome signatures of 6 immune cell subsets. The patients were drawn from subjects with a range of immune-related diseases.
IBD: data in pediatric inflammatory bowel disease(IBD) study
ASD: data in a large autism spectrum disorder (ASD) study
ct_TDR_bio_smry |
Cell-type-specific target TDR. |
TDR_bio_smry |
Average target TDR. |
ct_PWR_bio_smry |
Cell-type-specific target power. |
PWR_bio_smry |
Average target power. |
PWR_strata_ct_bio_smry |
Cell type specific target power by gene expression stratification. |
PWR_strata_bio_smry |
Average target power by gene expression stratification. |
ct_FDC_bio_smry |
Cell type specific target FDC. |
FDC_bio_smry |
Average target FDC. |
Shilin Yu <[email protected]> Guanqun Meng <[email protected]>
data(quickParaGSE60424) result <- simFromParam(sim_param="IAD",n_sim = 2,DE_pct = 0.05,n_gene = 1000, ss_group_set = c(8, 10), lfc_set = c(1, 1.5), lfc_target = 0.5, fdr_thred = 0.1)data(quickParaGSE60424) result <- simFromParam(sim_param="IAD",n_sim = 2,DE_pct = 0.05,n_gene = 1000, ss_group_set = c(8, 10), lfc_set = c(1, 1.5), lfc_target = 0.5, fdr_thred = 0.1)