similar to: list genes w/n a genomic fragment

Hi All, I have an CBS segmentation algorithm output for 10 tumor samples each from 2 different tumors. Now, I am in an urgent need to assign gene (followed by all genes present) that belong to a particular segment after I removed all the CNVs from segment data. The format of the data is: Sample Chromosome Start End Num_Probes Segment_Mean Sample1A-TA 1 51598 76187 15

Hi, I am given the following table: > head(hsa_refseq) chr genome region start stop nu strand nu.1 nu.2 gene_id 1 chr1 hg19_refGene CDS 67000042 67000051 0 + 0 gene_id NM_032291 2 chr1 hg19_refGene exon 66999825 67000051 0 + . gene_id NM_032291 3 chr1 hg19_refGene CDS 67091530 67091593 0 + 2 gene_id NM_032291 4 chr1 hg19_refGene exon

Dear useRs, I have a slightly complicated data structure and am stuck trying to extract what I need. I'm pasting an example of this data below. In some cases, there are duplicates in the "gene_id" column because there are two different "sample 1" values for a given "sample 2" value. Where these duplicates exist, I need to average the corresponding

> #sort the data by predicted probability > b.order<-bo.id.pred[(order(-predict)),] > b.order[1:20,] gene_id predict 43 637882902 0.07823997 53 638101634 0.66256490 61 639084581 0.08587504 41 637832824 0.02461066 25 637261662 0.11613879 22 637240022 0.06350477 62 639084582 0.02238538 63 639097718 0.06792841 44 637943079 0.04532625 80 640158389 0.06582658 3 637006517 0.57648451

Hi, r-community, This morning, I MET the following problem several times when I try to attach the data set. When I closed the current console and reopen the R console, the problem disappear. BUt with the time passed on, the problem occurs again. Can anyone help me with this? > attach(total) The following object(s) are masked from total ( position 3 ) : acid base cell_evalue

> head(en.id.pr) valid.gene_id b.pred rf.pred svm.pred 1521 2500151211 0 0 0 366 639679745 0 0 0 1965 2502081603 1 1 1 1420 644148030 1 1 1 1565 2500626489 1 1 1 1816 2501711016 1 1 1 > p.pred <- data.frame(en.id.pr, sum=apply(en.id.pr[,2:4], 1, sum)) #

Hi everyone, I am trying to find a way to filter a table; If I am given for example the following table: > head(intra) chr miRNA start end strand ACC hsa_ID region region_start region_end gene_id transcrip_id 1 chr1 miRNA 1102484 1102578 + ACC="MI0000342"; ID="hsa-mir-200b"; exon 1102484 1102578 NR_029639 NR_029639 2 chr1

for the script, please kindly see the script below. At line 10 and line 13, my problems occurs. The first one is I try to retrieve the gene official name from a column of a table. The pattern of official name is something starting with gene_name. For detail problems, please see the according lines. Any suggestions are appreciated example of matching source (extract the Nnat, sometime it would

Dear All, I need to fetch GO ontologies for Homo sapiens with their mappings to corresponding Uniprot identifiers. I would be using this information to compare result from a clustering algorithm with existing protein complexes. This would be a test to check how the clustering algorithm accurately captures GO terms with respect to the known protein complexes. Can anyone suggest a simple workflow

Dear useRs, I am glad to announce that the new R package "hypred", initial version 0.1, is now available on CRAN. "hypred" is a package for simulating high-density SNP data. Its main function, "hypredRecombine" is intended to be used as a "Software tool" in larger programs that simulate complex populations. The focus of the package is on producing

Dear useRs, I am glad to announce that the new R package "hypred", initial version 0.1, is now available on CRAN. "hypred" is a package for simulating high-density SNP data. Its main function, "hypredRecombine" is intended to be used as a "Software tool" in larger programs that simulate complex populations. The focus of the package is on producing

Dear Frank, I hope everything is well with you. ? ?First of all? ?, I? ?'d like thank you for your nice genomic simulation program (*hypred *package) and also for taking your time to answer my question. Dear Frank, I'm a PhD student of animal breeding science from Iran. Unfortunately although I try more time for running your package in order to simulate of a historical and reference

Great R people, I have noticed a strange behaviour in read.delim() and friends in the R 2.7.0 version. I will describe you the problem and also the solution I already found, just to be sure it is an expected behaviour and also to tell people, who may experience the same difficulty, a way to overcome it. And also to see if it is a proper behaviour or maybe a correction is needed. Here is the

I am running a server with Centos 4.4, using a combo of sendmail-clamav-spamassassin for mail, which has been working quite well up until last night. Around 3:00am, all mail going to all users was being rejected by spamassassin. My ~/mail/.procmailrc is as follows: ======================== MAILDIR=$HOME/mail # :0 H * ^X-Spam-Status:.*Yes { EXITCODE=67 :0: /dev/null } =========================

Hi there, I have a problem when I try to create a new user with swat. I use binaries of Samba 2.2.2 coming from www.samba.org, or from www.sunfreeware.com on Solaris 2.6 and on 8. I compiled the sources and the result is the same. When I click on "Add New User" (Server Password Management), I receive a nearly blank page: I can only see the samba gif on top and the user is not

I have an rsync cron script that copies one entire hard drive containing hard links to another drive. It does this by selecting the directories to be rsync'ed and copying them one by one. My rsync line is: rsync -uav source/ target and the problem is this would not copy hard links. If I change the rsync line to: rsync uaHv --delete source/ target will the files copied using the original

Genome Institute of Singapore (GIS) Biostatistics Group The Genome Institute of Singapore ( http://www.gis.a-star.edu.sg <http://www.gis.a-star.edu.sg> ) is looking for statisticians who are interested in a wide range of biomedical probems. We are especially interested in hearing from senior researchers interested in leading their own group. The institute is well-funded and engaged

Por si es de interés: Data Analysis for Genomics Data Analysis for Genomics will teach students how to harness the wealth of genomics data arising from new technologies, such as microarrays and next generation sequencing, in order to answer biological questions, both for basic cell biology and clinical applications. https://www.edx.org/course/harvardx/harvardx-ph525x-data-analysis-genomics-1401

We are seeking a talented and driven postdoctoral fellow to join the laboratory of Jake Michaelson, PhD, assistant professor in the department of psychiatry at the University of Iowa. The Michaelson lab investigates the effect of genetic variation on the development and function of the brain, particularly in the context of neurodevelopmental conditions such as autism and language impairment. Our

Hello, I am using R to look at whole-genome gene expression data. This means about 27,000 genes, each with a vector of numbers reflecting expression at different tissues and times. I need to do an all against all co-expression calculation (basically, just calculate Pearson's r for every gene-gene pair). I try to store the result of such a thing in a 27000x27000 matrix, but r seems not to like