similar to: get compressed data via a socket connection

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.

Dear useRs, the seqinR package contains utilities to import and analyze biological sequence data. For a general introduction see this document: http://pbil.univ-lyon1.fr/software/SeqinR//vignette.pdf Please do not use r-help for questions about seqinR or r-bugs for bug report about seqinR. Use instead the seqinR diffusion list: http://pbil.univ-lyon1.fr/software/SeqinR//mailing.php?lang=eng A

Dear useRs, the seqinR package contains utilities to import and analyze biological sequence data. For a general introduction see this document: http://pbil.univ-lyon1.fr/software/SeqinR//vignette.pdf Please do not use r-help for questions about seqinR or r-bugs for bug report about seqinR. Use instead the seqinR diffusion list: http://pbil.univ-lyon1.fr/software/SeqinR//mailing.php?lang=eng A

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

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

Good afternoon. I am a master student in University of Porto in Portugal. At this moment I’m starting to use R, so I have some doubts. The aim of my analysis is: calculate a pairwise FST matrix from fasta file and creat a principal component analyses with adegenet package (I use seqinr and ape package to read this file, then I convert this file into a genind object with DNA2genind function

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

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

Hello R developers, I am working on the "seqinr" package and I encounter a tricky problem using a C function. We defined a C fonction called "getzlibsock" which is dedicated to compressed socket connections. This function is using the R internal C function called "getConnection(int)" in order to get information about the socket previously opened with the

Dear all, I have a data which is in text file and i would like to import the data to R. From the manual, i?ve found the read.table command function is the most appropriate but when i wrote the command an error had occur. It say ?Error in read.table"C:/Users/user/Documents/cfa-1.txt", header = T, sep = "\t",skip=10) :more columns than column names?. Please help me with this as

Hi !  Facing problem with " getSequence" commend .  when only biomaRt package loaded the following example working well  >mart <- useMart("ensembl",dataset="hsapiens_gene_ensembl") >seq = getSequence(id="BRCA1", type="hgnc_symbol", seqType="peptide", mart = mart) show(seq) but when i have loaded the seqinr, i got problem

Dear useR's I have a small basic problem which I am hoping to get some help with. I have a data frame, testSeq_df, with 1 row and 500 columns. Each column is a character (a,c,g or t). I want this sequence to have 4 factors (a,c,g,t). When I try the following: for(i in 1:500){ if (length(levels(testSeq_df[,i]))==1) levels(testSeq_df[,i]) <-

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 would like to use 'write.fasta(sequences, names, nbchar = 60, file.out, open = "w")' to convert a DNA sequence in a text file to fasta format. How do I read the the text file to prepare the argument 'sequences' of the function. The DNA sequence in the text file is one line as below: ATCACACAACGACACTCACCCTGGACGCTCATC......... Thank you [[alternative HTML

I am trying to return an index for a data set by searching using filenames. The name may be ANG_AUT.N.0734C70411A-1_1sA_0734C70411A.fasta, but i'd just like to search it using the term "0734C70411" as the file may be 0734C70411A or 0734C70411C or 0734C70411D Any way to do this other than doing something like this. where 0734C70411A is part of matrix list[,8]