similar to: Overlay Gene Expression on SNP (copy number) data

Dear All, Thank you kindly for such detailed replies. I was looking to overlay data using algorithms so that i am able to tell which genes are differentially expressed due to changes in copy number. I did a pubmed search and found only 7 literature pieces all of which use in-house algorithms. I am yet to explore Gviz since it wouldn't work on R 2.14, would try it after upgrading to R 2.15.

Dear All, I have an exon array and did not find any differential gene expression between two samples. I was looking to perform correlation analysis on the same. Can anyone recommend any package that would do this for an affy exon array? Will SAM analysis give me correlated genes? Thanks and regards, Ekta The information contained in this electronic message and in any attachments to this message

Dear All, I am using the LIMMA package to create 2 contrasts for my data and then calculating the vennCounts of the decideTests from the contrast.fit to be able to create venn Diagrams. The code works fine but the summary(results) shows zeros for all i.e. no gene were up regulated or downregulated. This is not true for my data since toptable output shows Log fold change greater than > 2. I am

Dear Bioc and R List Users, I am having trouble analysing illumine data generated from BeadScan. I have .idat files and JPEG images. I realise that i need bead-level summary data to be able to begin quality control followed by normalization. Is there a way i can read .idat files for expression analysis or do i need to go back to BeadScan and generate .txt files/tiff files ? Appreciate any help

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

Hi, I was using Gviz package to create a boxplot. I understand that Gviz uses "panel.bwplot" to create the boxplot. Is there any way that I can remove the dashed line surrounding each pair of boxplots? Here is some sample code: ############# library(Gviz) thisdata <- matrix(sample(1:100,60),nrow=10,ncol=6) positions <- sample(1:100,6) limit1 <- min(positions)-1 limit2

Hi, I've noticed an issue with S4 methods and namespaces which only arises in particular, difficult to reproduce configurations. One example is the ggbio package in Bioconductor, which currently emits these warnings when its namespace is loaded: ---------------------- library(ggbio) Loading required package: ggplot2 Attaching package: ‘ggbio’ The following object(s) are masked from

Hi I want to simplify my problem into a prototype, how to write a function to get all the object in your current environment, e.g. .GlobalEnv, and print their mode? For example, if I have object a,b,c... in my environment, a=1;b='test';c=matrix(0,3,3)..... I want to write a function myfun(), when I run myfun(), this retrieve all the objects in my environment automatically and print

Hi dear all, In my laptop(ubuntu 11.10 64bit), I maintained a released R (2.14) and a developmental R, I can specify qplot(..., alpha = ) in R 2.14 , but when I try to use transparency in developmental R, I got a warning message and the plot is clearly not I want. minimal example: > qplot(data = mtcars, x = mpg, y = cyl, alpha = cyl) Warning message: In grid.Call.graphics(L_points, x$x, x$y,

Hi dear all, In my laptop(ubuntu 11.10 64bit), I maintained a released R (2.14) and a developmental R, I can specify qplot(..., alpha = ) in R 2.14 , but when I try to use transparency in developmental R, I got a warning message and the plot is clearly not I want. minimal example: > qplot(data = mtcars, x = mpg, y = cyl, alpha = cyl) Warning message: In grid.Call.graphics(L_points, x$x, x$y,

Hi dear all I am currently using Ubuntu 9.10 (karmic) and R 2.10, I can find a R release that support "karmic" in cran under ubuntu directory, I plan to update my laptop to Ubuntu 10.04 after April 30, since I have to use some new features or new version in it, e.g. QT 4.6, I don't know if there is going to be a R release that will support this newly released Ubuntu? Thanks in

Hi dear all, I am wondering if there is a function existing in R that did the quadratic bezier curve interpolation? I hope to generate a bezier curve based on three sets of points: two end of the line and a control point. Thanks in advance. Tengfei -- Tengfei Yin MCDB PhD student 1620 Howe Hall, 2274, Iowa State University Ames, IA,50011-2274 [[alternative HTML version deleted]]

I am really new with R Graphical user interfacefunctions. I am developing a software package to calculate pKa (biochemistry) but I want to make it look aesthetically pleasing and make it user friendly. I have heard that R has some GUI (Graphical user interface) and you can do some really cool stuff out there. What are the limitations and what are some resources for help. I have found a couple of

Hi, list- ?? I've got a large list of multi-column data and I'd like to filter it according to a threshold, such as, "show me the values that are above 0.4 for each line." For instance, if the line of data were: > line v1 v2 v3 v4 v5 v6 v7 v8 v9 1 -0.32 0.66 -0.35 -0.82 0.38 0.66 -0.02 -0.11 -0.64 I can do this > line[,line >= 0.4]

Hello, I have a data.frame: name    col1    col2    col3    col4 AA    23    54    0.999    0.78 BB    123    5    1    0.99 AA    203    98    0.79    0.99 I want to get mean value data.frame in terms of name: name    col1    col2    col3    col4 AA    113.0000  76.0000   0.8945   0.8850 BB    123.00   5.00   1.00   0.99 I tried to use by function: >aa<-by(test[,2:5], feature, mean)

I recently started using ubuntu 9.10 and I am using gedit editor and R plugin for writing R code. To install any package I need to do: $ install.packages() //window pop-up for mirror selection //then another window pop up for package selection. After this as long as I am not exiting, the function of the newly installed packages are available. After I exit (i use to put 'no' in 'save

Hi dear all I am currently using Ubuntu 9.10 (karmic) and R 2.10, I can find a R release that support "karmic" in cran under ubuntu directory, I plan to update my laptop to Ubuntu 10.04 after April 30, since I have to use some new features or new version in it, e.g. QT 4.6, I don't know if there is going to be a R release that will support this newly released Ubuntu? Thanks Best

Hello, I would like to plot specific SNPs with their exact locations on a chromosome. Based on my genotyping results I would like to separate these SNPs in three different categories: 1, 2 and 3 and use different colours to represent these categories. The script below generates the sample data. I can plot these with the image function using the following: val <- 1:3 samp <- sample(val,

Dear All, The following code extracts from NCBI very nice output for ONE allele of a SNP (often the allele with the second largest frequency - usually termed the minor allele). It gives an average minor allele frequency from all NCBI sources (which is what I want, except I'd like the addition of data for all the other alleles of one SNP) plus a table of minor allele frequencies from each

Hi, Does anybody know how to obtain the imputed SNP genotype probabilities from the snpStats package? I am interested in using an imputation method implemented in R to be further used in a simulation study context. I have found the snpStats package that seems to contain suitable functions to do so. As far as I could find out from the package vignette examples and its help, it gives the