To find out how many AIMS have above 0.95 values (high idas anc) or less than 0.05 (high melissa anc) values. And the same for Z chromosome.
sum(bld3_mp_f$mean > 0.95)
sum(bcr3_mp_f$mean > 0.95)
sum(bic3_mp_f$mean > 0.95)
sum(btb3_mp_f$mean > 0.95)
sum(frc3_mp_f$mean > 0.95)
sum(pin3_mp_f$mean > 0.95)
sum(psp3_mp_f$mean > 0.95)
sum(rdl3_mp_f$mean > 0.95)
sum(rnv3_mp_f$mean > 0.95)
sum(bcr3_mp_f$mean < 0.05)
sum(bic3_mp_f$mean < 0.05)
sum(bld3_mp_f$mean < 0.05)
sum(btb3_mp_f$mean < 0.05)
sum(frc3_mp_f$mean < 0.05)
sum(pin3_mp_f$mean < 0.05)
sum(psp3_mp_f$mean < 0.05)
sum(rdl3_mp_f$mean < 0.05)
#linkage
zrdl<-rdl3_mp_f[rdl3_mp_f$linkage == 23,]
length(which(zrdl$mean > 0.95)) #54
length(which(zrdl$mean < 0.05)) #1
zbcr<-bcr3_mp_f[bcr3_mp_f$linkage == 23,]
length(which(zbcr$mean > 0.95)) #54
length(which(zbcr$mean < 0.05)) #1
zbic<-bic3_mp_f[bic3_mp_f$linkage == 23,]
length(which(zbic$mean > 0.95)) #54
length(which(zbic$mean < 0.05)) #1
zbld<-bld3_mp_f[bld3_mp_f$linkage == "Z",]
length(which(zbld$mean > 0.95)) #54
length(which(zbld$mean < 0.05)) #1
zbtb<-btb3_mp_f[btb3_mp_f$linkage == 23,]
length(which(zbtb$mean > 0.95)) #54
length(which(zbtb$mean < 0.05)) #1
zfrc<-frc3_mp_f[frc3_mp_f$linkage == 23,]
length(which(zfrc$mean > 0.95)) #54
length(which(zfrc$mean < 0.05)) #1
zpin<-pin3_mp_f[pin3_mp_f$linkage == "Z",]
length(which(zpin$mean > 0.95)) #54
length(which(zpin$mean < 0.05)) #1
zpsp<-psp3_mp_f[psp3_mp_f$linkage == 23,]
length(which(zpsp$mean > 0.95)) #54
length(which(zpsp$mean < 0.05)) #1
zrnv<-rnv3_mp_f[rnv3_mp_f$linkage == 23,]
length(which(zrnv$mean > 0.95)) #54
length(which(zrnv$mean < 0.05)) #1
pop idas mel Zidas Zmel
bcr 211 175 48 44
bic 82 8 57 0
bld 124 53 67 0
btb 137 15 67 4
frc 118 53 66 1
pin 91 9 56 2
psp 122 23 72 4
rdl 123 15 54 1
rnv 136 12 64 3