Comments (4)
I'm working on the function to estimate round robin frequencies in place, but I'm having trouble passing ASAN at my comparison function for qsort.
The error it gives me:
mlg_counter.c:61:10: runtime error: member access within null pointer of type 'struct mask'
All I want to do is cast it into a char *
. @JonahBrooks, do you have any pointers (pun VERY intended)?
from poppr.
The issue above has been fixed and the results are fantastic as compared to doing the same thing in base R:
Setup:
# Function for counting the number of round robin mlgs per locus
rr <- function(x){
mask <- seq(ncol(x))
res <- vector(mode = "integer", length = ncol(x))
for (i in mask){
res[i] <- nrow(unique(x[, -i, drop = FALSE]))
}
return(res)
}
x <- matrix(sample(0:2, 6660, rep = TRUE), nrow = 666, ncol = 10)
library('microbenchmark')
library('poppr')
res <- microbenchmark(rr(x), .Call("mlg_round_robin", x, PACKAGE = "poppr"))
res
Unit: milliseconds
expr min lq mean median uq max neval
rr(x) 199.243444 220.686884 241.518477 229.980039 245.227290 569.22988 100
.Call("mlg_round_robin", x, PACKAGE = "poppr") 4.213625 4.396508 6.976769 4.925717 5.974037 50.60142 100
from poppr.
I've added a function to estimate round robin allele frequencies in the mlg-counter branch called rraf()
from poppr.
I've added both pgen()
and psex()
and I've documented Stacy's name on the functions, now I just need to add her name to the Authors@R field in the Description.
from poppr.
Related Issues (20)
- Genetic distance calculation yields NAs HOT 2
- Fix inefficient implementation of boot.bruvo HOT 3
- Interpretation of negative values in AMOVA HOT 3
- Error with reverse dependency checks for phangorn HOT 12
- genind2genalex() produces all zero genotypes with some SNP data. HOT 1
- bitwise.ia() degraded performance on windows HOT 1
- Clade rooting abnormally with updated HOT 4
- Error in bruvo.boot HOT 4
- CRAN test failure in R-devel on Debian
- Cannot re-level strata of a snpclone object HOT 1
- DUE 2021-09-21: Packages including omp.h after Rincludes.h
- Feature request: separation/removal of sf package requirement HOT 5
- Optional disabling of statistics in poppr table HOT 6
- threads not being respected and defaulting to 1 HOT 1
- PAIRWISE FST Calculation HOT 1
- Find All Duplicate MLGs HOT 4
- [read.genalex()]: Detect semicolon-separated files and throw error with correct command
- IS there a way to make a table indicated the multilocus genotypes for each isolate? where is that info stored? HOT 1
- Error Message from R using read.genalex HOT 1
- CRAN Errors to fix by 2023-03-22 HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from poppr.