Hello all: I have a protein sequence alignment in a data frame (align1, 72 x 236), where each row is a protein and each column a site in the alignment. AA is vector of amino acid symbols plus "-" (gap). I can calculate amino acid frequencies at each site by:
>align1.F <- matrix(0,nrow=22,ncol=236,dimnames=list(AA,seq(1:236))) >for(i in 1:236) > align1.F[names(summary(align1[[i]])),i] <- (summary(align1[[i]])/length(align1[[i]])) Is there a more efficient (i.e., without a loop) way to do this? Is there some way to use table or ftable to create an 22 x 236 table of amino acid frequencies from align1 and AA in one fell swoop? Thanks, Kurt Wollenberg, PhD Tufts Center for Vision Research New England Medical Center 750 Washington St, Box 450 Boston, MA, USA [EMAIL PROTECTED] 617-636-8945 (Fax) 617-636-9028 (Lab) The most exciting phrase to hear in science, the one that heralds new discoveries, is not "Eureka!" (I found it!) but "That's funny ..." --Isaac Asimov ********************** Confidentiality Notice\ **********************\ The inf...{{dropped}} ______________________________________________ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html