[BioC] Extract M values from probes set to cluster

Marcelo Laia marcelolaia at gmail.com
Tue Feb 13 21:53:06 CET 2007


Hi Jan,

I try your suggestion and some errors appear.

Could you (or someone) help me again?

Thank you very much

> fit.2 <- eBayes(fit)
> selected  <- (fit.2$p.value) <0.05
>
> selected[is.na(selected)]<-FALSE
>
> MA.selected <- MA.2$M[selected, ]
>
> M.uw<-unwrapdups(MA.2$M,ndups=2,spacing=32)
>
> M.ag<-apply(M.uw,1,function(x) (x[c(T,F)]+x[c(F,T)])/2)
>
> MA.selected<-M.ag[selected,]
Error: (subscript) logical subscript too long
> length(selected)
[1] 3456
> length(M.uw)
[1] 20736
> length(M.ag)
[1] 10368
> dim(MA.selected)
[1] 2632    3
> dim(M.uw)
[1] 3456    6
> dim(M.ag)
[1]    3 3456
>

Marcelo

2007/2/13, J.Oosting at lumc.nl <J.Oosting at lumc.nl>:
> First of all, you will have to get rid of the missing values in your
> selection vector
>
>  selected[is.na(selected)]<-FALSE
>
> Secondly, I think that you have duplicate spots on your array, with a
> total of 6912 spots. It depends on the layout on the chip how you can
> get to the M-values of both spots using the selected vector.
> unwrapdups() can help you here
> Something like this (untested)
> M.uw<-unwrapdups(MA.2$M,ndups=2,spacing=1) # ndups and spacing same as
> in lmFit() and duplicateCorrelation()
> M.ag<-apply(M.uw,1,function(x) (x[c(T,F)]+x[c(F,T)])/2)  # doesn't
> handle NA
> MA.selected<-M.ag[selected,]
>
> Jan
>
>
>
> -----Original Message-----
> From: bioconductor-bounces at stat.math.ethz.ch
> [mailto:bioconductor-bounces at stat.math.ethz.ch] On Behalf Of Marcelo
> Laia
> Sent: dinsdag 13 februari 2007 12:51
> To: Bioconductor
> Subject: Re: [BioC] Extract M values from probes set to cluster
>
> Dear Sean
>
> I thought that I did a mistake in selected. Then, I try... But no
> success.
>
> Here are the outputs:
>
> > selected  <- p.adjust(fit.2$p.value) <0.05
> > summary(selected)
> Mode   FALSE    TRUE    NA's
> logical    2565     130     761
> > length(selected)
> [1] 3456
> > MA.selected <- MA.2$M[selected, ]
> > dim(MA.selected)
> [1] 1782    3
>
> I thought that I did a mistake here:
>
> > selected  <- p.adjust(fit.2$p.value) <0.05
>                       ^^^^^^^^^^
>
> > selected  <- (fit.2$p.value) <0.05
> > length(selected)
> [1] 3456
> > MA.selected <- MA.2$M[selected, ]
> > dim(MA.selected)
> [1] 4154    3
> > fit.2 <- eBayes(fit)
>
> Thank you very much
>
> Laia
>


-- 
Marcelo Luiz de Laia
Ph.D Candidate
São Paulo State University (http://www.unesp.br/eng/)
School of Agricultural and Veterinary Sciences
Department of  Technology
Via de Acesso Prof. Paulo Donato Castellane s/n
14884-900   Jaboticabal - SP - Brazil
Phone: +55-016-3209-2675
Cell: +55-016-97098526



More information about the Bioconductor mailing list