dist.dna in ape package proportion of invariant sites

Hi,

I know there is a way to include a gamma (G) parameter in the dist.dna() function in R [gamma = TRUE], but is there a way to include an proportion of invariant sites (I) parameter?

Thank you

r dist.dna ape

• 1,932 views • link

updated 9.5 years ago by Biostar

• written 9.6 years ago by matthews.alix •

0 answers

No answers yet.

More posts like this

DESeq2: Variable Importance
written 5.2 years ago by gt •

Is there a way to obtain the relative importance of a variable in the DESeq2 model? I know you can run a likelihood ratio test …
mean evolutionary rate (MER) calculations
written 8.0 years ago by sansritisinha •

I used MEGA5.0 for computing MER for my dataset. the rates were computed based on JTT with an account of invariant sites model using 2 …
Pairwise distance matrix with nucleotide substitution model
written 8.4 years ago by luca.zoccarato.88 •

Dear community, I have a sequence alignment (full 16S rDNA) and I'd like to calculate the sequence pairwise distance matrix using a nucleotide substitution model. …
question about MrBayes software: evolutionary models
written 9.4 years ago by gobiidae.marin •

I'm a PhD student and read your tutorial about MrBayes softwar. for my thesis, I working with MrBayes software to analyzing data. your help and …
Equation for Log-normal rank-abundance distribution
written 10.0 years ago by songzewei •

I was checking the radfit {vegan} package in R, and here is what they have for the Log-normal rank abundance distribution: Function rad.lognormal fits a …
getting chromosome lengths from a vcf file
written 10.2 years ago by outlier95 •

I have a VCF file that includes variant *and* invariant sites for every locus. Is there an easy way to obtain the lengths for each …
dist.alignment vs dist.dna for DNA distance calculation
written 10.4 years ago by rimjhim.roy.ch

I am trying to get the proportion of sites that differ between each pair of sequences in my alignment. One approach would be: dist.dna from …
Pairwise distances DNA sequences (K80 model of DNA evolution) [Python]
written 10.7 years ago by erica.fary •

Hi everyone, I need to build phylogenetic tree based on aligned DNA sequences. I used first from Bio.Phylo the following functions: ``` distM = DistanceCalculator('ident').get_distance(my_aligned_sequences) …
qvalue package in R: Can it be used with very small proportion of true alternatives?
written 10.9 years ago by ckhh •

Hi, I'm trying to use 'qvalue' package in r to work out the true proportion of alternatives. I get pi0 is >1 with small lambda. …
How To Get Non-Snp Sites Or Invariant Sites Based On 1000 Genomes Project
written 13.1 years ago by Gangcai

<p>Hi everyone, I need to find the human genomic sites which are invariant among all the individuals (identified by sequencing in all individuals and also …