Bug fix

lvclark · lvclark · commit 1bea7c88aa2b · 2024-12-20T12:49:36.000-08:00
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,6 +1,6 @@
 Package: polysat
-Version: 1.7-7
-Date: 2022-08-20
+Version: 1.7-8
+Date: 2024-12-20
 Title: Tools for Polyploid Microsatellite Analysis
 Authors@R: c(person("Lindsay V.", "Clark", email = "Lindsay.Clark@seattlechildrens.org", role = c("aut", "cre"),
                     comment = c(ORCID = "0000-0002-3881-9252")),
diff --git a/NEWS b/NEWS
@@ -1,3 +1,7 @@
+Changes in version 1.7-8
+
+* Bug fixed in PIC function.
+
 Changes in version 1.7-7
 
 * Maintainer email address updated.
diff --git a/R/population_stats.R b/R/population_stats.R
@@ -896,7 +896,7 @@ PIC <- function(freqs, pops=row.names(freqs),
       } else {
         genomes <- freqs[,"Genomes"]
       }
-      meanfreq <- apply(freqs[,lcol], 2, function(x) weighted.mean(x, w = genomes))
+      meanfreq <- apply(freqs[,lcol, drop = FALSE], 2, function(x) weighted.mean(x, w = genomes))
       results["Overall",L] <- pic(meanfreq)
     }
   }

Original file line number	Diff line number	Diff line change
`@@ -896,7 +896,7 @@ PIC <- function(freqs, pops=row.names(freqs),`
`896`	`896`	`} else {`
`897`	`897`	`genomes <- freqs[,"Genomes"]`
`898`	`898`	`}`
`899`		`- meanfreq <- apply(freqs[,lcol], 2, function(x) weighted.mean(x, w = genomes))`
	`899`	`+ meanfreq <- apply(freqs[,lcol, drop = FALSE], 2, function(x) weighted.mean(x, w = genomes))`
`900`	`900`	`results["Overall",L] <- pic(meanfreq)`
`901`	`901`	`}`
`902`	`902`	`}`