SEQanswers

Go Back   SEQanswers > Introductions



Similar Threads
Thread Thread Starter Forum Replies Last Post
SPP -- ChIP-seq program installation?? ewilbanks Bioinformatics 10 06-19-2013 06:41 AM
Staden package question: is srf2fastq returning an error, or is this normal? mhayes Bioinformatics 4 12-12-2011 10:54 AM
Problems with SPP read.bams.tags EdinG Bioinformatics 0 10-12-2011 02:27 AM
how to install spp chip-seq package in linux? RVJ12 Bioinformatics 0 02-21-2011 11:01 AM
Install spp (chip-seq processing pipeline) bbl Bioinformatics 0 10-27-2010 07:02 AM

Reply
 
Thread Tools
Old 11-15-2009, 03:17 AM   #1
abmmki
Junior Member
 
Location: Africa

Join Date: Nov 2009
Posts: 5
Default Error -- SPP package, MSER -- could you please help

Hi,

I have some ChIp-seq data (Illumina v1.4; and ELAND extended mapped) and trying to analyze them using spp R package (http://www.nature.com/nbt/journal/v2.../nbt.1508.html) and package (http://compbio.med.harvard.edu/Supplements/ChIP-seq/). I am following the tutorial and it was going fine until MSER step which I really need. I got the following error:

mser <- get.mser(RLTA_chip.data, RLTA_input.data, step.size=1e5, test.agreement=0.99, n.chains=10, cluster=NULL, fdr=0.05, method=tag.wtd, whs=detection.window.halfsize)

excluding systematic background anomalies ... done
calculating statistical thresholds
FDR 0.05 threshold= Inf
chained subsampling using fraction 0.9510645 .

Error in ecdf(-mpd$re) : 'x' must have 1 or more non-missing values

In addition: Warning messages:
1: In min(npld$y[npld$fdr <= fdr]) :
no non-missing arguments to min; returning Inf
2: In min(npld$y[npld$fdr <= fdr]) :
no non-missing arguments to min; returning Inf
3: In min(npld$y[npld$fdr <= fdr]) :
no non-missing arguments to min; returning Inf
4: In is.na(mpd$re) :
is.na() applied to non-(list or vector) of type 'NULL'
5: In is.na(mpd$oe) :
is.na() applied to non-(list or vector) of type 'NULL'




> sessionInfo()
R version 2.9.0 (2009-04-17)
x86_64-unknown-linux-gnu

locale:
LC_CTYPE=en_US.UTF-8;LC_
NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=C;LC_MESSAGES=en_US.UTF-8;LC_PAPER=en_US.UTF-8;LC_NAME=C;LC_ADDRESS=C;LC_TELEPHONE=C;LC_MEASUREMENT=en_US.UTF-8;LC_IDENTIFICATION=C

attached base packages:
[1] stats graphics grDevices utils datasets methods base

other attached packages:
[1] nws_1.7.0.0 snow_0.3-3 spp_1.0 caTools_1.9 bitops_1.0-4.1

loaded via a namespace (and not attached):
[1] tools_2.9.0

----------------

Could you please help to understand where is the error........................and how can I solve this. (I have a data set where control has much much more tags then in sample IP................14 million versus 2 million.....................so, I want to use the mser)

Also in get.mser function: >> method=tag.wtd ............................I wanted to change it like method=tag.lwcc........................like following:



mser <- get.mser(RLTA_chip.data, RLTA_input.data, step.size=1e5, test.agreement=0.99, n.chains=10, cluster=NULL, fdr=0.05, method=tag.lwcc, whs=detection.window.halfsize, skip.control.normalization = F)

I got following error:


chained subsample step 0 :
finding background exclusion regions ... done
determining peaks on provided 1 control datasets:
using reversed signal for FDR calculations
bg.weight= 2.936861 processing chr1 in 7 steps [Error in lwcc(x, y, s, e, return.peaks = return.peaks, bg.x = bg.x, bg.y = bg.y, :
unused argument(s) (skip.control.normalization = TRUE)


I did not find any documentation about skip.control.normalization....................................tried to write tha "FALSE"................but it then gave me two un used arguments (TRUE and FALSE)...........................so, how to solve this problem ??


Thank you in advance for your generous help.

may be you can ans in my mail directly:

Khademul
abmmki@gmail.com
abmmki is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 02:39 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO