NEWS.md
read.pcadapt()
generates bed
files instead of pcadapt
files.
Computation of PCA is now based on R package RSpectra
.
Missing values are handled by specifying matrix-vector operations in RSpectra
that accounts for missing values.
Includes LD thinning to compute PCs.
No more dependency to R package RcppArmadillo.
For Pooled-seq data, use Mahalanobis distances based on PCA loadings, no more simulations of individual genotypes.
Switch from C/Lapack to Rcpp/RcppArmadillo.
pcadapt()
can take genotype matrices as input.
Modified code for binomial sampling.
pcadapt()
argument clean.files
is now deprecated.
pcadapt()
argument output.filename
is now deprecated.
read.pcadapt()
argument local.env
is now deprecated.
Latest update of R package vcfR
taken into account.
Option type = "vcfR"
has been added to read.pcadapt()
to overcome some conversion issues occurring with VCF files.
Argument transpose
is now deprecated. Read section A for more details.
get.pc()
has been added. For each SNP, it returns the most correlated principal component.Function read4pcadapt()
is now deprecated, it is now called read.pcadapt()
.
Using the pop
option when plotting scores now provides the color legend.
All analyses are now included in the R package. Users should not use the C software PCAdapt fast anymore.
Big datasets can be handled directly within the R session.
read4pcadapt()
now converts files to the pcadapt
format.
The first argument of pcadapt()
can be either a small genotype matrix or the output of read4pcadapt()
.
The Mahalanobis distance is now estimated from the z-scores rather than the loadings.
Make sure you have downloaded the latest version of the C software PCAdapt (last updated on February 11, 2016).
The scaling of the SNP before computing PCA has been changed. Instead of using standard deviation, we now use the square root of p(1−p) (haploid data) or of 2p(1−p) (diploid data) where p is the minimum allele frequency.
Bug fix: the genomic inflation factor has been corrected when K=1
.
Bug fix: a problem due to high proportion of missing data slowing the program has been fixed.
Argument "minmaf"
has been replaced with "min.maf"
.