m_qtlmap_haplotype_V1

[ Top ] [ Modules ]

NAME

    m_qtlmap_haplotype_V1

DESCRIPTION

NOTES

SEE ALSO

haplotype_V1

[ Top ] [ m_qtlmap_haplotype_V1 ] [ Subroutines ]

NAME

    haplotype_V1

DESCRIPTION

NOTES

SOURCE

 53    subroutine haplotype_V1()
 54         integer           :: stat,ll,c
 55         integer ,dimension(nchr)       :: valnpo
 56         if (maxval(nmk) > MAX_MARKER) then
 57           call stop_application("You can not used this version of haplotype calculation.")
 58         end if
 59 
 60         ! BUFFER STRUCTURES...
 61         !-------------------------
 62         ll = 2**maxval(nmk)
 63 
 64         allocate (h(nchr,maxval(nmk),ll),STAT=stat)
 65         call check_allocate(stat,'h [m_qtlmap_haplotype_V1]')
 66         allocate (prohp(nchr,ll,np),STAT=stat)
 67         call check_allocate(stat,'prohp [m_qtlmap_haplotype_V1]')
 68         allocate (prohm(nchr,ll,nm),STAT=stat)
 69         call check_allocate(stat,'prohm [m_qtlmap_haplotype_V1]')
 70         allocate (prot(nchr,maxval(nmk),nd,4),STAT=stat)
 71         call check_allocate(stat,'prot [m_qtlmap_haplotype_V1]')
 72         allocate (ptfin(nchr,maxval(nmk),(nd*4),4),STAT=stat)
 73         call check_allocate(stat,'ptfin [m_qtlmap_haplotype_V1]')
 74         !HAPLOTYPE_DATA --> DONNEE PERSITENTE A L APPLI
 75         allocate (genotyp(nchr,maxval(nmk),(size(numero)),2),STAT=stat)
 76         call check_allocate(stat,'genotyp [m_qtlmap_haplotype_V2]')
 77         genotyp=nmanque
 78         allocate (ngenom(nchr,nm+1),STAT=stat)
 79         call check_allocate(stat,'ngenom [m_qtlmap_haplotype_V2]')
 80         allocate (phasp(nchr,np),STAT=stat)
 81         call check_allocate(stat,'phasp [m_qtlmap_haplotype_V2]')
 82         allocate (phasm(nchr,nm),STAT=stat)
 83         call check_allocate(stat,'phasm [m_qtlmap_haplotype_V2]')
 84         phasp = .false.
 85         phasm = .false.
 86 
 87         call combine()
 88         call ancetre()
 89         call gammapf(prot)
 90         call pdegp()
 91         call pdegm()
 92 
 93         call log_mess('Second dim of pdd:'//str(ngend(1,size(ngend,2))),DEBUG_DEF)
 94         do c=1,nchr
 95           valnpo(c)=get_npo(c)
 96         end do
 97 
 98         allocate( pdd(nchr,maxval(ngend),4,maxval(valnpo)) )
 99         pdd=0.d0
100 
101 
102         call pded()
103 
104         deallocate(prot)
105         deallocate(ptfin)
106         deallocate(prohp)
107         deallocate(prohm)
108         deallocate(h)
109 
110    end subroutine