[BioC] SnapCGH problem

Daniel Brewer daniel.brewer at icr.ac.uk
Thu Feb 8 16:36:02 CET 2007


I having trouble setting up Agilent 4100a array data to be used with
snapCGH.  This is what I do:

#Input file
targets <- readTargets()
stuff <- read.maimages(targets$FileName,source="agilent")

#Setup layout
stuff$genes$Block <- 1
names(stuff$genes)[2] <- "Column"
stuff$printer <- getLayout(stuff$genes)

#setup clone information
locData <- read.table(file="4100aLocationShort.csv", header=TRUE, sep="\t")
names(locData) <- c("ProbeName","Chromosome","Start","End")
stuff$genes <- merge(stuff$genes,locData, all.x=TRUE)
stuff$genes$Position <- stuff$genes$Start
stuff$genes$Chr <- stuff$genes$Chromosome

#Setup design (from snapCGH guide)
stuff$design <- c(-1,-1)

#Normalisation within
normed <- normalizeWithinArrays(stuff, method="loess")

#Get rid of "genes" with no location information
normed$genes <- normed$genes[!(is.na(normed$genes$Chr)) &
normed$genes$Chr != "",]

#ProcessCGH
normed2 <- processCGH(normed, method.of.averaging = mean,ID="ProbeName")

The result is:
> Averaging duplicated clones
> Processing chromosome  NA 
..
> Processing chromosome  NA 
> Warning messages:
> 1: > not meaningful for factors in: Ops.factor(MA$genes$Chr, maxChromThreshold) 
> 2: < not meaningful for factors in: Ops.factor(MA$genes$Chr, minChromThreshold) 
> 3: <= not meaningful for factors in: Ops.factor(uniq.chrom, maxChrom)

I was wondering whether you could help me solve this.  I can think of
two reasons why this might be occurring,
1) I have not deleted the data corresponding to the genes I have
removed.  Not sure how to do that.
2) My genes table is not set up correctly.  Here is a sample of the table:
    ProbeName Row Column ProbeUID ControlType
873   1000182  57     33     5912           0
874     10003   1     54       27           0
                                                                   GeneName
873                    actin related protein 2/3 complex, subunit 3 (21 kD)
874 Human (clone pAT 464) potential lymphokine/cytokine mRNA, complete cds.
    SystematicName Description Block Chromosome     Start       End
Position
873       BG483858      pSPORT     1         12 109335427 109350878
109335427
874         M25315       pBlue     1         17  31439709  31441605
31439709
    Chr
873  12
874  17


I appreciate your help

Daniel

-- 
**************************************************************

Daniel Brewer, Ph.D.

Institute of Cancer Research

Email: daniel.brewer at icr.ac.uk

**************************************************************

The Institute of Cancer Research: Royal Cancer Hospital, a charitable Company Limited by Guarantee, Registered in England under Company No. 534147 with its Registered Office at 123 Old Brompton Road, London SW7 3RP.

This e-mail message is confidential and for use by the addre...{{dropped}}



More information about the Bioconductor mailing list