# NOT RUN {
# Monocentrics - Beginning with two data.frames with column OTU
require(idiogramFISH)
dfOfChrSize$OTU<-1
dfOfMarks2$OTU<-1
# check missing data
dfOfMarks2[which(is.na(dfOfMarks2$markSize) & dfOfMarks2$chrRegion %in% c("p","q") ),]
# markSize and markDistCen missing
# add markDistCen
dfOfMarks2$markDistCen <- ifelse( is.na(dfOfMarks2$markDistCen) &
dfOfMarks2$chrRegion %in% c("p","q"),
0,
dfOfMarks2$markDistCen)
# add markSize
XshortSize <- dfOfChrSize[which(dfOfChrSize$chrName=="X"),]$shortArmSize
dfOfMarks2[which(dfOfMarks2$chrName=="X" & dfOfMarks2$chrRegion=="p") ,]$markSize<-XshortSize
# add column chrSize if missing
dfOfChrSizeWithChrSize <- armRatioCI(dfOfChrSize)
# data.frame of chr. size to list
listOfdfChr <- idiogramFISH:::dfToListColumn(dfOfChrSizeWithChrSize)
ti<-posCalc(dfOfMarks2,listOfdfChr)
as.list(ti[[1]][,1])
posCalc(dfOfMarks2,listOfdfChr, result="list")
# holocentrics example
dfChrSizeHolo$OTU<-2
dfMarkPosHolo$OTU<-2
# data.frame of chr. size to list
listOfdfChrHolo <- idiogramFISH:::dfToListColumn(dfChrSizeHolo)
ti2<-posCalc(dfMarkPosHolo,listOfdfChrHolo)
ti2
posCalc(dfMarkPosHolo,listOfdfChrHolo, result="list")
as.list(ti2[[1]][,1] )
# }
Run the code above in your browser using DataLab