Hi all,

I want to construct relatedness among individuals and have a look at the
following script.

#########################
rm(list=ls())

N=5
id   = c(1:N)
dad = c(0,0,0,3,3)
mom  = c(0,0,2,1,1)
sex  = c(2,2,1,2,2) # 1= M and 2=F

   A=diag(nrow = N)
   for(i in 1:N)    {
      for(j in i:N)      {
         ss = dad[j]
         dd = mom[j]
         sx = sex[j]
          if( ss > 0 && dd > 0 )
            {
              if(i == j)
                   { A[i,i] = 1 + 0.5*A[ss,dd] }
                 else
                  { A[i,j] = A[i,ss] + 0.5*(A[i,dd])
                    A[j,i] = A[i,j] }
            }

      } #inner for loop
     } # outer for loop
  A

If the sex is male ( sex=1)  then I want to set A[i,i]=0.5*A[ss,dd]
If it is female ( sex=2) then A[i,i] = 1 + 0.5*A[ss,dd]


How do I do it ?

I tried several cases but it did not work from me. Your assistance is
highly  appreciated  in advance

Thanks

	[[alternative HTML version deleted]]

