search for: homozyg

...mp;&((Re(Tot_log_sum))>log(5/4)/log(2))&&((Re(Tot_log_sum))<log(7/4)/log(2))) {return(3)} if(((Im(Tot_log_sum))==0)&&((Re(Tot_log_sum))>log(7/4)/log(2))) {return(4)} } When I try to run it in the loop: for (j in 1:length(BB_mean_ref)){ # find copy number states (0 for homozygous deletion, 1 for hemizygous deletion, 2 for neutral, 3 for single amplification, 4 for multiple amplification) state_log_sum[j]=CN_state_log_sum(Tot_log_sum[j]) state_sum_log[j]=CN_state_sum_log(Tot_sum_log[j]) } I get the error message: Error in Im(Tot_sum_log) != 0 || Re(Tot_sum_log) <- 2...

Sequenom has an odd format of calling a SNP genotype gg [1] "C" "GA" "A" "C" "C" "AG" "C" "C" "T" "G" homozygous A is called A and heterozygous is called AT The genetics package cannot handle the fact that some genotypes are declared with 2 letter while other are declared with only 1. Consequently the genotype() or makeGenotypes() will not work. I need to either find a clever way that the genetics package...

Dear R mailing list, I have a dataset with genotypes from trios and I would like to infer haplotypes for each mother, father and child. The package that I could find that can do this is tdthap. But when the mother is homozygous (e.g., 2/2) the haplotype is called as not possible to infer (0); I would prefer for it to call the genotype (2). From what I understand it is doing what I would like for the father (example below). Can anyone provide me with some information about this tdthap behaviour? And is there any ot...

Hello, I am using Mechanize to post a form to a website. When I do this by hand in my browser the response takes about 35s to come back (it''s a long page full of tables and graphics). When I do this with Mechanize, the server starts to respond and then appears to hang. The obvious conclusion is that my code is wrong but I am reasonably sure that I haven''t altered it

...)) [1] 1.542047 >DirichNormHWE(nvec=c(2,10,88), c(1,1))/DirichNormSat(nvec=c(2,10,88), c(1,1,1)) [1] 1.542047 Could you explain why the HWETriangBF2 is returining completely different values of Bayes Factor?? -- View this message in context: http://r.789695.n4.nabble.com/HWEBayes-swapping-the-homozygotes-genotype-frequencies-tp3886313p3886313.html Sent from the R help mailing list archive at Nabble.com.