similar to: [BioC] How to do RMA without summary to probeset level?

Hi, I am a Ph.D. student from Québec, Canada. I’m a beginner with R and Bioconductor. Until now the only experience I have is in analyzing microarray data using affy and limma packages. Now I am trying to analyze Rat Gene 10 st arrays and I would like to run RMA analysis and Smyth moderated t test on those arrays. Since no cdf official package is available for those arrays, after reading many

Dear Help, After loading the pd.Citrus library and checking the DataFrame, I ran > the R code for: > > 1) 'oligo' > > > > {> library(pd.citrus) > Loading required package: RSQLite > Loading required package: DBI > > data(pmSequence) > > > show(pmSequence) > DataFrame with 341730 rows and 2 columns > fid sequence > <integer>

Hi all! While I am trying to read .cel files with oligo package: afbatch=read.celfiles(list.celfiles()) I get an error: Package 'pd.mirna.1.0.2xgain' was not found in the BioConductor repository How can I overcome this? Thank you in advance

Dear all, I've noticed by trying to download gz files from here : https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM907811 At the bottom one can download GSM907811.CEL.gz . If I download this manually and try oligo::read.celfiles("GSM907811.CEL.gz") everything works fine. (oligo is a bioConductor package) However, if I download using download.file("

On Wed, 21 Oct 2020 16:15:22 -0500 Ana Marija <sokovic.anamarija at gmail.com> wrote: > Hello, > > I have a data frame with one column: > > > remove > > V1 > > 1 ABAFT_g_4RWG569_BI_SNP_A10_35096 > 2 ABAFT_g_4RWG569_BI_SNP_B12_35130 > 3 ABAFT_g_4RWG569_BI_SNP_E09_35088 > 4 ABAFT_g_4RWG569_BI_SNP_E12_35136 > 5

Hello, To remove the file extension it's much easier to use base R filename <- tools::file_path_sans_ext(basename(celFiles)) Hope this helps, Rui Barradas ?s 22:41 de 21/10/20, Rui Barradas escreveu: > Hello, > > This is probably because basename keeps the file extension, try instead > > > filename <- sub("(^[^\\.]*)\\..+$", "\\1",

Hello, This is probably because basename keeps the file extension, try instead filename <- sub("(^[^\\.]*)\\..+$", "\\1", basename(celFiles)) celFiles[!filename %in% as.character(remove$V1)] Hope this helps, Rui Barradas ?s 22:15 de 21/10/20, Ana Marija escreveu: > Hello, > > I have a data frame with one column: > >> remove > >

Hello, I have a data frame with one column: > remove V1 1 ABAFT_g_4RWG569_BI_SNP_A10_35096 2 ABAFT_g_4RWG569_BI_SNP_B12_35130 3 ABAFT_g_4RWG569_BI_SNP_E09_35088 4 ABAFT_g_4RWG569_BI_SNP_E12_35136 5 ABAFT_g_4RWG569_BI_SNP_F11_35122 6 ABAFT_g_4RWG569_BI_SNP_F12_35138 7 ABAFT_g_4RWG569_BI_SNP_G07_35060 8 ABAFT_g_4RWG569_BI_SNP_G12_35140 I want to remove these 8

Hello: I am trying to work with a couple of microarray data sets, using platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 1 minor 8.1 year 2003 month 11 day 21 language R In the shortcut for invoking R I have set

Hi, I am relatively new to R and Bioconductor and am trying to filter the topTable that I generated of differentially expressed genes from my normlized eset file comprised of ~ 40 HG-133A Affy microarrays . I would like to see if particular probesets are represented in this list. Alternatively I would like to generate a topTable of differentially expressed genes using only specified probesets

Hi Guys, I am having a real hard time trying to figure out for microarry. Here is my code One-Sample t-Test dim(data.sub) [1] 10000 140 ##there are 10000 probesets and 140 columns hist(data.sub) ## Histogram. Identify if the probesets are normal distributed q<-rnorm(10000) ##generate 10000 random, normal distributed values qqplot(data.sub,q)) ##Show the plot of the probeset

On 05/02/2018 03:21 PM, Joris Meys wrote: > Dear all, > > I've noticed by trying to download gz files from here : > https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSM907811 > > At the bottom one can download GSM907811.CEL.gz . If I download this > manually and try > > oligo::read.celfiles("GSM907811.CEL.gz") > > everything works fine. (oligo

Use mode="wb" when you download the file. See https://github.com/HenrikBengtsson/Wishlist-for-R/issues/30. R core, and others, is there a good argument for why we are not making this the default download mode? It seems like a such a simple fix to such a common "mistake". Henrik On Thu, May 3, 2018, 00:44 Joris Meys <jorismeys at gmail.com> wrote: > Dear all, >

Using the correct mode absolutely solves it. Apologies for not trying the obvious. Cheers Joris On Thu, May 3, 2018 at 2:10 PM, Martin Morgan <martin.morgan at roswellpark.org > wrote: > > > On 05/02/2018 03:21 PM, Joris Meys wrote: > >> Dear all, >> >> I've noticed by trying to download gz files from here : >>

Dear all, I've been diving a bit deeper into this per request of Tomas Kalibra, and found the following : - the lock on the file is only after trying to read it using oligo, so that's not a R problem in itself. The problem is independent of extrenal packages. - using Windows' fc utility and cygwin's cmp utility I found out that every so often the download.file() function inserts

Hi, I have a question concerning the analysis of some affymetrix chips. I downloaded some of the data from GEO GSE11324 (see below). In doing so I'm stuck after I identified the probesets with significant changes. I have problems in assigning probeset specific gene names as well as getting the genomic coordinates. Furthermore I have no clue how to deal with the fact, that most genes have

Hi everyone, Firstly, I'm sorry for the cross-post. I re-read the posting guide and it appears my question is more related to r-devel.... (i think :-) ) I downloaded the source code available at: http://cran.fhcrc.org/src/base/R-2/R-2.4.0.tar.gz to a linux machine (Linux 2.6.11-1.1369_FC4smp #1 SMP). I successfully configured and compiled it, which means I'm able to run it from the

Hello, I want to use expresso for preprocessing the hgu133a-spikein data from affycompII. But there is an error: > library(affy) > path <- "z:/Microarray/hgu133a-spikein/rawdata" > celFile <- list.celfiles(path=path,full.names=TRUE); > affyBatch <- ReadAffy(filenames=celFile[1:6]); > eset1 <-

Hi everyone, I downloaded the source code available at: http://cran.fhcrc.org/src/base/R-2/R-2.4.0.tar.gz to a linux machine (Linux 2.6.11-1.1369_FC4smp #1 SMP). I successfully configured and compiled it, which means I'm able to run it from the R-2.4.0/bin directory. I want to do a system wide installation (via 'sudo make install') and help2man fails, ie, when I execute the

dear members, i would like to write a variable in a plot title (main="") but i don't know the right syntax:(...i tried a lot of different ways without success. here my example: y=30 z=33 for (i in 10:length(tissue)) { png(filename = tissues[i], width = 1024, height = 768, pointsize = 12, bg = "white") gene.graph("ENSG00000115252", rma.affy, gps=list(1:3,