Translation of a number-crunching program from Justbasic to R - repost

I have written a program in JustBasic (a free dialect of Basic) that inputs a comma delimited matrix of numbers, and uses a computationally-intensive procedure to compute a coefficient alpha, a statistic measuring the internal consistency of the numbers. I am looking for a programming freelancer who is very familiar with the statistical computing package R, who can rewrite my program in R. The original program is very highly commented, and I believe it should be easy to follow, even if you don't program in JustBasic. The procedure in the original code consists in taking the number of "items" in the "scale" (i.e. the number of variables that form the rows of the matrix) and randomly sorting them into two halves. Then the program computes the mean of each of these two subsets, for each case, treating missing values by basing the mean only on those values that are nonmissing. Then when done with all the cases, the program computes a Pearson correlation between the means of the two subsets. Then the program randomly divides the items into two subsets again, and does the same thing all over again, and continues for many random subsets. The program "steps up" the correlation coeeficients using the a variation on the Spearman-Brown formula. The program you write will be in the public domain or the GNU equivalent of the public domain, so your work will be "for hire" and you won't retain any copyright to it. I think that someone who is very familiar with writing software in R should be able to do this task quickly. The original code is not very long.