similar to: samr - extract genes from siggenes.table

Hi, here's my siggenes.table$genes.up snippet. Two class unpaired SAMR analysis. "Row" "Gene ID" "Gene Name" "Score(d)" "Numerator(r)" "Denominator(s+s0)" "Fold Change" "q-value(%)" "1" "25" "RPL15P22" "RPL15P22" "-1.44115338424578" "-18"

Dear All: Significance Analysis of Microarrays(SAM) As we know sam do multiple t.test as following ## Default S3 method: t.test(x, y = NULL, alternative = c("two.sided", "less", "greater"),mu = 0, paired = FALSE, var.equal = FALSE,conf.level = 0.95, ...) var.equal: a logical variable indicating whether to treat the two variances as being equal. If 'TRUE'

Dear R users, I'am using marray and Limma packages to analyze genepix output. 1) how can I filter bad spots from my data (data contains 3 types of bad spots). my experiment contains 12 samples and the bad spot are not associated to the same probes 2) how can I remove control probes from my data ? I'm sorry, i'm new with R and I can't find answer in packages doc. best regards,

Hello Everybody, I am trying to perform SAM with the samr package. I am using the following code: sink ("R005") library(siggenes) library(samr) library(nnet) A <- as.matrix(read.table("D:\samrgenes1000.txt")) B <- as.matrix(read.table("D:\genenames1000.txt")) y1 <- c(rep(1,20),rep(2,6)) #there are 20 chips of one kind and 6 of the other kind. datasam =

Hi, I have a problem on my mac when trying in R to produce png images. I am getting this warnings with the ArrayQualityMetrics package: > arrayQualityMetrics(rma_fatBody, outdir="normData", force =T) The report will be written into directory 'normData'. KernSmooth 2.23 loaded Copyright M. P. Wand 1997-2009 (loaded the KernSmooth namespace) libpng warning: Application built

Dear all, I am running R 2.4.1. > library(siggenes); > library(multtest); > cl<-rep(c(0,1),c(3,3)); > sub<-exprs(AffyExpData[,c(1:3,7:9)]); > gn<-geneNames(AffyRAwData); > sam.out<-sam(sub,cl,rand=123,gene.names=gn); We're doing 20 complete permutations > sam.out SAM Analysis for the Two-Class Unpaired Case Assuming Unequal Variances Delta p0

Not sure if this is an R or BioC question so will cross-post. I cannot extract R code from vignettes because system.file("doc", "my.bioc.package") returns "". I got this code directly off of http://www.bioconductor.org/docs/vignettes.html. A specific example using siggenes follows but I have tested this using multiple packages, including those for which I can

Hi, I am using SAM (from siggenes_1.2.11 package) method to select genes from a microarray data set. After installing the latest R2.4.0 on my computer, to my surprise the results are totally different from that calculated using R2.1.1. Even the example code doesn't work the same way under these two versions of R. Does anybody know what is going on? Thanks for any suggestions.

Dear RUsers, I am new to R. I am learning how to use R. I am a PC user and run R on windows. I would appreciate if some one could guide me on a few questions I have: 1) I have 4 cel files (2 replicates for NORM and Disease resp). I have been able to run siggenes on this dataset where I have 4 labels in the class file groupsnhi.cl op-> (0,0,1,1) and my data has been read into datrmanhi after

Hi,when I perform SAM on my array data(siggenes)I have some problems in retrieving the separate lists of up regulated and down regulated genes. When I write: fold<-function(x){ gruppi<-split(x,controllo) geni1<-abs(mean(gruppi[[2]])-mean(gruppi[[1]])) return(geni1) } fold<-esApply(expr.contr.tratt.4,1,fold)

Hi, I was working on a classification problem using the pamr package. I used the pamr.adaptthresh() function to find the optimal accuracy of the classifier. I must not be doing it right, since it doesn't return the threshold values for optimum classification. For example,if I run it on a dataset, I get the following result using pamr.adaptthresh(): predicted true (1)

Hi, Sorry about the earlier formatting errors... I was working on a classification problem using the pamr package. I used the pamr.adaptthresh() function to find the optimal accuracy of the classifier. I must not be doing it right, since it doesn't return the threshold values for optimum classification. For example,if I run it on a dataset, I get the following result using

Dear RUsers, I am new to R. I am learning how to use R. I am a PC user and run R on windows. I would appreciate if some one could guide me on a few questions I have: 1) I have 4 cel files (2 replicates for NORM and Disease resp). I have been able to run siggenes on this dataset where I have 4 labels in the class file groupsnhi.cl op-> (0,0,1,1) and my data has been read into datrmanhi after

Suppose I have two var x and y,now I want to fits a natural cubic spline in x to y,at the same time create new var containing the smoothed values of y. How can I get it?

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

Dear R-help fellows good afternoon. I am struggling in the attempt to impose some graphical conditions (changing point symbols, colors, etc) to biplot function (I am using it to visualize the results of princomp) but I can't apparently manage to change anything but the axis... and I have been browsing manuals and vignettes without finding any explicit suggestions on how to operate... Can

Hallo, I just installed all needed packages for my project on my PC. But I cannot load all at one time. I now want to load limma. How can I realize the following plan: I want to install for example limma inclusive all needed other sub packages (add-on). Can anyone tell me the corresponding command? Thanks, Corinna Here is the result of the command library(): Pakete in Library

Dear BioConductor and R fellow users I apologize in advance for double posting, but I am not sure which list would actually be best fit for this message. I am experiencing a weird error with my R installation on Ubuntu 10.04.4 (LTS) 64bit: When I run R on the terminal everything goes smoothly: $R R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation

Hello list! I have a proble trying to perform a SAM analysis using the function samr from the samr package. I have put the option *center.arrays=TRUE *in order to scale all the experiments to median=0. I would like to retrieved the scaled data but it seems that samr does not return it...Does anyone have any idea on this? Thanks a lot!!! E. [[alternative HTML version deleted]]

Hi R-help, I am trying to use 'samr' for 10 pre and post paired samples to test whether post is different from pre (i.e., the location shift for the delta of (post-pre)). However, I got an error message saying > samr.obj<-samr(d, resp.type="Two class paired", nperms=100, random.seed=100) perm= 1 Error in !logged2 : invalid argument type Does anyone know what this