Displaying 1 result from an estimated 1 matches for "sitesid".
Did you mean:
siteid
2007 Jan 21
2
efficient code. how to reduce running time?
...all(test) || all(!test) )
{
print (paste("perfect correlation with", sites[i]));
print (paste("removing alleles", sites[j]));
newsites <- newsites[newsites!=sites[j]];
}
}
}
sites <- newsites;
print(sites);
# QTN: the site nearest right to S/4
sitesid = floor(length(sites)/4);
QTNid = sites[sitesid];
QTN = snp[,QTNid];
print (paste("QTN:", names(snp[QTNid])));
print (QTN);
# remove QTN from sites
sites <- sites [ sites != QTNid ];
print(sites);
print (paste("Number of usable SNPs:", length(sites)));
# p: allele frequenc...