Entering edit mode
NDowell
▴
20
@ndowell-4997
Last seen 10.2 years ago
Thanks Jim.
Yes, I had to do this a week or so ago with affxparser (and it usually
works) but I am getting the error below.
Dan's quick response that there is not Mac or Windows binary explains
the errors.
My fault for not reading to bottom of the Rsubread page.
I will have to look elsewhere.
Sorry for bothering ya'll before reading the complete web page.
Noah
> source("http://www.bioconductor.org/biocLite.R")
> biocLite("Rsubread", type= "source")
BioC_mirror: 'http://www.bioconductor.org'
Using R version 2.14, BiocInstaller version 1.2.1.
Installing package(s) 'Rsubread'
trying URL 'http://www.bioconductor.org/packages/2.9/bioc/src/contrib/
Rsubread_1.4.2.tar.gz'
Content type 'application/x-gzip' length 21848822 bytes (20.8 Mb)
opened URL
==================================================
downloaded 20.8 Mb
* installing *source* package Rsubread ...
** libs
*** arch - i386
gcc-4.2 -arch i386 -std=gnu99
-I/Library/Frameworks/R.framework/Resources/include
-I/Library/Frameworks/R.framework/Resources/include/i386
-I/usr/local/include -DMAKE_FOR_EXON -fPIC -g -O2 -c R_wrapper.c -o
R_wrapper.o
gcc-4.2 -arch i386 -std=gnu99
-I/Library/Frameworks/R.framework/Resources/include
-I/Library/Frameworks/R.framework/Resources/include/i386
-I/usr/local/include -DMAKE_FOR_EXON -fPIC -g -O2 -c aligner.c -o
aligner.o
aligner.c:301: error: expected specifier-qualifier-list before
pthread_spinlock_t
aligner.c:305: error: expected declaration specifiers or ... before
pthread_spinlock_t
aligner.c: In function run_search:
aligner.c:349: error: input_lock undeclared (first use in this
function)
aligner.c:349: error: (Each undeclared identifier is reported only
once
aligner.c:349: error: for each function it appears in.)
aligner.c:350: warning: implicit declaration of function
pthread_spin_lock
aligner.c:415: warning: implicit declaration of function
pthread_spin_unlock
aligner.c: In function run_search_thread:
aligner.c:659: error: struct gene_thread_data_transport has no
member named init_lock
aligner.c:661: error: struct gene_thread_data_transport has no
member named input_data_lock
aligner.c:661: error: too many arguments to function run_search
aligner.c: In function run_search_index:
aligner.c:719: warning: passing argument 11 of run_search makes
integer from pointer without a cast
aligner.c:719: error: too many arguments to function run_search
aligner.c:725: error: pthread_spinlock_t undeclared (first use in
this function)
aligner.c:725: error: expected ; before data_lock
aligner.c:726: error: expected ; before init_lock
aligner.c:736: error: struct gene_thread_data_transport has no
member named input_data_lock
aligner.c:736: error: data_lock undeclared (first use in this
function)
aligner.c:737: error: struct gene_thread_data_transport has no
member named init_lock
aligner.c:737: error: init_lock undeclared (first use in this
function)
aligner.c:743: warning: implicit declaration of function
pthread_spin_init
aligner.c:755: warning: implicit declaration of function
pthread_spin_destroy
aligner.c: At top level:
gene-algorithms.h:27: warning: inline function
add_gene_vote_weighted declared but never defined
gene-algorithms.h:26: warning: inline function add_gene_vote
declared but never defined
make: *** [aligner.o] Error 1
ERROR: compilation failed for package Rsubread
* removing /Library/Frameworks/R.framework/Versions/2.14/Resources/li
brary/Rsubread
The downloaded packages are in
/private/var/folders/14/14lHU6sYHXWDo48Jc3Lq++lgIs+/-Tmp-/Rtm
p9gion7/downloaded_packages
Warning message:
In install.packages(pkgs = pkgs, repos = repos, ...) :
installation of package Rsubread had non-zero exit status
On Dec 5, 2011, at 3:15 PM, James W. MacDonald wrote:
> Hi Noah,
>
> On 12/5/2011 4:09 PM, Noah Dowell wrote:
>> Hello All,
>>
>> I am trying to use the Rsubread package which is supposedly
available but I am getting the following warning message.
>
> It is only available as a source package. If you are set up to build
packages from source, you can simply add the argument type="source" to
your biocLite() call.
>
> Otherwise,
>
> http://cran.r-project.org/bin/macosx/RMacOSX-FAQ.html#Installation-
of-source-packages
>
> Best,
>
> Jim
>
>
>> Thanks for your help.
>>
>> Best,
>>
>> Noah
>>
>>
>>> source("http://www.bioconductor.org/biocLite.R")
>> BiocInstaller version 1.2.1, ?biocLite for help
>>> biocLite("Rsubread")
>> BioC_mirror: 'http://www.bioconductor.org'
>> Using R version 2.14, BiocInstaller version 1.2.1.
>> Installing package(s) 'Rsubread'
>> Warning message:
>> In getDependencies(pkgs, dependencies, available, lib) :
>> package Rsubread is not available (for R version 2.14.0)
>>
>>
>>> sessionInfo()
>> R version 2.14.0 (2011-10-31)
>> Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
>>
>> locale:
>> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>>
>> attached base packages:
>> [1] stats graphics grDevices utils datasets methods
base
>>
>> other attached packages:
>> [1] BiocInstaller_1.2.1
>>
>> loaded via a namespace (and not attached):
>> [1] tools_2.14.0
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor@r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
http://news.gmane.org/gmane.science.biology.informatics.conductor
>
> --
> James W. MacDonald, M.S.
> Biostatistician
> Douglas Lab
> University of Michigan
> Department of Human Genetics
> 5912 Buhl
> 1241 E. Catherine St.
> Ann Arbor MI 48109-5618
> 734-615-7826
>
> **********************************************************
> Electronic Mail is not secure, may not be read every day, and should
not be used for urgent or sensitive issues
[[alternative HTML version deleted]]
Dear Wei,
i think you might consider the following three reasons why it would be interesting to have Rsubread available in Mac OS X systems:
i've used Rsubread to teach my students about mapping reads without leaving the R console by using a small sample of reads (e.g. 1e5, 1e6), if eventually one makes a teaching course where students come with laptops, a significant fraction of them come now with mac. this could apply to all the BioC courses that take place throughout the world :)
apple not only sells laptops but also high-performance computing solutions. those users would also benefit of having Rsubread running in Mac OS X.
Mac OS X runs on top of a unix system, so it should not be all too difficult to have and maintain ANSI C code running in both, a linux and a Mac OS X system.
cheers, robert.
Dear Wei,
Thanks for your response. I was going to use Rsubread with a small(ish) data set on a desktop with a little more memory than the standard laptop. I like the approach behind Rsubread and I wanted to do some comparisons with other aligners on small data sets before making a choice.
Of course I understand as a maintainer you have to make some careful choices where and what to put your effort into.
I should have a server up and running soon (at my new institute; not UCLA) and will consider using Rsubread at that time.
Best, Noah
FWIW, I've had pretty good results with some older 36-cycle Illumina RNAseq data using Rsubread, just aligning it on my laptop (granted it's a dual-core i7 with 16GB of RAM and 2 SSDs). For some of the RNAseq data from CD34+CD38- cells and neutrophils that Andrew Smith was kind enough to send, it took about 30 minutes, and I was able to do other things at the same time. Again, that might have something to do with my using 7.4GB for the index, having another 8.6GB for other processes, and running Linux. But mostly the RAM, I think.
YMMV... Rsubread has some terrific features, I am torn between using it all the time, using BowTie/TopHat/Cufflinks/cummeRbund all the time, or both. Last time I just used both, but with our RNAseq/QC pipeline, if I want to re-align against hg18 I either have to use bedtools and liftOver, re- run the pipeline, or use Rsubread. The latter is more convenient, even if the former does have some more tools for removing PCR dupes, looking at alternative splicing, etc.
Either one runs fine on our servers, of course, and both go a lot faster with 24 cores and 48GB of RAM :-)
One other thing Noah, do you have Xcode and the R tools for Mac OSX installed? I'm using my wife's Macbook Air at this particular instant, and I'm tempted to try compiling Rsubread, but I don't really have time to debug something for fun right now. Still -- it would be interesting to see if that resolved your issues. It's not like laptops are being produced with less RAM and smaller processors as time goes by, and there are old-but-still-very-useful RNAseq datasets out there. And Rsubread is FAST so it's by far the most sensible choice for a laptop.
Dear Tim, Robert and Noah,
I think we will eventually make Rsubread available on Mac OS X system (but may not on Windows). Hopefully it will become available on the next bioc release. But Yang Liao, who wrote the entire C program for the read alignment, has to submit his Ph.D thesis February next year, and also we are now busy with publishing Rsubread. But I'm quite sure it should be available sometime next year.
Rsubread uses the spinlock library to enable multithreaded running of the read alignment, which is not supported on Mac. We'll have to change it to the mutex library which is supported by both Mac and Linux for portability (however, we found spinlock is more efficient than mutex:-( ).
The align() function in Rsubread can map junction reads from the RNA- seq data in addition to the exonic reads, although it won't give you the information of where the junction locations are (it gives the mapping location of the junction read in one of the exons it spans). The subjunc() function in Rsubread package can however be used to find the exact exon junction locations in the reference genome.
Cheers,
Wei
TIm,
Thanks for providing your experiences on using different aligners. I used R/Bioconductor for tiling and expression arrays as a grad student but now I am using next generation sequencing as a postdoc so I am getting up to speed with some new packages like Rsubread.
Yes, I have the Xcode developer tools installed and have installed from source recently for other packages with no errors.
Best, Noah
Hi all,
Wei can correct me if I am wrong, but I believe it's not as simple as building Rsubread from source on Mac (or Windows). If it were, Bioconductor would provide a binary package on those platforms.
I believe the problem is that certain constructs (e.g. phtread) are used in such a way that is not portable to Mac OS X or Windows,
Dan
Hi Dan,
Yes, that's right. pthread is a construct which is not portable between linux and Mac. But there might be other such constructs as well. We will need to have a close look at them.
Thanks, Wei