CGHpower Results

For citing, please use (manuscript submitted).

Dataset Information

Name: Simulation (30 arrays)

Group A: Group 1 (15)
group1.microarray025, group1.microarray011, group1.microarray024, group1.microarray007, group1.microarray026, group1.microarray016, group1.microarray022, group1.microarray007.1, group1.microarray018, group1.microarray013, group1.microarray020, group1.microarray004, group1.microarray015, group1.microarray011.1, group1.microarray025.1

Group B: Group 2 (15)
group2.microarray019, group2.microarray025.2, group2.microarray026.1, group2.microarray004.1, group2.microarray015.1, group2.microarray021, group2.microarray005, group2.microarray001, group2.microarray001.1, group2.microarray025.3, group2.microarray008, group2.microarray017, group2.microarray001.2, group2.microarray012, group2.microarray008.1

Unused: (0)

Power Calculations Estimates of average power

Number of probes:40397
Number of regions:461
Estimate of gamma:0.45

Sample SizeAverage power
83.7%
1511%
2317.7%
3023.6%
6040.5%
9050.1%
12056.7%
15061.1%
18064.1%
21067%
24069.2%
27070.7%
30072.2%
33073.6%

Average power plotted as a function of sample size.

aCGH Profiles and Chromosomal Regions

aCGH Profiles
Copy number profiles produced with CGHcall. The mean probability of losses is shown in red, and the values can be red from the Y axis. Mean probability of gains is in green, and the values are 1 - the value from the Y axis.
Chromosomal Regions
Output from CGHregions. Chromosomes are plotted individually, and each bump represents a breakpoint between regions. The loss/gain frequencies are shown in red/green.

Diagnostic Plots and Parameter Estimation

To help interpreting the these plots, comparison can be made to the evaluation data sets.
Estimates of G
The two estimators of G should not be in too much of a disagreement compared to each other. If the difference is severe, the quality of parameter estimation is questionable, and so is the reliability of the power calculations.
Distribution of p values
The density of the p values should increase for small values and the function should be convex.
Estimation of gamma
Gamma is the proportion of non-differentially behaving regions. The plot shows the distance between the two estimators of G, and the value of gamma is chosen so that this difference is minimized. The function should have a minimum somewhere in the middle. A minimum very close to 0 is a sign of problems in parameter estimation.
Estimation of lambda
Effect size is the difference between the groups. The function can have one ore more peaks, depending on the particular data set.
Skewness
Skewness of the RWLRs is plotted and that of a normal distribution is superimposed in red. RWLRs are assumed to be approximately normally distributed, so a large deviation might negatively affect the performance of the method.
Kurtosis
Kurtosis of the RWLRs is plotted and that of a normal distribution is superimposed in red. RWLRs are assumed to be approximately normally distributed, so a large deviation might negatively affect the performance of the method.