similar to: Converting EnSeMBL Probe names into Gene Name

Hello, I am looking for an explanation and/or fix for a discrepancy in the behaviour of the R lm.influence() function [ version R 1.5.0 (2002-04-29) ] and the same function in Splus [ Splus version 5.1 release 1, running on SGI IRIX 6.2]. The discrepancy is of concern because I am migrating some Splus scripts to R and need to ensure consistency of results. Specifically, when I fit a glm()

Hello all, I have a data file table.txt which i have attached. I am trying to pass the columns as arguments to a function "totnorm" where i am displaying a total normalization plot. The function is given below: totnorm<-function(x,y){scale<-sum(x)/sum(y);xlab<-colnames(x);ylab<-colnames(y);x1<-x[[1]];y1<-scale*y[[1]];plot(x1,y1,xlab=xlab,ylab=ylab,col=6, col.lab=4);}

Hello all, I have a data file table.txt which i have attached. I am trying to pass the columns as arguments to a function "totnorm" where i am displaying a total normalization plot. The function is given below: totnorm<-function(x,y){scale<-sum(x)/sum(y);xlab<-colnames(x);ylab<-colnames(y);x1<-x[[1]];y1<-scale*y[[1]];plot(x1,y1,xlab=xlab,ylab=ylab,col=6, col.lab=4);}

Dear all, i have the attached data called "new" as data.frame. First I have only three columns called Var1, Var2 and Freq and with bind I attached the column test for a color specification (TEST DATA below). With this plot function (require packages lattice) dotplot(reorder(Var1, rep(score, cl.count)) ~ Freq | Var2, data = DATA, origin = 0, type = c("p",

Hi everybody, I have a problem when trying to do the quality control with the packages simpleaffy and affyQCReport with the drosophila chip 2.0 At first I got the messeage, that the *.qcdef file is not there. I followed the instructions in tha manual and created the file like that: array drosophila2cdf alpha1 0.05 alpha2 0.065 spk bioB AFFX-r2-Ec-bioB-3_at spk bioC AFFX-r2-Ec-bioC-3_at spk bioD

Hello all, I am trying to perform ANOVA on my sample data given below to see if any gene(column 1 stands for gene names) is differentially expressed after subjecting it to the 6 different experiments(columns 2 to 7 are experiments). Gene 14A_U133A_Detection 14B_U133A_Signal 88A_U133A_Signal 88B_U133A_Signal 183A_U133A_Signal 183B_U133A_Signal AFFX-BioB-5_at 403 409.3 611.5

I have a data frame which has the following data. data<-read.table("table.txt",header=TRUE) data X14A_U133A_StatPairs X14A_U133A_Detection X14B_U133A_Signal 1 AFFX-BioB-5_at 403.0 409.3 2 AFFX-BioB-M_at 757.3 574.4 3 AFFX-BioB-3_at 284.4 327.3 4 AFFX-BioC-5_at

Dear all, i try to plot with ggplot2. Therefor I have an matrix with 3 colums. With cbind I add an additional column called "col". I need this column "col" because in a later step and want to specify here some plot details which I will get from another analysis If I want to plot with this code, I have the problem that the legend is wrong. Blue changed to green and green to

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

Hello, I'm using recursive SVM script (rSVM - http://www.stanford.edu/group/wonglab/RSVMpage/R-SVM.html ) on some microarray data. The data to be input are log2, as numeric matrix w/ attributes -- str(svm_num_mat) num [1:10, 1:12340] 13.1 13.1 13.1 13.1 13.0 ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:10] "rma_log2_con_sample_1"

Hi, Reading help("Documentation"), I'm led to believe that a help call like: ?myFun(x, sqrt(wt)) Will search for help on the appropriate method in the case that myFun is generic. This isn't working for me. Here is an example using the Biobase package: ## If Biobase is not installed source("http://bioconductor.org/biocLite.R") biocLite("Biobase")

Hi, I use write.table() to write a file to an external xls file. the column names left-shift one position in output file. I check with col.names() row.names(), the file is fine. How to prevent the shifting? I71 I111 I304 I307 I305 I306 I114 I72 AFFX-BioB-5_at 6.66435 6.787807 5.335962 5.250163 6.47423 5.882104 5.965109 6.591687195 AFFX-BioB-M_at 6.163227 5.965427 4.665569 2.743531 6.097244

On 10/11/2006 2:48 PM, Seth Falcon wrote: > Hi, > > Reading help("Documentation"), I'm led to believe that a help call > like: > > ?myFun(x, sqrt(wt)) > > Will search for help on the appropriate method in the case that myFun > is generic. This isn't working for me. Here is an example using the > Biobase package: > > ## If Biobase is

I am very new to R. I was trying to load some publicly available Expression data in to R. I used the following commands mydata<-read.table("dataALLAMLtrain.txt", header=TRUE, sep ="\t",row.names=NULL) It reads data without any error Now if I use edit(mydata) It shows only 3916 entries, whereas the actual file contains 7129 entries) My data is something like Gene Description

Hi: I have clustered microarray gene expression data and trying to map between microarray probe, gene, pathway, gene ontology, and homology for a set of (affy) microarray probes. Is there any package in R which facilitates this? I am looking at bioconductor, but till now could not find a solution. A link to some worked example would be appreciated. Thanks and regards. John [[alternative HTML

Dear R-helpers & bioconductor Sorry for cross-posting, this concerns R-programming stuff applied on Bioconductor context. Also sorry for this long message, I try to be complete in my request. I am trying to write a subset method for a specific class (ExpressionSet from Bioconductor) allowing selection more flexible than "[" method . The schema I am thinking for is the following:

Hi, I have a list of IPI gene IDs. I want to find out whether there is a package which can map the gene ontology to these IPIs, and plot the pie chart to demonstrate the molecular function distributions. The input is like the following gene IPI IDs: IPI:IPI00008860.1|SWISS-PROT:Q9BXJ4-1|TREMBL:Q542Y2|ENSEMBL:ENSP00000231338;EN

Hi Is it possible to open a connection to an FTP site such that I can read the directory listing? Eg: URL <- url("ftp://ftp.ensembl.org", open="r") Error in url("ftp://ftp.ensembl.org", open="r") : unable to open connection URL <- url("ftp://ftp.ensembl.org") open(URL) Error in open.connection(URL) : unable to open connection I think

Hi, I can run the code some days ago . But cant run now.  Problem 1: Output is ok ensembl = useDataset("hsapiens_gene_ensembl",mart=ensembl) utr5 = getSequence(chromosome=3, start=185514033, end=185535839, type="entrezgene",seqType="5utr", mart=ensembl)  Output :                                                                                                5utr

Steffen, When the new biomaRt tries to load it errors out because I do not have RMySQL installed. There is not a Windows binary for RMySQL and it does contain C code that I do not know how to build. I do not use the MySQL option in biomaRt. Does RMySQL need to be a required dependency? Below is my screen output and sessionINfo. require(biomaRt) Loading required package: biomaRt Loading required