KonradZych / phenotypes2genotypes

Package for creating genetic map using gene expression data.
6 stars 5 forks source link

Error from Cross.DeNovo #58

Closed DannyArends closed 11 years ago

DannyArends commented 11 years ago

using this call:

 crossSkeleton   <- cross.denovo(population, n.chr=5, map="genetic")

results in:

 saving data into cross object done in: 28.657 seconds.
 Error in if (rowMinimums[i] > 0.5) { : 
 missing value where TRUE/FALSE needed
 In addition: Warning messages:
 1: In summary.cross(cross) :
   Some chromosomes > 1000 cM in length; there may be a problem with the genetic map.
   (Perhaps it is in basepairs?)
 2: In est.rf(cross) : Alleles potentially switched at markers 
   8837 34796 108535 173733 209835 211485 366014 366023 410899 554496 627262 712657 780566 839440 880003 906005 907074 960552 1016044 1120339 1124617 1125275 1129220 1177346 1177398 1178055 1330721 1369026 1396700 1464959 1517513 1724058 13884 62430 139634 268463 273538 535000 587917 599335 637041 853855 871168 890249 950031 952539 955105 1044391 1051033 1144602 1188243 1195939 1196683 1248470 1263007 1295550 1320059 1327503 1421858 1428926 1450675 1559769 1587269 1616628 1618488 1631010 1642692 1688633 1782628

This can be traced back to:

 ~/Github/Rpackages/phenotypes2genotypes/R (master) $ grep "rowMinimums" *.R
 assignLinkageGroups.R:  rowMinimums    <- apply(dataRf, 1, min)
 assignLinkageGroups.R:    if(rowMinimums[i] > 0.5){
 assignLinkageGroups.R:      cat("dropping", rownames(dataRf)[i], "min:", 

rowMinimums[i])