More posts like this
-
Convert R script to python(implementation of the UPGMA OR neighbor joining algorithm and test succe…
written by anasjamshedWe have MAF file for alignment which need to be converted. The MAF file was converted into .fasta format for handling with ease in the …
-
How to find differential gene expressions
written by Raz •Hi, I am a mathematician and new to mRNA analysis. I am analyzing some mRNA seq data. My data set contains four-time points. in each …
-
Difference between "LogNormalize" normalization method in Seurat and z score normalization
written by f.morgana •Hi all, Perhaps a silly question but the default normalization method in Seurat is "LogNormalize". Instead of automatically running the code, I would like to …
-
What is the difference between deferentially expressed genes and deferentially expressed transcrip…
written by WUSCHELI know the Biology behind `deferentially expressed transcripts`. But I do not understand **what does it mean / difference between** `deferentially expressed genes` and `deferentially …
-
Distance based tree reconstruction
written by Jess •I want to generate a tree using manhattan distance matrix (normalized to 0-1). This is purely genetic diversity study and no evolutionary changes are involved. …
-
Zygosity of variants/allele in VCF
written by genie66 •I am new to bioinformatics and trying to understand the zygosity of variant/allele in VCF data. Couldn't find proper bioinformatics explanation anywhere, can someone please …
-
Count vs Read
written by mhyunjunkang •Hi everyone, I have a very basic question, but I still don't understand. What is difference between "count" and "read" in RNA-seq data? If I …
-
SSPACE option (-k and -a) for scaffolding
written by PicasaHi, I am trying to play with SSPACE's parameters, especially the -k and -a option. According to the manual: -k Minimum number of links (read …
-
dist.dna, dist.ml and phylogenetic tree
written by l.roca •Hi, I wish to prepare a a phylogenetic tree based on an aligned fasta file. Is the script below correct? ``` data <- read.dna("F:/AlignedSeqs.fasta", format …
-
Error Model For Sequence Simulation
written by dli<p>Is there a simple example for the error mode while simulating next generation sequencing reads? i know that wgsim in samtools and many other tools …
Wikipedia has a pretty straightforward explanation/illustration of NJ -- http://en.wikipedia.org/wiki/Neighbor_joining. Can you make your question more specific by saying what you are still missing having read this?