[R] illumina array - remove probes in discontinued genes

Hi I'm working with the illumina HT-12v4 array and after using LIMMA some of the most significant probes that come up are in genes that have been discontinued on NCBI (ie LOC732450). Is there a way to filter my data of these probes prior to statisitcal analysis? Thanks, .kripa

[R] illumina probeID to entrezID

Hi, how do you convert the illumina ProbeID to entrezID? I've used LIMMA to identify top genes but i'm not sure how to apply illuminaHumanv4.db to get the EntrezID. Thanks in advance, .kripa [[alternative HTML version deleted]] _

[R] GLMNET warning msg

Hi I'm getting the following warning msg after ?cv.glmnet and I'm wondering what it means... dim(x) 10 12000; dim(y) 10; #two groups case=1 and control=0 cv.glmnet(x, y) Warning message: Option grouped=FALSE enforced in cv.glmnet, since < 3 observations per fold Thanks, .kripa

[R] PAMR - adding additional labels

Hello, I recently started using PAM to analyze my microarrays (http://cran.r-project.org/web/packages/pamr/pamr.pdf). Firstly does anyone know of a specific forum or help section for this package? Currently with ?pamr.from.excel you can only have one line of labels, ie control vs treatmen

Re: [R] pulling out pairs from data frame

0 M > A.K. > > > From: Kripa R > To: arun > Cc: R help > Sent: Tuesday, August 13, 2013 10:56 AM > Subject: RE: [R] pulling out pairs from data frame > > > > > Oops! Ok So I have this file: > > SampleName Ind

Re: [R] pulling out pairs from data frame

%in% reps) > ádat1 > # áSameName Individual Age Gender Newcol > #1 á á á á1 á á á á á4 á80 á á áM á á á1 > #2 á á á á2 á á á á 15 á56 á á áF á á á1 > #3 á á á á3 á á á á á1 á75 á á áF á á á0 > #4 á á á á4 á á á á 15 á56 á á áF á á á1 > #5 á á á á5 á á á á á2 á58 á á áF á á á0 &g

[R] pulling out pairs from data frame

Hello everyone, I'm having trouble pulling out paired samples from a data set... I have the following: reps<-c(4,15) #the variable reps is a list of all paired samples >data SameName Individual Age Gender 1 4 80 M

[R] order function

Hi I'm having a simple issue with the order function. When I use the following code my data is not ordered correctly > output[order(output[,3]),] Namebeta pval 881 9.09303277751237 0.000100253350483199 74026.40553461638365 0.00010228641631914 4879

[R] order function

Hi I'm having a simple issue with the order function. When I use the following code my data is not ordered correctly > output[order(output[,3]),] Namebeta pval 881 9.09303277751237 0.000100253350483199 74026.40553461638365 0.00010228641631914 4879 -