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]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.