File: testaugmentation.R

package info (click to toggle)
qtl 1.23-16-1
  • links: PTS, VCS
  • area: main
  • in suites: wheezy
  • size: 11,580 kB
  • sloc: ansic: 11,586; cpp: 3,085; ruby: 193; sh: 171; makefile: 1
file content (30 lines) | stat: -rw-r--r-- 720 bytes parent folder | download | duplicates (9)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
# Test augmentation with MQM
#
# Note: the full version of this test has moved to ./contrib/bin/rtest, 
#       as it takes a long time to run. The full version can be run with:
#
#   cd contrib/bin
#   rm CMakeCache.txt ; cmake -DTEST_R=TRUE
#   make testR

library(qtl)

set.seed(1000)              

version = mqm_version()
cat("R/qtl=",version$RQTL)
cat("R-MQM=",version$RMQM)
cat("MQM=",version$MQM)

testaugmentation <- function(cross, ...){
  crossML <- mqmaugment(cross, ...)

  res1 <- mqmscan(crossML,logtransform=TRUE)
  list(res1)
}

data(listeria)
r <- testaugmentation(listeria)
if(!round(r[[1]][3,3],3)==0.307) stop("Listeria ML dataaugmentation error")

cat("testaugmentation.R, tests succesfully run!")