Displaying 1 result from an estimated 1 matches for "nisland".
Did you mean:
island
2007 Dec 19
1
FW: cgh package
....table() for the islands
location, and the best score, ideally, the ones corresponding to the
island(s) plotted by the sw.plot().
Here is the exemple code of what I used:
x <- sw.threshold(intens[,8], function(x) median(x) + .2 * mad(x), -1)
r <- sw.rob(intens[,8])
sw <- sw(x,max.nIslands = NULL, trace = F)
sw.perm <- sw.perm.test(x, max.nIslands = NULL, nIter= 1e4)
sw.plot(intens[,8], location = coords[,2],threshold.func = function(x)
median(x) + .2 * mad(x), sign = -1, expected = NULL, legend =
TRUE,rob=r,main=paste("chrom 9 loss p=",sw.perm[1]))
sessionInfo()
R...