[Bioperl-l] data into R from Perl arrays

Roy Chaudhuri roy.chaudhuri at gmail.com
Mon Apr 2 06:33:06 EDT 2012


Alternatively you could go for a Perl-only approach using something like 
GD::Graph::Histogram.

Cheers,
Roy.

On 02/04/2012 09:17, Adam Witney wrote:
>
> The quickest way to do this specific example is probably to just save
> your R commands to a file (eg R_commands.R) and then run some kind of
> system/exec/backtick function in your perl script to invoke R,
> something like (untested):
>
> #!/usr/bin/perl use strict; use warnings; system( 'R --file
> R_commands.R' );
>
> Alternatively if you want perl and R to be able to interact and pass
> data back and forth, you could use something like Statistics::R
>
> http://search.cpan.org/~fangly/Statistics-R-0.27/lib/Statistics/R.pm
>
> HTH
>
> adam
>
> On 2 Apr 2012, at 06:21, José Ramón Blas Pastor wrote:
>
>> Hi,
>>
>> a very simple doubt, but I do not know how to manage this.
>>
>> I want to plot a histogram for all data in 'datos.txt'.
>>
>> a) by using R: datos<-scan("datos.txt") pdf("xh.pdf") hist(datos)
>> dev.off()
>>
>>
>> b) How could I invoke R inside Perl to do the same??
>> #!/usr/bin/perl open(DAT,"datos.txt"); while (<DAT>) { chomp;
>> push(@datos,$_); } #now I want a histogram of values in @datos
>>
>> Thanks!!
>>
>> JR
>>
>> -- José Ramón Blas - PhD Dept. Biochemistry - Medicine School
>> University of Castilla-La Mancha C Almansa, 14 02006 Albacete
>> (Spain)
>>
>> Phone: +34 967599200 ext. 2958
>>
>> _______________________________________________ Bioperl-l mailing
>> list Bioperl-l at lists.open-bio.org
>> http://lists.open-bio.org/mailman/listinfo/bioperl-l
>
>
> _______________________________________________ Bioperl-l mailing
> list Bioperl-l at lists.open-bio.org
> http://lists.open-bio.org/mailman/listinfo/bioperl-l



More information about the Bioperl-l mailing list