diff --git a/DESCRIPTION b/DESCRIPTION index b1ba150..c45ecb9 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,11 +1,11 @@ Package: sciClone Title: Detect subclones from genomic sequence data -Version: 1.1.0 +Version: 1.1.1 Author: Chris Miller, Brian White, Nathan Dees Description: Identifies clonal entities from sequenced populations of cells, using read count information. Maintainer: Chris Miller LazyLoad: yes License: Apache License 2.0 | file LICENSE Depends: IRanges, bmm, rgl, RColorBrewer, ggplot2, grid, plotrix, - methods, NORMT3, MKmisc, TeachingDemos -Packaged: 2016-03-20 11:07:15 PM; cmiller + methods, NORMT3, MKmisc, TeachingDemos, dplyr +Packaged: 2023-10-11 05:03:14 PM; cmiller diff --git a/R/sciClone.R b/R/sciClone.R index 1c12804..06f18c2 100644 --- a/R/sciClone.R +++ b/R/sciClone.R @@ -9,6 +9,8 @@ sciClone <- function(vafs, copyNumberCalls=NULL, regionsToExclude=NULL, copyNumberMargins=0.25, maximumClusters=10, annotation=NULL, doClusteringAlongMargins=TRUE, plotIntermediateResults=0){ + suppressPackageStartupMessages(library(dplyr)) + if(verbose){print("checking input data...")} #how many samples do we have? @@ -227,8 +229,9 @@ sciClone <- function(vafs, copyNumberCalls=NULL, regionsToExclude=NULL, vafs.1d.merged.cn2 = cbind(vafs.merged.cn2.orig,cluster=marginalClust[[i]]$cluster.assignments) vafs.1d.merged = merge(vafs.merged.orig,vafs.1d.merged.cn2, by.x=c(1:length(vafs.merged.orig)), by.y=c(1:length(vafs.merged.orig)),all.x=TRUE) + ##sort by chr, st - vafs.1d.merged = vafs.1d.merged[order(vafs.1d.merged[,1,drop=FALSE], vafs.1d.merged[,2,drop=FALSE]),] + vafs.1d.merged = dplyr::arrange(vafs.1d.merged, chr, st) vafs.1d[[i]] = vafs.1d.merged } }