similar to: R

Dear R developers, I am having a slightly weird issue with a S4 method defined in my package adegenet 1.2-4, with R 2.11.0. As far as I know, the problem is new, and the code implementing the method has not changed for more than a year and worked well so far. The problem is the following. I define, in the package, a method "[" for the S4 class 'genind'. The method's

Hello, I am working with a matrix of multilocus genotypes for ~180 individual snail samples, with substantial missing data. I am trying to calculate the pairwise genetic distance between individuals using the stats package 'dist' function, using euclidean distance. I took a subset of this dataset (3 samples x 3 loci) to test how euclidean distance is calculated: 3x3 subset used

Hi all, I am learning R by "doing". And this is my first post. I want to use R: 1- to fetch a DNA sequence from a databank (see bellow) and 2- store it as FASTA file. The problem: neither an error is prompted nor the fasta file is created. Testing the code (see bellow), I notice that everything works until the *"write.dna" *command - which is not creating the fasta file.

Dear R developers I am currently working on the seqinR package. The seqinR package allows a remote access to biological databases via a socket connection. We are using the functions socketConnection, writeLines and readLines to open the socket, send request to the server and receive response from the server respectively. Recently, a new function implemented in the socket server allows

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..

Dear R-helpers: thank your for your attention. i am a newer to R and i am doing some protein category classification based on the amino acid sequence.while i have some questions urgently. 1. any packages for analysis amino acid sequence 2. given two sequences "AAA" and "BBB",how can i combine them into "AAABBB" 3. based on "AAABBB",how can i get some

Hi - When I'm trying to read in a text file into a labeled character array, the memory stamp/footprint of R will exceed 4 gigs or more. I've seen this behavior on Mac OS X, Linux for AMD_64 and X86_64., and the R versions are 2.4, 2.4 and 2.2, respectively. So, it would seem that this is platform and R version independant. The file that I'm reading contains the upstream regions

I suspect you meant WD <- "~/Documents/Scripting/R_Studio/Sequences/" but I am entirely unfamiliar with the packages you are using, and know nothing about what is on your hard drive. For future reference: A) Read the Posting Guide. This is a plain text email list, and your html formatting gets removed leaving a mess that is not always readable. B) Most frequent users of R

Hi all, I have a sequence file (fasta format) and want to calculate the rho statistics for dinucleotide abundance value on my data.. the code which I use is (using seqinr library and current working directory) seq_info<-read.fasta("gene.txt") rho(seq_info[1],2) but it yields only the dinucleotides, not their rho values, i.e, > rho(seq_info[1],2) aa ac ag at ca cc cg ct ga gc

Hi! I am trying to analyse my data using amova (http://www.oga-lab.net/RGM2/func.php?rd_id=pegas:amova): My input to R is a DNA sequence file, format=fasta dna<- read.dna("XX.fasta", format="fasta") #left other options as default d<- dist.dna(dna, model="raw") g<- read.table("XXX.design") Load necessary libraries: library(pegas)

Hi, I want to analyse DNA sequence data (mtDNA) in R as in calculate Fst, Heterozygosity and such summary statistics. Package Adagenet converts the DNA sequence into retaining only retaining the polymorphic sites and then calcuates Fst.. but is there any other way to do this? I mean analyse the DNA sequence as it is.. and calculate the statistics? Thanks! [[alternative HTML version deleted]]

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]]

Buenas, Me encuentro con el mismo problema, de que me dice que el argumento del if no es un "valor ausente donde TRUE/FALSE es necesario" Este es mi codigo de pruebas. readseq <- "aaaaaaaaaaa", "aaa", "aa") auxiliar <- count(readseq[j],i+2) aux_a <- auxiliar["listaa"] if(aux_a > 0){ matrizgraf3[i][k] = matrizgraf3[i][k] + 1 listaa

Estimado Gemma Ruiz Olalla De curioso miré una búsqueda en internet, y se la comparto, de esta observe la página 8 http://adegenet.r-forge.r-project.org/files/tutorial-basics.pdf Estoy de acuerdo con Jorge y Carlos, pero por las dudas, es solo una codificación, yo sin mirar más no interpretaría que el heterocigota tiene un "peso" de 0.5, porque si hay sobredominancia, dominancia, etc

Buenas tardes, Estamos intentando hacer un estudio sobre diversidad en gatos callejeros ('stray cats') de la base de datos 'nancycats' (librería 'adegenet') de R. Sin embargo, no entendemos muy bien los datos. ¿Alguien puede explicarnos a qué corresponde la variable respuesta? Entendemos que las filas son los 237 gatos (observaciones), y las columnas se corresponden con

Dear list, here is a problem I met when trying to use a method for an S4 object, without loading the package in which the method was defined. I do not know if this is a bug, or a mistake of mine. Normally, I think the package in which the appropriate method is defined is loaded automatically when calling the method. The problem is that if the package is indeed loaded automatically, the

Dear all, It is my understanding that using 'importFrom' in the NAMESPACE of a package allows to avoid conflicts between different packages defining objects with identical names. However, I can still see conflicts while loading the package using 'library'. Here is a toy example, with a package 'foo' importing 'as.igraph' from the igraph package, and 'nj'

Dear R users, seqinR 1.0-5 has been released yesterday on CRAN, so that the source code of the package should be available on all CRAN mirrors within the next 24h. The updated package vignette is here: http://pbil.univ-lyon1.fr/software/SeqinR/seqinr_1_0-5.pdf User level visible changes are: o A new function dotPlot() is now available.

Dear R users, seqinR 1.0-5 has been released yesterday on CRAN, so that the source code of the package should be available on all CRAN mirrors within the next 24h. The updated package vignette is here: http://pbil.univ-lyon1.fr/software/SeqinR/seqinr_1_0-5.pdf User level visible changes are: o A new function dotPlot() is now available.

I am trying to use the cca/rda/capscale functions in vegan to analyse genetic distance data ( provided as a dist object calculated using dist.genpop in package adegenet) with geographic distance partialled out ( provided as a distance object using dist function in veganthis method is attempting to follow the method used by Geffen et al 2004 as suggested by Legendre and . FORTIN (2010). I