similar to: merge on criteria

Hi, i have two files (file1.txt and file2.txt) which i would like to merge, based on certain criteria, i.e. it combines data based on matching geneID and exons. i have used the merge option, but it does not give me the desired outcome. merged.txt shows the result i would like. *File1. txt* ** AffyProbe ProbeType Flag GeneSymbol GeneID Exons Chrom Strand Affytart AffyEnd 1

Hi, I have the following function: > kurtosis <-function(x) (mean((x-mean(x))^4))/(sd(x)^4) #x is a vector and data > print(mydata) V1 V2 V3 V4 V5 1 1007_s_at DDR1 2865.1 2901.3 1978.3 2 1053_at RFC2 103.6 81.6 108.0 3

hell all: I have a vector as follows: > head(res) 1007_s_at.value 1053_at.value 117_at.value 121_at.value 1255_g_at.value 0.225801033 0.009747404 0.709517954 0.008825740 0.496859178 1294_at.value 0.005091231 after I convert the res into a data frame I got the following: resx<- data.frame(res) > head(resx) res

Hi, I have the following data file to be parsed and captured as a data frame: __DATA__ #GDS_ID GENE_NAME GENE_DESCRIPTION GENE_FUNCTION 1007_s_at | DDR1 | discoidin domain receptor tyrosine kinase 1 | protein-coding 1053_at | RFC2 | replication factor C (activator 1) 2, 40kDa | protein-coding 117_at | HSPA6 | heat shock 70kDa protein 6 (HSP70B') | protein-coding __END__ In particular it is

hi all i started recently using R and i found myself stuck when i try to analyze microarray data. i use the "affy" package to obtain the intensities of the probes, i have two CTRs and two treated. HG.U133A.Experiment1.CEL HG.U133A.Experiment2.CEL HG.U133A_Control1.CEL HG.U133A_Control2.CEL 1007_s_at 2156.23115 467.75615 364.60615 362.11865

Hi, I have the following data frame: > print(mydatframe) __DATAFRAME__ V1 V2 V3 1 1007_s_at DDR1 discoidin domain receptor tyrosine kinase 1 2 1053_at RFC2 replication factor C (activator 1) 2, 40kDa 3 117_at HSPA6 heat shock 70kDa protein 6 (HSP70B') __END__ Is there a way to create a hash with V2 as Key and V3 as its value? - Gundala Viswanath Jakarta - Indonesia

Dear mailing list, I have a C function that gives me a wrong result when I run it under Windows 7. This is the code under Linux (RHEL5): > library(phenoTest) > data(epheno) > sign <- sample(featureNames(epheno))[1:20] > score <- getFc(epheno)[,1] > head(score) 1007_s_at 1053_at 117_at 121_at 1255_g_at 1294_at -1.183019 1.113544 1.186186 -1.034779 -1.044456

Hi Everyone, I am very new to R. When i run the code yesterday, it was working fine. I was able to find gene annotation. somehow, today, when i try to run it again, it is giving me errors message that i don't understand. It says that it cannot open file. what file is it looking for? ###################################################### ##make sure to install and load annotationtool in R

Dear all, Here is my code which am using to combine 5th column from different data sets. Here is the function to do my job genesymbol.append.file <-NULL gene.column <- NULL readGeneSymbol <- function(files,genesymbol.column=5){ for(i in fnames){ temp <- read.table(fnames,header=T,sep="\t",stringsAsFactors=F,quote="\"")

Hi all, My name is Amy, I am a masters student in Bioinformatics at North Carolina State University. I am working on a project and I am trying to use the lumi R package for microarray data analysis. I have shown the sample code here and have questions about modifying the sample code for my own data. lumi package in R, example.lumi, the sample data has 8000 features and 4 samples I have

Hi all, My name is Amy, I am a masters student in Bioinformatics at North Carolina State University. I am working on a project and I am trying to use the lumi R package for microarray data analysis. I have shown the sample code here and have questions about modifying the sample code for my own data. lumi package in R, example.lumi, the sample data has 8000 features and 4 samples I have

require(gplots) data(rtPCR) overplot( RQ ~ Conc..ug.ml. | Test.Substance, data=rtPCR, subset=Detector=="ProbeType 7" & Conc..ug.ml. > 0, same.scale=TRUE, log="xy", f=3/4, main="Detector=ProbeType 7", xlab="Concentration (ug/ml)", ylab="Relative Gene Quantification" ) # Error in lowess.default(mf[[-response]], mf[[response]], f = f,

R-users, I have the following piece of code which I am trying to run on a dataframe (aga2) with about a half million records.  While the code works, it is extremely slow.  I've read some of the help archives indicating that I should allocate space to the p1 and ags1 vectors, which I have done, but this doesn't seem to improve speed much.  Would anyone be able to provide me with advice on

Hi all, I originally posted this to the bioconductor group, but maybe it's better suited to the r-help... I'm using som() to partition samples of gene expression data into clusters. The point is to classify control vs. experimental cases (sample clustering). The original matrix was 22283 x 8. The 8 samples have 4 controls and 4 experimentals. I transposed the matrix so that its dim

Hi, everyone When I ran this cript, There is Error in substring(tmp.subject, tmp.end[ex] + 1, tmp.start[ex + 1] - 1) : invalid substring argument(s) Could someone figure out what the problem is? for(i in 1:length(genebody[,1])){ tmp.id<-as.vector(genebody[i,1]) # get gene id tmp.subject<-as.vector(genebody[i,2]) # get gene sequence

Hi list, I have a question about using *read.table()* to read in a txt file. Basically, it consists of 16346 rows, 6 columns (no header). The code I used is: exprSet <- read.table('process_all4_GSA2.txt', row.names = 1,header =FALSE) and I got an error message: > exprSet <- read.table('process_all4_GSA2.txt', row.names = 1,header =FALSE) Error in

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

Hi all, Suppose I have a vector like this: [1] "STAT1" "STAT1" "STAT1" "STAT1" "GAPDH" "GAPDH" "GAPDH" "ACTB" "ACTB" [10] "ACTB" "DDR1" "RFC2" "HSPA6" "PAX8" "GUCA1A" "UBE1L" "THRA" "PTPN21" [19]

Hi, I was wondering I'm going about this in the correct way. I need to test if there are coding sequences or exons in hg19 which match a string of 100bp "D" i.e. [A,G or T]. However I'm getting a strange result. I get a hit on chr7, using the 100bp search however when I search with 60bp sequence of "D" I don't get any hits. library("BSgenome")

HI, this is my problem I want to subset this file df, using only unique df$exon printing the line once even if df$exon appear several times: unique(df$exon) will show me the unique exons If I try to print only the unique exon lines with df[unique(df$exon),] -this doesn't print only the unique ones :( could you help? thanks Nat exon size chr start