[aroma.affymetrix] CbsModel or writeRegions error

classic Classic list List threaded Threaded
6 messages Options
Reply | Threaded
Open this post in threaded view
|

[aroma.affymetrix] CbsModel or writeRegions error

Emanuel Gonçalves
Hi all,

I'm processing a large set of SNP6 samples (1020) and performing CBS segmentations afterwards. It's a very lengthy process and the annoying part is that it systematically crashes on the last sample, see code and output below:

# Run CBS paired segmentation with pooled normal samples
cbs
<- CbsModel(cesN, cesN1.ref)
fit
(cbs, verbose=verbose)
print(cbs)

# Export CBS regions
pathname
<- writeRegions(cbs, verbose=verbose)


20170831 06:57:51| Pathname: cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YKG-1,chr25,6c653767ed90d29ad5254a760d8f3952.xdr
20170831 06:57:51| Already done. Skipping.
20170831 06:57:51|Array #1018 ('YKG-1') of 1019 on chromosome 25...done
20170831 07:01:07|Genomic-signal tags:
20170831 07:01:07|Reference tags: 6c653767ed90d29ad5254a760d8f3952
20170831 07:01:08|Array #1019 ('YMB-1-E') of 1019 on chromosome 1...
20170831 07:01:08| Pathname: cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YMB-1-E,chr01,6c653767ed90d29ad5254a760d8f3952.xdr
20170831 07:01:08| Already done. Skipping.
Error: empty (zero-byte) input file
Execution halted

At first I thought it was a problem of the sample, but that is not the case as I could process it separately. In fact, at this stage all the samples seem to be processed. Also strange is that the script works and exports the regions for subsets of the data-set.

Any suggestion in how to export these would be much appreciated?

Thank you,

P.S: I'm using the multithreaded option of aroma

# - Setup configurations
# options(mc.cores = 10)
print(future::availableCores())

# Multithreaded
future
::plan('multiprocess')

# Increase ram
setOption
(aromaSettings, 'memory/ram', 600.0)

# Logs
log
<- verbose <- Arguments$getVerbose(-4, timestamp=TRUE)

# Reduce decimal places to minimize space
options
(digits=4)


--
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example.
 
 
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to [hidden email]
To unsubscribe and other options, go to http://www.aroma-project.org/forum/

---
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [hidden email].
For more options, visit https://groups.google.com/d/optout.
Reply | Threaded
Open this post in threaded view
|

Re: [aroma.affymetrix] CbsModel or writeRegions error

Henrik Bengtsson-2
Hi,

does the error occur when you run:

  fit(cbs, verbose=verbose)

or did that complete successfully and you get the error while running:

  pathname <- writeRegions(cbs, verbose=verbose)

If you run interactively, what does traceback() output if called
immediately after the error occurs?

Also, see if

fit(cbs, arrays = 1019:1020, chromosomes = c(1, 25), verbose = verbose)

or alternatively,

  pathname <- writeRegions(cbs, arrays = 1019:1020, chromosomes = c(1,
25), verbose=verbose)

gives the same error.  Then you can get to traceback() a bit sooner.

/Henrik

PS. I assume that you already know that rerunning the commands will
not redo the actual analysis; already processed samples will be
skipped - though with 1020 samples it might still take some time.


On Thu, Aug 31, 2017 at 1:26 AM, Emanuel Gonçalves
<[hidden email]> wrote:

> Hi all,
>
> I'm processing a large set of SNP6 samples (1020) and performing CBS
> segmentations afterwards. It's a very lengthy process and the annoying part
> is that it systematically crashes on the last sample, see code and output
> below:
>
> # Run CBS paired segmentation with pooled normal samples
> cbs <- CbsModel(cesN, cesN1.ref)
> fit(cbs, verbose=verbose)
> print(cbs)
>
> # Export CBS regions
> pathname <- writeRegions(cbs, verbose=verbose)
>
>
>> 20170831 06:57:51| Pathname:
>> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YKG-1,chr25,6c653767ed90d29ad5254a760d8f3952.xdr
>> 20170831 06:57:51| Already done. Skipping.
>> 20170831 06:57:51|Array #1018 ('YKG-1') of 1019 on chromosome 25...done
>> 20170831 07:01:07|Genomic-signal tags:
>> 20170831 07:01:07|Reference tags: 6c653767ed90d29ad5254a760d8f3952
>> 20170831 07:01:08|Array #1019 ('YMB-1-E') of 1019 on chromosome 1...
>> 20170831 07:01:08| Pathname:
>> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YMB-1-E,chr01,6c653767ed90d29ad5254a760d8f3952.xdr
>> 20170831 07:01:08| Already done. Skipping.
>> Error: empty (zero-byte) input file
>> Execution halted
>
>
> At first I thought it was a problem of the sample, but that is not the case
> as I could process it separately. In fact, at this stage all the samples
> seem to be processed. Also strange is that the script works and exports the
> regions for subsets of the data-set.
>
> Any suggestion in how to export these would be much appreciated?
>
> Thank you,
>
> P.S: I'm using the multithreaded option of aroma
>
> # - Setup configurations
> # options(mc.cores = 10)
> print(future::availableCores())
>
> # Multithreaded
> future::plan('multiprocess')
>
> # Increase ram
> setOption(aromaSettings, 'memory/ram', 600.0)
>
> # Logs
> log <- verbose <- Arguments$getVerbose(-4, timestamp=TRUE)
>
> # Reduce decimal places to minimize space
> options(digits=4)
>
>
> --
> --
> When reporting problems on aroma.affymetrix, make sure 1) to run the latest
> version of the package, 2) to report the output of sessionInfo() and
> traceback(), and 3) to post a complete code example.
>
>
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group with website http://www.aroma-project.org/.
> To post to this group, send email to [hidden email]
> To unsubscribe and other options, go to http://www.aroma-project.org/forum/
>
> ---
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to [hidden email].
> For more options, visit https://groups.google.com/d/optout.

