I'm learning functional programming and wonder if someone could help me convert some R code into functional F#. I have used matrices for storage and for loops in R, but realize that there must be better options in F#. The only keypoint with the code below is that A and D should be matrices. I also include some toy data.

R code:

pedtab <- matrix(scan("Pedfile.txt"),ncol=3,byrow=TRUE)

nind <- nrow(pedtab)

sire <- pedtab[,2]

dam <- pedtab[,3]

Aprim <- matrix(0,ncol=nind+1,nrow=nind+1)

Dprim <- matrix(0,ncol=nind+1,nrow=nind+1)

for(c in 2:(nind+1)) {

for(d in 2:c) {

if (c == d) {

Aprim[c,d] <- 1+0.5*Aprim[sire[c-1]+1, dam[c-1]+1]

Aprim[d,c] <- Aprim[c,d]

diag(Dprim) <- 1

}

else {

Aprim[c,d] <- 0.5*Aprim[sire[c-1]+1, d]+0.5*Aprim[dam[c-1]+1, d]

Aprim[d,c] <- Aprim[c,d]

Dprim[c,d] <- ((Aprim[sire[c-1]+1, sire[d-1]+1]*Aprim[dam[c-1]+1, dam[d-1]+1])+

(Aprim[sire[c-1]+1, dam[d-1]+1]*Aprim[dam[c-1]+1, sire[d-1]+1]))/4

}

}

}

A <- Aprim[2:(nind+1),2:(nind+1)]

D <- Dprim[2:(nind+1),2:(nind+1)]

"Pedfile.txt"

1 0 0

2 0 0

3 0 0

4 0 0

5 0 0

6 0 0

7 0 0

8 0 0

9 0 0

10 1 4

11 1 4

12 1 4

13 1 5

14 1 5

15 1 5

16 2 6

17 2 6

18 2 6

19 2 7

20 2 7

21 2 7

22 3 8

23 3 8

24 3 8

25 3 9

26 3 9

27 3 9