some functions didn't work in my dataset

31 views
Skip to first unread message

changs...@gmail.com

unread,
Nov 25, 2020, 1:27:27 AM11/25/20
to R/qtl2 discussion
Dear Karl,

Sorry to bother you, I got some problems with my rislef4 datasets. 

1. zip_datafiles() didn't work without any warning message.
2. genoprob_to_alleleprob() return a empty list without any warning message.

Here is my code:
```
library(qtl2)

setwd("C:/Users/wpf/Desktop/project/MAGIC4w/")

zip_datafiles("magic4w.json", quiet = F)

rice <- read_cross2("magic4w.json", quiet = F)

pmap <- rice$pmap

prcl <- calc_genoprob(rice, map = pmap, error_prob = 0.0001,
                      map_function = "c-f", lowmem = F,
                      quiet = F, cores = 1)

apr <- genoprob_to_alleleprob(prcl, quiet = F)
```

my R.version

platform       x86_64-w64-mingw32          
arch           x86_64                      
os             mingw32                     
system         x86_64, mingw32             
status                                     
major          4                           
minor          0.2                         
year           2020                        
month          06                          
day            22                          
svn rev        78730                       
language       R                           
version.string R version 4.0.2 (2020-06-22)
nickname       Taking Off Again  

If you need my dataset for testing, I will email you.

Thanks, advance!

Yours Sincerely,
Pengfei Wang  

Karl Broman

unread,
Nov 25, 2020, 1:36:19 AM11/25/20
to R/qtl2 discussion
genoprob_to_alleleprob() is only used for heterozygous populations. In an riself4 population, the genotype probabilities are the same as the allele probabilities, so there would be no reason to all this function.

I can’t say anything about the zip_datafiles() without seeing the actual warning messages. 

karl

changs...@gmail.com

unread,
Nov 25, 2020, 1:47:11 AM11/25/20
to R/qtl2 discussion
I truly appreciate your timely help. So the premise of "riself" and "risib" is infinite selfing and sib mating without heterozygous sites, is it?

Karl Broman

unread,
Nov 25, 2020, 9:18:09 AM11/25/20
to R/qtl2 discussion
R/qtl2 assumes no heterozygosity.

karl

changs...@gmail.com

unread,
Nov 25, 2020, 9:37:54 AM11/25/20
to R/qtl2 discussion
Thank you. Wish you a good day!

Pengfei Wang
Reply all
Reply to author
Forward
0 new messages