--
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example.


You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to [hidden email]
To unsubscribe and other options, go to http://www.aroma-project.org/forum/

---
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [hidden email].
For more options, visit https://groups.google.com/d/optout.
Reply | Threaded
Open this post in threaded view
|

Re: [aroma.affymetrix] CbsModel or writeRegions error

Emanuel Gonçalves
Hi Henrik,

does the error occur when you run:

  fit(cbs, verbose=verbose)

Yes, because I never got the "print(cbs)" output
 
or did that complete successfully and you get the error while running:

  pathname <- writeRegions(cbs, verbose=verbose)

If you run interactively, what does traceback() output if called
immediately after the error occurs?

Unfortunately, I can't run this interactively because even though the samples are already preprocessed it takes almost a day to load everything.

Also, see if

fit(cbs, arrays = 1019:1020, chromosomes = c(1, 25), verbose = verbose)

or alternatively,

  pathname <- writeRegions(cbs, arrays = 1019:1020, chromosomes = c(1,
25), verbose=verbose)

gives the same error.  

I changed the code accordingly:

fit(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)

It ran without throwing any error (I tried before the same with a subset of ~100 samples and it also worked). 

It seems to be a specific problem when all the 1020 samples are ran together. I've placed a try/catch with the traceback().

# Run CBS paired segmentation with pooled normal samples
cbs
<- CbsModel(cesN, cesN1.ref)
out <- tryCatch(
 
{
    fit
(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)
 
},
  error
=function() {
    message
('Try/catch: error')
    traceback
()
 
},
  warning
=function() {
    message
('Try/catch: warning')
    traceback
()
 
},
 
finally={
    message
('Try/catch: finally')
    traceback
()
 
}
)
print(cbs)

I started the script with all the samples (without "arrays=1019:1020") as before.

What surprised me is that the script continued to the writeRegions call and it started to export all the samples.

pathname <- writeRegions(cbs, verbose=verbose)

> # Export CBS regions
> pathname <- writeRegions(cbs, verbose=verbose)
20170831 21:53:10|Array #1 ('201T') of 1019...
20170831 21:53:10| Extracting regions from all fits...
20170831 21:53:10|  Obtaining CN model fits (or fit if missing)...
20170831 21:57:15|  Obtaining CN model fits (or fit if missing)...done
20170831 21:57:15|  Extracting regions for chromosome #1...
20170831 21:57:15|  Extracting regions for chromosome #1...done
20170831 21:57:15|  Extracting regions for chromosome #2...
20170831 21:57:15|  Extracting regions for chromosome #2...done
20170831 21:57:15|  Extracting regions for chromosome #3...
20170831 21:57:15|  Extracting regions for chromosome #3...done

Does it mean it will export all the samples processed so far? Can I export everything without re-calling the fit function from CbsModel?

Thanks a lot, 
 
Then you can get to traceback() a bit sooner.

/Henrik

PS. I assume that you already know that rerunning the commands will
not redo the actual analysis; already processed samples will be
skipped - though with 1020 samples it might still take some time.


On Thu, Aug 31, 2017 at 1:26 AM, Emanuel Gonçalves
<<a href="javascript:" target="_blank" gdf-obfuscated-mailto="ve6ut0lsBwAJ" rel="nofollow" onmousedown="this.href=&#39;javascript:&#39;;return true;" onclick="this.href=&#39;javascript:&#39;;return true;">emanuelv...@...> wrote:

> Hi all,
>
> I'm processing a large set of SNP6 samples (1020) and performing CBS
> segmentations afterwards. It's a very lengthy process and the annoying part
> is that it systematically crashes on the last sample, see code and output
> below:
>
> # Run CBS paired segmentation with pooled normal samples
> cbs <- CbsModel(cesN, cesN1.ref)
> fit(cbs, verbose=verbose)
> print(cbs)
>
> # Export CBS regions
> pathname <- writeRegions(cbs, verbose=verbose)
>
>
>> 20170831 06:57:51| Pathname:
>> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YKG-1,chr25,6c653767ed90d29ad5254a760d8f3952.xdr
>> 20170831 06:57:51| Already done. Skipping.
>> 20170831 06:57:51|Array #1018 ('YKG-1') of 1019 on chromosome 25...done
>> 20170831 07:01:07|Genomic-signal tags:
>> 20170831 07:01:07|Reference tags: 6c653767ed90d29ad5254a760d8f3952
>> 20170831 07:01:08|Array #1019 ('YMB-1-E') of 1019 on chromosome 1...
>> 20170831 07:01:08| Pathname:
>> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YMB-1-E,chr01,6c653767ed90d29ad5254a760d8f3952.xdr
>> 20170831 07:01:08| Already done. Skipping.
>> Error: empty (zero-byte) input file
>> Execution halted
>
>
> At first I thought it was a problem of the sample, but that is not the case
> as I could process it separately. In fact, at this stage all the samples
> seem to be processed. Also strange is that the script works and exports the
> regions for subsets of the data-set.
>
> Any suggestion in how to export these would be much appreciated?
>
> Thank you,
>
> P.S: I'm using the multithreaded option of aroma
>
> # - Setup configurations
> # options(mc.cores = 10)
> print(future::availableCores())
>
> # Multithreaded
> future::plan('multiprocess')
>
> # Increase ram
> setOption(aromaSettings, 'memory/ram', 600.0)
>
> # Logs
> log <- verbose <- Arguments$getVerbose(-4, timestamp=TRUE)
>
> # Reduce decimal places to minimize space
> options(digits=4)
>
>
> --
> --
> When reporting problems on aroma.affymetrix, make sure 1) to run the latest
> version of the package, 2) to report the output of sessionInfo() and
> traceback(), and 3) to post a complete code example.
>
>
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group with website <a href="http://www.aroma-project.org/" target="_blank" rel="nofollow" onmousedown="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGgFm4_bwjITDRdF-uORwY5SRzKjA&#39;;return true;" onclick="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGgFm4_bwjITDRdF-uORwY5SRzKjA&#39;;return true;">http://www.aroma-project.org/.
> To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="ve6ut0lsBwAJ" rel="nofollow" onmousedown="this.href=&#39;javascript:&#39;;return true;" onclick="this.href=&#39;javascript:&#39;;return true;">aroma-af...@googlegroups.com
> To unsubscribe and other options, go to <a href="http://www.aroma-project.org/forum/" target="_blank" rel="nofollow" onmousedown="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2Fforum%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHDHkQUDcBfAJzIqTbjSXZqyWgITA&#39;;return true;" onclick="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2Fforum%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHDHkQUDcBfAJzIqTbjSXZqyWgITA&#39;;return true;">http://www.aroma-project.org/forum/
>
> ---
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="ve6ut0lsBwAJ" rel="nofollow" onmousedown="this.href=&#39;javascript:&#39;;return true;" onclick="this.href=&#39;javascript:&#39;;return true;">aroma-affymetr...@googlegroups.com.
> For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href=&#39;https://groups.google.com/d/optout&#39;;return true;" onclick="this.href=&#39;https://groups.google.com/d/optout&#39;;return true;">https://groups.google.com/d/optout.

