similar to: new package RFLPtools

Hello, I'm struggling with understanding the output on the ancestral.pars() command from the phangorn package, I'm new to doing phylogenetic analyses using R. I used it on nucleotide data, and it works fine, I'm just not sure how to read the output. The output is phyDat class, and outputs a matrix for each node/leaf in the tree. I figured out that the matrix columns represent the four

Dear all, I am working with taxonomic data, represented as a list of classes, orders, families, genera and finally species. > class(mydata) [1] "data.frame" > mode(mydata) [1] "list" > names(mydata) [1] "tclass" "torder" "tfamily" "tgenus" "tspecies" > length(mydata$tclass) [1] 161590 The first 10

Dear all: I have a character object x with ' (single-quote) character. x <- c('"hydrolase activity","actin binding","3',5'-cyclic-nucleotide phosphodiesterase activity") I want to write a function that will identify ' and replaces with \' myf <- function(term){ if (grep("'",term)) {

Job Title : Postdoctoral Position in Biostatistics/Bioinformatics Background The joint research unit between bioM?rieux, international leader in in vitro diagnostic, and Hospices Civils de Lyon, second-largest University Hospital Network in France, develops a research program dedicated to the study of injury induced immunosuppression. This program aims at (i) describing and understanding the

Dear useRs, The seqinR package is a library of utilities to retrieve and analyse biological sequences. A new version of seqinR, seqinR 1.0-7, has been released on CRAN. Here is a summary of changes: o A new *experimental* function extractseqs() to download sequences thru zlib compressed sockets from an ACNUC server is released. Preliminary tests suggest that working with about 100,000

Dear useRs, The seqinR package is a library of utilities to retrieve and analyse biological sequences. A new version of seqinR, seqinR 1.0-7, has been released on CRAN. Here is a summary of changes: o A new *experimental* function extractseqs() to download sequences thru zlib compressed sockets from an ACNUC server is released. Preliminary tests suggest that working with about 100,000

I have following codes for ggplots. The legends are given in the plot do not match with the values specified in the codes given below. Your helps highly appreciated. Greg library(ggplot2) p <- ggplot(a,aes(x=NO_BMI_FI_beta ,y=FI_beta ,color= Super.Pathway))+ theme_bw() +theme(panel.border=element_blank()) + geom_point(size=3) p2<-p+scale_color_manual(name="Super.Pathway",

Hello everyone! I have two questions regarding cluster analysis with the package 'cluster' : 1. I have a dissimilarity matrix in csv format, and I would like to use it to generate a dendogram using 'hclust'. How do I make that matrix compatible with the function? I think it has to be similar to the objects generated by the 'dist' function. Any ideas? 2. Does

Hi there, I am looking for R-packages that can help me visualize properties on nucleotide sequences. I want to display sequences in the 1-100K base range as lines and plot features above and below those lines. Any ideas would be welcome. Thanks, Bernd

Hi, I am very new to R and wanted to know if there is a package that, given very long nucleotide sequences, searches and identifies short (7-10nt) motifs.. I would like to look for enrichment of certain motifs in genomic sequences. I tried using MEME (not an R package, I know), but the online version only allows sequences up to MAX 60000 nucleotides, and that's too short for my needs..

Hi there, I am looking for R-packages that can help me visualize properties on nucleotide sequences. I want to display sequences in the 1-100K base range as lines and plot features above and below those lines. Any ideas would be welcome. Thanks, Bernd [[alternative HTML version deleted]]

I have a matrix of similarity scores that I want to convert into a matrix of dissimilarity scores so that I can apply some clustering methods to the data. That is, high values in my matrix signify similarity and low values (zero being the lowest) signify no similarity. What functions/options in R or its packages are available for making this kind of transformation of a matrix?

Ape is an R package for "analyses of phylogenetics and evolution". The first version (0.1) has been released on 27 August 2002 and is available on CRAN. >From the 'Description' file of version 0.1: Ape provides functions for reading, and plotting phylogenetic trees in parenthetic format (standard Newick format), analyses of comparative data in a

I think you mean between column 1 and 2 of A? Why is 36003918 not included? It is clearly between 35838396 and 36151202 in the first row of A. My earlier solution should work fine. Just create a new matrix AX that has the columns switched so that the start is always column 1 and use that to identify the ones you want to select. That way you are not modifying B. This will be faster than checking

Hello, I have a problem to treat my data. I seek the orders being able to treat under sampling: I have X samples divided into 10. How to take, in a random way, under sample from the 1st sample, and in addition, one under sample of the 2nd sample, and so on to X to calculate the average taxonomic richness of the selection. Then I would like to know how to renew the experiment by taking

Hi there, I am just starting in R and this might be a very basic question. I applied one on the examples of pairs() to my own data. The examples shows scatter plots on one side of the matrix and the correlation coefficients on the other which works well. I then modified it slightly because I want different colors in the plots. Now i am getting the error of unused arguments, because the color

I've a data frame with this structure: 'data.frame': 1987 obs. of 11 variables: $ site : Factor w/ 24 levels "B(W)","BC-1",..: 1 1 2 2 2 1 1 1 ... $ sampdate : Date, format: "2000-07-18" "2000-07-18" ... $ tclass : Factor w/ 8 levels "Annelida","Arachnida",..: 1 5 5 5 5 ... $ torder : Factor

Dear list, This question is about Hidden Markov Model. Given a transition matrix, an emission matrix and a sequence of observed symbols (actually, nucleotide sequences, A, T, C and G), I hope to predict the sequence of state by Viterbi algorithm. I searched R repository for related packages. msm package has function viterbi.msm (as well as very good document), but it only works for

All Is there a package or library that will, given a nucleotide sequence 1. calculate the extinction coefficient at 260 nm for (Beer-Lambert's law) 2. calculate molecular weight 3. return it's complementary sequence I was able to find several packages that can do similar calculations for an amino acid sequence for proteins but none for nucleic acids. Any pointers, etc. would be

Hello, I have two tab-delimited files which I would like to combine. In the first one I have gene IDs (Unique) on column 1 and than various experimental results from microarray analysis (see attached files list1 ) the second arrays have the same genes IDs (more and in a different order, some are double) (see attached files list2 ) What I would like to do is to search in the second list for gene