Package ‘EigenCorr’ August 11, 2011 Version 0.2 Date 2011-08-10 License GPL version 2 or newer Description Compute p-values of EigenCorr1, EigenCorr2 and Tracy-Widom to select principal components for adjusting population stratification. Title EigenCorr Author Seunggeun, Lee Maintainer Seunggeun, Lee Depends R (>= 2.0.0) Suggests URL http://www.r-project.org

R topics documented: EigenCorr . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . EigenCorr_Data . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Read_PCA . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Index

1 3 3 5

1

2

EigenCorr

Compute scores and p-values of EigenCorrs

EigenCorr

Description EigenCorr computes p-values of EigenCorr1, EigenCorr2 and Tracy-Widom, which can be used for selecting principal components to adjust population stratification. Usage EigenCorr( Y,PC,EVAL,N.Pval=100,N.Simu=10000) Arguments Y

an output phenotype vector.

PC

a matrix of principal components. The kth column should be the kth principal components. The order of rows and columns should be matched with Y and EVAL vectors. It need not to have all principal components.

EVAL

an vector of all eigenvalues. It should have all eigenvalues!!

N.Pval

an integer, the number of principal components whose p-values will be computed. If the number of principal components is smaller than N.Pval, p-values of all principal components (except last 2) will be computed (default=100).

N.Simu

an integer, the number of simulations to generate NULL distribution (default=10000).

Value EigenCorr returns an dataframe object with 5 columns. Gamma Correlations between Y and PCs EigenCorr_Score EigenCorr scores. EigenCorr score is defined as multiplication of square correlation and eigenvalue. EigenCorr1

P-values of EigenCorr1. They are computed based on the simulated null distribution.

EigenCorr2

P-values of EigenCorr2.

TW

P-values of eigenvalues based on the Tracy-Widom distribution.

Author(s) Seunggeun Lee References Lee, S., Wright, F.A. and Zou, F. (2010) Control of population stratification by correlation-selected principal components. Biometrics, in press.

EigenCorr_Data

3

Examples library(EigenCorr) data(EigenCorr_Data)

# EigenCorr example dataset

Y<-EigenCorr_Data$Y PC<-EigenCorr_Data$PC EVAL<-EigenCorr_Data$EVAL EigenCorr( Y,PC,EVAL,N.Pval=50)

EigenCorr_Data

Example data for EigenCorr

Description This is an example dataset for EigenCorr. It contains outcome phenotypes, principal components and eigenvalues. Usage EigenCorr_Data Format EigenCorr_Data contains the following objects: Y a vector of the outcome phenotypes PC a matrix of the principal components EVAL a vector of the eigenvalues

Read_PCA

Read EigenStrat and ShrukenPCA output files

Description Read_PCA reads EigenStrat and ShrukenPCA main output file. This file is formated as: 1) First line indicates number of principal components (k), 2) Next k lines have top k eigenvalues, and 3) N additional lines for k principal components where N is the number of samples. Usage Read_PCA(filename)

4

Read_PCA

Arguments filename

an EigenStrat or ShrukenPCA output file.

Value Read_PCA returns a list object. K

an integer, the number of principal components.

EVAL

a eigen value vector.

PC

a matrix of K principal components. The kth column is the kth principal component.

Author(s) Seunggeun Lee

Index EigenCorr, 1 EigenCorr_Data, 3 Read_PCA, 3

5

Package 'EigenCorr'

Aug 11, 2011 - The kth column should be the kth principal components. The order of rows and ... Example data for EigenCorr. Description. This is an example ...
Download PDF
82KB Sizes 1 Downloads 257 Views
Report

Recommend Documents