--
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example.
 
 
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to [hidden email]
To unsubscribe and other options, go to http://www.aroma-project.org/forum/

---
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [hidden email].
For more options, visit https://groups.google.com/d/optout.
Reply | Threaded
Open this post in threaded view
|

Re: [aroma.affymetrix] CbsModel or writeRegions error

Henrik Bengtsson-2
Hmm... assuming you're on a *nix-like system, see if there are any
empty (zero-size) RDS files is:

ls -la cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/*.rds

The original error message suggests that there could exist such a file.

Also, you wrote "it systematically crashes on the last sample"; did
you run it multiple times and exact same error occurred at the exact
same place?  The reason why I'm asking these questions, is to rule out
certain parts of the code.  But, a traceback() would be the absolutely
most helpful information.

BTW, your mentioning 1020 samples, but the output says 1019 - probably
not important but better to make sure we're on the same page.

/Henrik

On Thu, Aug 31, 2017 at 2:42 PM, Emanuel Gonçalves
<[hidden email]> wrote:

> Hi Henrik,
>
>> does the error occur when you run:
>>
>>   fit(cbs, verbose=verbose)
>
>
> Yes, because I never got the "print(cbs)" output
>
>>
>> or did that complete successfully and you get the error while running:
>>
>>   pathname <- writeRegions(cbs, verbose=verbose)
>>
>> If you run interactively, what does traceback() output if called
>> immediately after the error occurs?
>
>
> Unfortunately, I can't run this interactively because even though the
> samples are already preprocessed it takes almost a day to load everything.
>
>> Also, see if
>>
>> fit(cbs, arrays = 1019:1020, chromosomes = c(1, 25), verbose = verbose)
>>
>> or alternatively,
>>
>>   pathname <- writeRegions(cbs, arrays = 1019:1020, chromosomes = c(1,
>> 25), verbose=verbose)
>>
>> gives the same error.
>
>
> I changed the code accordingly:
>
> fit(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)
>
> It ran without throwing any error (I tried before the same with a subset of
> ~100 samples and it also worked).
>
> It seems to be a specific problem when all the 1020 samples are ran
> together. I've placed a try/catch with the traceback().
>
> # Run CBS paired segmentation with pooled normal samples
> cbs <- CbsModel(cesN, cesN1.ref)
> out <- tryCatch(
>   {
>     fit(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)
>   },
>   error=function() {
>     message('Try/catch: error')
>     traceback()
>   },
>   warning=function() {
>     message('Try/catch: warning')
>     traceback()
>   },
>   finally={
>     message('Try/catch: finally')
>     traceback()
>   }
> )
> print(cbs)
>
> I started the script with all the samples (without "arrays=1019:1020") as
> before.
>
> What surprised me is that the script continued to the writeRegions call and
> it started to export all the samples.
>
> pathname <- writeRegions(cbs, verbose=verbose)
>
>> > # Export CBS regions
>> > pathname <- writeRegions(cbs, verbose=verbose)
>> 20170831 21:53:10|Array #1 ('201T') of 1019...
>> 20170831 21:53:10| Extracting regions from all fits...
>> 20170831 21:53:10|  Obtaining CN model fits (or fit if missing)...
>> 20170831 21:57:15|  Obtaining CN model fits (or fit if missing)...done
>> 20170831 21:57:15|  Extracting regions for chromosome #1...
>> 20170831 21:57:15|  Extracting regions for chromosome #1...done
>> 20170831 21:57:15|  Extracting regions for chromosome #2...
>> 20170831 21:57:15|  Extracting regions for chromosome #2...done
>> 20170831 21:57:15|  Extracting regions for chromosome #3...
>> 20170831 21:57:15|  Extracting regions for chromosome #3...done
>
>
> Does it mean it will export all the samples processed so far? Can I export
> everything without re-calling the fit function from CbsModel?
>
> Thanks a lot,
>
>>
>> Then you can get to traceback() a bit sooner.
>>
>> /Henrik
>>
>> PS. I assume that you already know that rerunning the commands will
>> not redo the actual analysis; already processed samples will be
>> skipped - though with 1020 samples it might still take some time.
>>
>>
>> On Thu, Aug 31, 2017 at 1:26 AM, Emanuel Gonçalves
>> <[hidden email]> wrote:
>> > Hi all,
>> >
>> > I'm processing a large set of SNP6 samples (1020) and performing CBS
>> > segmentations afterwards. It's a very lengthy process and the annoying
>> > part
>> > is that it systematically crashes on the last sample, see code and
>> > output
>> > below:
>> >
>> > # Run CBS paired segmentation with pooled normal samples
>> > cbs <- CbsModel(cesN, cesN1.ref)
>> > fit(cbs, verbose=verbose)
>> > print(cbs)
>> >
>> > # Export CBS regions
>> > pathname <- writeRegions(cbs, verbose=verbose)
>> >
>> >
>> >> 20170831 06:57:51| Pathname:
>> >>
>> >> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YKG-1,chr25,6c653767ed90d29ad5254a760d8f3952.xdr
>> >> 20170831 06:57:51| Already done. Skipping.
>> >> 20170831 06:57:51|Array #1018 ('YKG-1') of 1019 on chromosome 25...done
>> >> 20170831 07:01:07|Genomic-signal tags:
>> >> 20170831 07:01:07|Reference tags: 6c653767ed90d29ad5254a760d8f3952
>> >> 20170831 07:01:08|Array #1019 ('YMB-1-E') of 1019 on chromosome 1...
>> >> 20170831 07:01:08| Pathname:
>> >>
>> >> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YMB-1-E,chr01,6c653767ed90d29ad5254a760d8f3952.xdr
>> >> 20170831 07:01:08| Already done. Skipping.
>> >> Error: empty (zero-byte) input file
>> >> Execution halted
>> >
>> >
>> > At first I thought it was a problem of the sample, but that is not the
>> > case
>> > as I could process it separately. In fact, at this stage all the samples
>> > seem to be processed. Also strange is that the script works and exports
>> > the
>> > regions for subsets of the data-set.
>> >
>> > Any suggestion in how to export these would be much appreciated?
>> >
>> > Thank you,
>> >
>> > P.S: I'm using the multithreaded option of aroma
>> >
>> > # - Setup configurations
>> > # options(mc.cores = 10)
>> > print(future::availableCores())
>> >
>> > # Multithreaded
>> > future::plan('multiprocess')
>> >
>> > # Increase ram
>> > setOption(aromaSettings, 'memory/ram', 600.0)
>> >
>> > # Logs
>> > log <- verbose <- Arguments$getVerbose(-4, timestamp=TRUE)
>> >
>> > # Reduce decimal places to minimize space
>> > options(digits=4)
>> >
>> >
>> > --
>> > --
>> > When reporting problems on aroma.affymetrix, make sure 1) to run the
>> > latest
>> > version of the package, 2) to report the output of sessionInfo() and
>> > traceback(), and 3) to post a complete code example.
>> >
>> >
>> > You received this message because you are subscribed to the Google
>> > Groups
>> > "aroma.affymetrix" group with website http://www.aroma-project.org/.
>> > To post to this group, send email to [hidden email]
>> > To unsubscribe and other options, go to
>> > http://www.aroma-project.org/forum/
>> >
>> > ---
>> > You received this message because you are subscribed to the Google
>> > Groups
>> > "aroma.affymetrix" group.
>> > To unsubscribe from this group and stop receiving emails from it, send
>> > an
>> > email to [hidden email].
>> > For more options, visit https://groups.google.com/d/optout.
>
> --
> --
> When reporting problems on aroma.affymetrix, make sure 1) to run the latest
> version of the package, 2) to report the output of sessionInfo() and
> traceback(), and 3) to post a complete code example.
>
>
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group with website http://www.aroma-project.org/.
> To post to this group, send email to [hidden email]
> To unsubscribe and other options, go to http://www.aroma-project.org/forum/
>
> ---
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to [hidden email].
> For more options, visit https://groups.google.com/d/optout.

--
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example.


You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to [hidden email]
To unsubscribe and other options, go to http://www.aroma-project.org/forum/

---
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [hidden email].
For more options, visit https://groups.google.com/d/optout.
Reply | Threaded
Open this post in threaded view
|

Re: [aroma.affymetrix] CbsModel or writeRegions error

Emanuel Gonçalves
Hmm... assuming you're on a *nix-like system, see if there are any
empty (zero-size) RDS files is:

ls -la cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/*.rds

The original error message suggests that there could exist such a file.

Indeed, I was wondering about that. I only have *.xdr files, none had 0 bites. Also there are 25 files per sample.

Also, you wrote "it systematically crashes on the last sample"; did
you run it multiple times and exact same error occurred at the exact
same place?  

Yes, I tried three times and always with the same error.

BTW, your mentioning 1020 samples, but the output says 1019 - probably 
not important but better to make sure we're on the same page. 

In fact it's 1020 samples in total, but I excluded one sample when I got the first error (it was the very last sample). Then I re-run a sub-set of the last 20/30 samples containing the problematic ones, but I had no error and the regions were exported correctly. 

The reason why I'm asking these questions, is to rule out
certain parts of the code.  But, a traceback() would be the absolutely
most helpful information.

I only had this error  when running all the samples together in the CbsModel. Unfortunately, because it takes so long to load all the samples through CbsModel it's hard to debug. I'm running now the script with a try/catch.

Is there an alternative way to export the regions without running the CbsModel again? I have the previous script running (where only the last sample was ran in CbsModel) and it seems to be exporting all the processed samples.

Thanks a bunch,
 

/Henrik

On Thu, Aug 31, 2017 at 2:42 PM, Emanuel Gonçalves
<<a href="javascript:" target="_blank" gdf-obfuscated-mailto="0qSte597BwAJ" rel="nofollow" onmousedown="this.href=&#39;javascript:&#39;;return true;" onclick="this.href=&#39;javascript:&#39;;return true;">emanuelv...@...> wrote:

> Hi Henrik,
>
>> does the error occur when you run:
>>
>>   fit(cbs, verbose=verbose)
>
>
> Yes, because I never got the "print(cbs)" output
>
>>
>> or did that complete successfully and you get the error while running:
>>
>>   pathname <- writeRegions(cbs, verbose=verbose)
>>
>> If you run interactively, what does traceback() output if called
>> immediately after the error occurs?
>
>
> Unfortunately, I can't run this interactively because even though the
> samples are already preprocessed it takes almost a day to load everything.
>
>> Also, see if
>>
>> fit(cbs, arrays = 1019:1020, chromosomes = c(1, 25), verbose = verbose)
>>
>> or alternatively,
>>
>>   pathname <- writeRegions(cbs, arrays = 1019:1020, chromosomes = c(1,
>> 25), verbose=verbose)
>>
>> gives the same error.
>
>
> I changed the code accordingly:
>
> fit(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)
>
> It ran without throwing any error (I tried before the same with a subset of
> ~100 samples and it also worked).
>
> It seems to be a specific problem when all the 1020 samples are ran
> together. I've placed a try/catch with the traceback().
>
> # Run CBS paired segmentation with pooled normal samples
> cbs <- CbsModel(cesN, cesN1.ref)
> out <- tryCatch(
>   {
>     fit(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)
>   },
>   error=function() {
>     message('Try/catch: error')
>     traceback()
>   },
>   warning=function() {
>     message('Try/catch: warning')
>     traceback()
>   },
>   finally={
>     message('Try/catch: finally')
>     traceback()
>   }
> )
> print(cbs)
>
> I started the script with all the samples (without "arrays=1019:1020") as
> before.
>
> What surprised me is that the script continued to the writeRegions call and
> it started to export all the samples.
>
> pathname <- writeRegions(cbs, verbose=verbose)
>
>> > # Export CBS regions
>> > pathname <- writeRegions(cbs, verbose=verbose)
>> 20170831 21:53:10|Array #1 ('201T') of 1019...
>> 20170831 21:53:10| Extracting regions from all fits...
>> 20170831 21:53:10|  Obtaining CN model fits (or fit if missing)...
>> 20170831 21:57:15|  Obtaining CN model fits (or fit if missing)...done
>> 20170831 21:57:15|  Extracting regions for chromosome #1...
>> 20170831 21:57:15|  Extracting regions for chromosome #1...done
>> 20170831 21:57:15|  Extracting regions for chromosome #2...
>> 20170831 21:57:15|  Extracting regions for chromosome #2...done
>> 20170831 21:57:15|  Extracting regions for chromosome #3...
>> 20170831 21:57:15|  Extracting regions for chromosome #3...done
>
>
> Does it mean it will export all the samples processed so far? Can I export
> everything without re-calling the fit function from CbsModel?
>
> Thanks a lot,
>
>>
>> Then you can get to traceback() a bit sooner.
>>
>> /Henrik
>>
>> PS. I assume that you already know that rerunning the commands will
>> not redo the actual analysis; already processed samples will be
>> skipped - though with 1020 samples it might still take some time.
>>
>>
>> On Thu, Aug 31, 2017 at 1:26 AM, Emanuel Gonçalves
>> <[hidden email]> wrote:
>> > Hi all,
>> >
>> > I'm processing a large set of SNP6 samples (1020) and performing CBS
>> > segmentations afterwards. It's a very lengthy process and the annoying
>> > part
>> > is that it systematically crashes on the last sample, see code and
>> > output
>> > below:
>> >
>> > # Run CBS paired segmentation with pooled normal samples
>> > cbs <- CbsModel(cesN, cesN1.ref)
>> > fit(cbs, verbose=verbose)
>> > print(cbs)
>> >
>> > # Export CBS regions
>> > pathname <- writeRegions(cbs, verbose=verbose)
>> >
>> >
>> >> 20170831 06:57:51| Pathname:
>> >>
>> >> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YKG-1,chr25,6c653767ed90d29ad5254a760d8f3952.xdr
>> >> 20170831 06:57:51| Already done. Skipping.
>> >> 20170831 06:57:51|Array #1018 ('YKG-1') of 1019 on chromosome 25...done
>> >> 20170831 07:01:07|Genomic-signal tags:
>> >> 20170831 07:01:07|Reference tags: 6c653767ed90d29ad5254a760d8f3952
>> >> 20170831 07:01:08|Array #1019 ('YMB-1-E') of 1019 on chromosome 1...
>> >> 20170831 07:01:08| Pathname:
>> >>
>> >> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YMB-1-E,chr01,6c653767ed90d29ad5254a760d8f3952.xdr
>> >> 20170831 07:01:08| Already done. Skipping.
>> >> Error: empty (zero-byte) input file
>> >> Execution halted
>> >
>> >
>> > At first I thought it was a problem of the sample, but that is not the
>> > case
>> > as I could process it separately. In fact, at this stage all the samples
>> > seem to be processed. Also strange is that the script works and exports
>> > the
>> > regions for subsets of the data-set.
>> >
>> > Any suggestion in how to export these would be much appreciated?
>> >
>> > Thank you,
>> >
>> > P.S: I'm using the multithreaded option of aroma
>> >
>> > # - Setup configurations
>> > # options(mc.cores = 10)
>> > print(future::availableCores())
>> >
>> > # Multithreaded
>> > future::plan('multiprocess')
>> >
>> > # Increase ram
>> > setOption(aromaSettings, 'memory/ram', 600.0)
>> >
>> > # Logs
>> > log <- verbose <- Arguments$getVerbose(-4, timestamp=TRUE)
>> >
>> > # Reduce decimal places to minimize space
>> > options(digits=4)
>> >
>> >
>> > --
>> > --
>> > When reporting problems on aroma.affymetrix, make sure 1) to run the
>> > latest
>> > version of the package, 2) to report the output of sessionInfo() and
>> > traceback(), and 3) to post a complete code example.
>> >
>> >
>> > You received this message because you are subscribed to the Google
>> > Groups
>> > "aroma.affymetrix" group with website <a href="http://www.aroma-project.org/" target="_blank" rel="nofollow" onmousedown="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGgFm4_bwjITDRdF-uORwY5SRzKjA&#39;;return true;" onclick="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGgFm4_bwjITDRdF-uORwY5SRzKjA&#39;;return true;">http://www.aroma-project.org/.
>> > To post to this group, send email to [hidden email]
>> > To unsubscribe and other options, go to
>> > <a href="http://www.aroma-project.org/forum/" target="_blank" rel="nofollow" onmousedown="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2Fforum%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHDHkQUDcBfAJzIqTbjSXZqyWgITA&#39;;return true;" onclick="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2Fforum%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHDHkQUDcBfAJzIqTbjSXZqyWgITA&#39;;return true;">http://www.aroma-project.org/forum/
>> >
>> > ---
>> > You received this message because you are subscribed to the Google
>> > Groups
>> > "aroma.affymetrix" group.
>> > To unsubscribe from this group and stop receiving emails from it, send
>> > an
>> > email to aroma-affymetr...@googlegroups.com.
>> > For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href=&#39;https://groups.google.com/d/optout&#39;;return true;" onclick="this.href=&#39;https://groups.google.com/d/optout&#39;;return true;">https://groups.google.com/d/optout.
>
> --
> --
> When reporting problems on aroma.affymetrix, make sure 1) to run the latest
> version of the package, 2) to report the output of sessionInfo() and
> traceback(), and 3) to post a complete code example.
>
>
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group with website <a href="http://www.aroma-project.org/" target="_blank" rel="nofollow" onmousedown="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGgFm4_bwjITDRdF-uORwY5SRzKjA&#39;;return true;" onclick="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNGgFm4_bwjITDRdF-uORwY5SRzKjA&#39;;return true;">http://www.aroma-project.org/.
> To post to this group, send email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="0qSte597BwAJ" rel="nofollow" onmousedown="this.href=&#39;javascript:&#39;;return true;" onclick="this.href=&#39;javascript:&#39;;return true;">aroma-af...@googlegroups.com
> To unsubscribe and other options, go to <a href="http://www.aroma-project.org/forum/" target="_blank" rel="nofollow" onmousedown="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2Fforum%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHDHkQUDcBfAJzIqTbjSXZqyWgITA&#39;;return true;" onclick="this.href=&#39;http://www.google.com/url?q\x3dhttp%3A%2F%2Fwww.aroma-project.org%2Fforum%2F\x26sa\x3dD\x26sntz\x3d1\x26usg\x3dAFQjCNHDHkQUDcBfAJzIqTbjSXZqyWgITA&#39;;return true;">http://www.aroma-project.org/forum/
>
> ---
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to <a href="javascript:" target="_blank" gdf-obfuscated-mailto="0qSte597BwAJ" rel="nofollow" onmousedown="this.href=&#39;javascript:&#39;;return true;" onclick="this.href=&#39;javascript:&#39;;return true;">aroma-affymetr...@googlegroups.com.
> For more options, visit <a href="https://groups.google.com/d/optout" target="_blank" rel="nofollow" onmousedown="this.href=&#39;https://groups.google.com/d/optout&#39;;return true;" onclick="this.href=&#39;https://groups.google.com/d/optout&#39;;return true;">https://groups.google.com/d/optout.

--
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example.
 
 
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to [hidden email]
To unsubscribe and other options, go to http://www.aroma-project.org/forum/

---
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [hidden email].
For more options, visit https://groups.google.com/d/optout.
Reply | Threaded
Open this post in threaded view
|

Re: [aroma.affymetrix] CbsModel or writeRegions error

Henrik Bengtsson-2
I've been trying to figure where the error could show up and I have
some "poor" guesses.  In order to figure it out better, please install
the developer's version of R.utils, which I just updated:

   source('http://callr.org/install#HenrikBengtsson/R.utils@develop')

If the error occurs anywhere related to the aroma framework or some of
my underlying packages, the error message should then be much more
informative - particularly, it'll show which the problematic file is.

Unfortunately, the above would require you to rerun your script (but
it should still skip already processed files).  If it is the case that
one of the RDS files is corrupt, then a shortcut could be to do:

 path <- "cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/"
 files <- dir(path = path, pattern = "*[.]rds$", full.names = TRUE)
 for (file in files) {
    message("Reading file: ", file)
    res <- R.utils::loadObject(file)
 }

That should identify the problematic file, iff such exists.

Without understanding the problem is and / or not having access to all
your *.rds files, it's hard to say if it's safe to export the segments
or not using writeRegions(), e.g. if there's one problematic file,
then it could be that you'll be missing all segments for that
particular sample and chromosome.

Hope this help

Henrik

On Thu, Aug 31, 2017 at 3:52 PM, Emanuel Gonçalves
<[hidden email]> wrote:

>> Hmm... assuming you're on a *nix-like system, see if there are any
>> empty (zero-size) RDS files is:
>>
>> ls -la
>> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/*.rds
>>
>> The original error message suggests that there could exist such a file.
>
>
> Indeed, I was wondering about that. I only have *.xdr files, none had 0
> bites. Also there are 25 files per sample.
>
>> Also, you wrote "it systematically crashes on the last sample"; did
>> you run it multiple times and exact same error occurred at the exact
>> same place?
>
>
> Yes, I tried three times and always with the same error.
>
>> BTW, your mentioning 1020 samples, but the output says 1019 - probably
>> not important but better to make sure we're on the same page.
>
>
> In fact it's 1020 samples in total, but I excluded one sample when I got the
> first error (it was the very last sample). Then I re-run a sub-set of the
> last 20/30 samples containing the problematic ones, but I had no error and
> the regions were exported correctly.
>
>> The reason why I'm asking these questions, is to rule out
>> certain parts of the code.  But, a traceback() would be the absolutely
>> most helpful information.
>
>
> I only had this error  when running all the samples together in the
> CbsModel. Unfortunately, because it takes so long to load all the samples
> through CbsModel it's hard to debug. I'm running now the script with a
> try/catch.
>
> Is there an alternative way to export the regions without running the
> CbsModel again? I have the previous script running (where only the last
> sample was ran in CbsModel) and it seems to be exporting all the processed
> samples.
>
> Thanks a bunch,
>
>>
>>
>> /Henrik
>>
>> On Thu, Aug 31, 2017 at 2:42 PM, Emanuel Gonçalves
>> <[hidden email]> wrote:
>> > Hi Henrik,
>> >
>> >> does the error occur when you run:
>> >>
>> >>   fit(cbs, verbose=verbose)
>> >
>> >
>> > Yes, because I never got the "print(cbs)" output
>> >
>> >>
>> >> or did that complete successfully and you get the error while running:
>> >>
>> >>   pathname <- writeRegions(cbs, verbose=verbose)
>> >>
>> >> If you run interactively, what does traceback() output if called
>> >> immediately after the error occurs?
>> >
>> >
>> > Unfortunately, I can't run this interactively because even though the
>> > samples are already preprocessed it takes almost a day to load
>> > everything.
>> >
>> >> Also, see if
>> >>
>> >> fit(cbs, arrays = 1019:1020, chromosomes = c(1, 25), verbose = verbose)
>> >>
>> >> or alternatively,
>> >>
>> >>   pathname <- writeRegions(cbs, arrays = 1019:1020, chromosomes = c(1,
>> >> 25), verbose=verbose)
>> >>
>> >> gives the same error.
>> >
>> >
>> > I changed the code accordingly:
>> >
>> > fit(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)
>> >
>> > It ran without throwing any error (I tried before the same with a subset
>> > of
>> > ~100 samples and it also worked).
>> >
>> > It seems to be a specific problem when all the 1020 samples are ran
>> > together. I've placed a try/catch with the traceback().
>> >
>> > # Run CBS paired segmentation with pooled normal samples
>> > cbs <- CbsModel(cesN, cesN1.ref)
>> > out <- tryCatch(
>> >   {
>> >     fit(cbs, arrays=1019:1020, chromosomes=c(1, 25), verbose=verbose)
>> >   },
>> >   error=function() {
>> >     message('Try/catch: error')
>> >     traceback()
>> >   },
>> >   warning=function() {
>> >     message('Try/catch: warning')
>> >     traceback()
>> >   },
>> >   finally={
>> >     message('Try/catch: finally')
>> >     traceback()
>> >   }
>> > )
>> > print(cbs)
>> >
>> > I started the script with all the samples (without "arrays=1019:1020")
>> > as
>> > before.
>> >
>> > What surprised me is that the script continued to the writeRegions call
>> > and
>> > it started to export all the samples.
>> >
>> > pathname <- writeRegions(cbs, verbose=verbose)
>> >
>> >> > # Export CBS regions
>> >> > pathname <- writeRegions(cbs, verbose=verbose)
>> >> 20170831 21:53:10|Array #1 ('201T') of 1019...
>> >> 20170831 21:53:10| Extracting regions from all fits...
>> >> 20170831 21:53:10|  Obtaining CN model fits (or fit if missing)...
>> >> 20170831 21:57:15|  Obtaining CN model fits (or fit if missing)...done
>> >> 20170831 21:57:15|  Extracting regions for chromosome #1...
>> >> 20170831 21:57:15|  Extracting regions for chromosome #1...done
>> >> 20170831 21:57:15|  Extracting regions for chromosome #2...
>> >> 20170831 21:57:15|  Extracting regions for chromosome #2...done
>> >> 20170831 21:57:15|  Extracting regions for chromosome #3...
>> >> 20170831 21:57:15|  Extracting regions for chromosome #3...done
>> >
>> >
>> > Does it mean it will export all the samples processed so far? Can I
>> > export
>> > everything without re-calling the fit function from CbsModel?
>> >
>> > Thanks a lot,
>> >
>> >>
>> >> Then you can get to traceback() a bit sooner.
>> >>
>> >> /Henrik
>> >>
>> >> PS. I assume that you already know that rerunning the commands will
>> >> not redo the actual analysis; already processed samples will be
>> >> skipped - though with 1020 samples it might still take some time.
>> >>
>> >>
>> >> On Thu, Aug 31, 2017 at 1:26 AM, Emanuel Gonçalves
>> >> <[hidden email]> wrote:
>> >> > Hi all,
>> >> >
>> >> > I'm processing a large set of SNP6 samples (1020) and performing CBS
>> >> > segmentations afterwards. It's a very lengthy process and the
>> >> > annoying
>> >> > part
>> >> > is that it systematically crashes on the last sample, see code and
>> >> > output
>> >> > below:
>> >> >
>> >> > # Run CBS paired segmentation with pooled normal samples
>> >> > cbs <- CbsModel(cesN, cesN1.ref)
>> >> > fit(cbs, verbose=verbose)
>> >> > print(cbs)
>> >> >
>> >> > # Export CBS regions
>> >> > pathname <- writeRegions(cbs, verbose=verbose)
>> >> >
>> >> >
>> >> >> 20170831 06:57:51| Pathname:
>> >> >>
>> >> >>
>> >> >> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YKG-1,chr25,6c653767ed90d29ad5254a760d8f3952.xdr
>> >> >> 20170831 06:57:51| Already done. Skipping.
>> >> >> 20170831 06:57:51|Array #1018 ('YKG-1') of 1019 on chromosome
>> >> >> 25...done
>> >> >> 20170831 07:01:07|Genomic-signal tags:
>> >> >> 20170831 07:01:07|Reference tags: 6c653767ed90d29ad5254a760d8f3952
>> >> >> 20170831 07:01:08|Array #1019 ('YMB-1-E') of 1019 on chromosome 1...
>> >> >> 20170831 07:01:08| Pathname:
>> >> >>
>> >> >>
>> >> >> cbsData/gdsc_all,ACC,ra,-XY,BPN,-XY,AVG,A+B,FLN,-XY,paired/GenomeWideSNP_6/YMB-1-E,chr01,6c653767ed90d29ad5254a760d8f3952.xdr
>> >> >> 20170831 07:01:08| Already done. Skipping.
>> >> >> Error: empty (zero-byte) input file
>> >> >> Execution halted
>> >> >
>> >> >
>> >> > At first I thought it was a problem of the sample, but that is not
>> >> > the
>> >> > case
>> >> > as I could process it separately. In fact, at this stage all the
>> >> > samples
>> >> > seem to be processed. Also strange is that the script works and
>> >> > exports
>> >> > the
>> >> > regions for subsets of the data-set.
>> >> >
>> >> > Any suggestion in how to export these would be much appreciated?
>> >> >
>> >> > Thank you,
>> >> >
>> >> > P.S: I'm using the multithreaded option of aroma
>> >> >
>> >> > # - Setup configurations
>> >> > # options(mc.cores = 10)
>> >> > print(future::availableCores())
>> >> >
>> >> > # Multithreaded
>> >> > future::plan('multiprocess')
>> >> >
>> >> > # Increase ram
>> >> > setOption(aromaSettings, 'memory/ram', 600.0)
>> >> >
>> >> > # Logs
>> >> > log <- verbose <- Arguments$getVerbose(-4, timestamp=TRUE)
>> >> >
>> >> > # Reduce decimal places to minimize space
>> >> > options(digits=4)
>> >> >
>> >> >
>> >> > --
>> >> > --
>> >> > When reporting problems on aroma.affymetrix, make sure 1) to run the
>> >> > latest
>> >> > version of the package, 2) to report the output of sessionInfo() and
>> >> > traceback(), and 3) to post a complete code example.
>> >> >
>> >> >
>> >> > You received this message because you are subscribed to the Google
>> >> > Groups
>> >> > "aroma.affymetrix" group with website http://www.aroma-project.org/.
>> >> > To post to this group, send email to [hidden email]
>> >> > To unsubscribe and other options, go to
>> >> > http://www.aroma-project.org/forum/
>> >> >
>> >> > ---
>> >> > You received this message because you are subscribed to the Google
>> >> > Groups
>> >> > "aroma.affymetrix" group.
>> >> > To unsubscribe from this group and stop receiving emails from it,
>> >> > send
>> >> > an
>> >> > email to [hidden email].
>> >> > For more options, visit https://groups.google.com/d/optout.
>> >
>> > --
>> > --
>> > When reporting problems on aroma.affymetrix, make sure 1) to run the
>> > latest
>> > version of the package, 2) to report the output of sessionInfo() and
>> > traceback(), and 3) to post a complete code example.
>> >
>> >
>> > You received this message because you are subscribed to the Google
>> > Groups
>> > "aroma.affymetrix" group with website http://www.aroma-project.org/.
>> > To post to this group, send email to [hidden email]
>> > To unsubscribe and other options, go to
>> > http://www.aroma-project.org/forum/
>> >
>> > ---
>> > You received this message because you are subscribed to the Google
>> > Groups
>> > "aroma.affymetrix" group.
>> > To unsubscribe from this group and stop receiving emails from it, send
>> > an
>> > email to [hidden email].
>> > For more options, visit https://groups.google.com/d/optout.
>
> --
> --
> When reporting problems on aroma.affymetrix, make sure 1) to run the latest
> version of the package, 2) to report the output of sessionInfo() and
> traceback(), and 3) to post a complete code example.
>
>
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group with website http://www.aroma-project.org/.
> To post to this group, send email to [hidden email]
> To unsubscribe and other options, go to http://www.aroma-project.org/forum/
>
> ---
> You received this message because you are subscribed to the Google Groups
> "aroma.affymetrix" group.
> To unsubscribe from this group and stop receiving emails from it, send an
> email to [hidden email].
> For more options, visit https://groups.google.com/d/optout.

--
--
When reporting problems on aroma.affymetrix, make sure 1) to run the latest version of the package, 2) to report the output of sessionInfo() and traceback(), and 3) to post a complete code example.


You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group with website http://www.aroma-project.org/.
To post to this group, send email to [hidden email]
To unsubscribe and other options, go to http://www.aroma-project.org/forum/

---
You received this message because you are subscribed to the Google Groups "aroma.affymetrix" group.
To unsubscribe from this group and stop receiving emails from it, send an email to [hidden email].
For more options, visit https://groups.google.com/d/optout.