Comments (3)
Do you have an example that demonstrates this? Because it seems to be using it appropriately.
Consider this example:
library(qtl)
data(hyper)
hyper <- calc.genoprob(hyper[1:2,])
hyper$pheno[,3] <- sqrt(hyper$pheno[,1])
set.seed(10347067)
out1 <- scantwopermhk(hyper, pheno.col=1, n.perm=3)
set.seed(10347067)
out3 <- scantwopermhk(hyper, pheno.col=3, n.perm=3)
set.seed(10347067)
out1b <- scantwopermhk(hyper, pheno.col=1, n.perm=3)
all(unlist(out1) == unlist(out1b))
all(unlist(out1) != unlist(out3))
from qtl.
Apologies for a less elegant example, but:
library(qtl)
load("goughF2_simple_v4.RData")
reduced_f2 <- calc.genoprob(f2, step = 1, error.prob = 0.0002, map.function = 'c-f')
reduced_f2 <- reduce2grid(reduced_f2)
reduced_f2 <- subset(reduced_f2, ind = reduced_f2$pheno$sex == "M")
nrow(reduced_f2$pheno)
out1 <- scantwopermhk(reduced_f2, pheno.col = 'wk5', assumeCondIndep = TRUE, n.perm = 1)
out2 <- scantwopermhk(reduced_f2, pheno.col = 2, assumeCondIndep = TRUE, n.perm = 1)
summary(out1);summary(out2)
wk1 (1 permutations)
full fv1 int add av1 one
5% 7.7 5.83 5.05 5.18 3.9 2.35
10% 7.7 5.83 5.05 5.18 3.9 2.35
wk1 (1 permutations)
full fv1 int add av1 one
5% 8.69 7.6 5.58 5.09 2.61 2.76
10% 8.69 7.6 5.58 5.09 2.61 2.76
from qtl.
Ah; it's a problem in the pasting over of the appropriate column name. Got it fixed in the devel branch.
from qtl.
Related Issues (20)
- scanone perms with 4-way cross giving all 0's HOT 3
- Change read.cross warning about length of arg alleles
- pull.rf gives error with duplicate marker names
- GH Action check failing on MacOS HOT 2
- stepwiseqtl has problems with X chr and sex covariate HOT 1
- Problem in scantwo permutations HOT 1
- addcovar arg for cim()
- addtoqtl should clear attributes HOT 1
- Error: "long vectors (argument 5) are not supported in .C" HOT 2
- switch.alleles should give error with results of markerlrt HOT 1
- reduce2grid error HOT 2
- locationXO and countXO for BCsFt
- addint bug when X chr QTL and missing phenotypes HOT 2
- summary.scanone error with threshold as a vector
- Improper header in pheno dataframe with single column pheno file HOT 3
- Warning on CRAN
- Additional warnings on CRAN
- Error handling of pheno.col in cim()
- Compilation error on R-devel HOT 3
- Calloc/Realloc can no longer be used? HOT 1
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 qtl.