Displaying 5 results from an estimated 5 matches for "homozygous".
2009 Nov 04
3
Cannot Change Function (PR#14041)
...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 :...
2006 Dec 31
1
Genotype importing from Sequenom
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 ac...
2009 Jan 22
1
infer haplotypes phasing trios tdthap
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 other...
2006 Dec 07
6
Response To Form Submission Hanging
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
2011 Oct 08
1
HWEBayes, swapping the homozygotes genotype frequencies
I evaluated the Bayes factor in the k=2 allele case with a "triangular"
prior under the null as in the example in the help file:
HWETriangBF2(nvec=c(88,10,2))
[1] 0.4580336
When I swap the n11 entry and n22 entry of nvec, I received totally
different Bayes factor:
>
> HWETriangBF2(nvec=c(2,10,88))
[1] 5.710153
>
In my understanding, defining the genotype frequency as