Package 'EigenCorr'
Aug 11, 2011 - License GPL version 2 or newer. Description Compute p-values of EigenCorr1, EigenCorr2 and Tracy-Widom to select principal components for adjusting population stratification. Title EigenCorr. Author Seunggeun, Lee . Maintainer Seunggeu
Package 'EigenCorr' - umich.edu and www-personal
Aug 11, 2011 - Maintainer Seunggeun, Lee . Depends R (>= 2.0.0). Suggests. URL http://www.r-project.org. R topics documented:.
Package 'TeachingSampling' February 14, 2012 Type Package Title ...
Feb 14, 2012 - Creates a matrix of domain indicator variables for every single unit in ... y Vector of the domain of interest containing the membership of each ...
Package 'MethodEvaluation' - GitHub
Feb 17, 2017 - effects in real data based on negative control drug-outcome pairs. Further included are .... one wants to nest the analysis within the indication.
Package 'CohortMethod' - GitHub
Jun 23, 2017 - in an observational database in the OMOP Common Data Model. It extracts the ..... Create a CohortMethod analysis specification. Description.
Package 'hcmr' - GitHub
Effective green time to cycle length ratio. P ... Two-Lane Highway - Base Percent Time Spent Following .... Passenger-Car Equivalent of Recreational Vehicles:.
Package 'CaseCrossover' - GitHub
Apr 21, 2017 - strategies for picking the exposure will be tested in the analysis, a named list of .... A data frame of subjects as generated by the function ...
Package 'SelfControlledCaseSeries' - GitHub
Mar 27, 2017 - 365, minAge = 18 * 365, maxAge = 65 * 365, minBaselineRate = 0.001,. maxBaselineRate = 0.01 .... Modeling and Computer Simulation 23, 10 .... function ggsave in the ggplot2 package for supported file formats. Details.
Package 'RMark'
Dec 12, 2012 - The RMark package is a collection of R functions that can be used as an interface to MARK for analysis of capture-recapture data. Details.
The PythonTeX package
It would be nice for the print statement/function,6 or its equivalent, to automatically return its output within the LATEX document. For example, using python.sty it is .... If you are installing in TEXMFLOCAL, the paths will have an additional local
package management.key - GitHub
Which version of Faker did our app depend on? If we run our app in a year and on a different machine, will it work? If we are developing several apps and they each require different versions of Faker, will our apps work? Page 6. Gem Management with B
Package No.2 Rs. 2,25000 Package No.1 -
Mentioned as Education Partner in AIESEC Chennai. Newsletter. ✓. Logo Visibility on Facebook for the year 2012. ✓. Detailed database of all youth approached ...
Package 'cmgo' - GitHub
Aug 21, 2017 - blue all Voronoi segments, b) in red all segments fully within the channel polygon, c) in green all ..... if [TRUE] the plot will be saved as pdf.
Package 'EmpiricalCalibration' - GitHub
study setup. This empirical null distribution can be used to compute a .... Description. Odds ratios from a case-control design. Usage data(caseControl). Format.
Package 'OhdsiRTools' - GitHub
April 7, 2017. Type Package. Title Tools for Maintaining OHDSI R Packages. Version 1.3.0. Date 2017-4-06. Author Martijn J. Schuemie [aut, cre],. Marc A.
Package 'FeatureExtraction' - GitHub
deleteCovariatesSmallCount = 100, longTermDays = 365, ..... Description. Uses a bag-of-words approach to construct covariates based on free-text. Usage.
Package 'EvidenceSynthesis' - GitHub
Mar 19, 2018 - This includes functions for performing meta-analysis and forest plots. Imports ggplot2 (>= 2.0.0),. gridExtra, meta,. EmpiricalCalibration. License Apache License 2.0. URL https://github.com/OHDSI/EvidenceSynthesis. BugReports https://
Package 'RNCEP'
A numeric argu- ment passed to the when2stop list indicates a distance from the end.loc in kilometers at which to stop the simulation. The simulation will end ...... This provides an indication of the precision of an interpolated result described in
Package 'forecast'
Oct 4, 2011 - Depends R (>= 2.0.0), graphics, stats, tseries, fracdiff, zoo. LazyData yes .... Largely wrappers for the acf function in the stats package. The main ...
Package 'IcTemporalPatternDiscovery' - GitHub
Nov 25, 2015 - exposureOutcomePairs = data.frame(outcomeId = c(196794, ... strategies for picking the exposure will be tested in the analysis, a named list of.
Package 'sirt'
Oct 7, 2016 - Rasch type models using the generalized logistic link function (Stukel, 1988) ...... The item parameters can be reparametrized as ai = exp [(-δi + τi)/2] and bi = exp [(δi + τi)/2]. .... Optional file name for sinking the summary in
Package 'CDM'
Dec 5, 2016 - A vector, a matrix or a data frame of the estimated parameters for the fitted model. ..... ECPE dataset from the Templin and Hoffman (2013) tutorial of specifying cognitive ...... An illustration of diagnostic classification modeling in
Package 'miceadds'
Dec 14, 2016 - DESCRIPTION miceadds package ... The miceadds package also includes some functions R utility functions ..... Applied missing data analysis.