[R] t.test using RSPerl

From: Wagle, Mugdha <Mugdha.Wagle_at_STJUDE.ORG>
Date: Wed 15 Jun 2005 - 01:42:30 EST


Hi,  

I've just started using R and RSPerl. I have some code as follows:  

&R::initR("--no-save");
&R::call("t.test", (\@array1, \@array2));
 

where @array1 and @array2 are both 1-dimensional arrays in Perl having 54675 elements each. On execution the output is as follows:  

Calling R function name `t.test', # arguments: 3 1) Arg type 3
Got a reference to a value 10
Here now!2) Arg type 3
Got a reference to a value 10
Here now!Calling R

t.test(c(0, 6.24280675278087, 6.35175793656943, 5.76925805661511,
7.0789316246711, 7.4636498661157, 8.13730810691084, 8.78203131644273,
9.64502765609435, 9.95631242346133, 5.83129579495516, 6.8798700754926,
7.31814159140937.......(REST OF THE ARRAY ELEMENTS).....
4.91632461462501, 3.38099467434464,
3.91800507710569, 3.23867845216438, 3.38439026334577, 4.64918707140487,
3.23474917402449, 3.62966009445396, 3.36729582998647, 3.91999117507732
))
Performed the call, result has length 9

My question is : with other functions such as sum and log10, the actual values of the result are displayed. Here the call seems to have worked but the output is not what you get when running t.test directly on the R command prompt..  

data: data4[2] and data4[3]
t = 0.2186, df = 109.847, p-value = 0.8274 alternative hypothesis: true difference in means is not equal to 0 95 percent confidence interval:
 -3722.830 4645.723
sample estimates:
mean of x mean of y
 6185.139 5723.693  

which is what I had expected, after seeing the outputs in the case of simpler functions like sum. Could anyone please tell me how I can obtain the output I expect(i.e. the same as the command line output....giving values of t, p-value and the means)?  

Thank you very much for the help!!  

Sincerely,
Mugdha Wagle
Hartwell Center for Bioinformatics and Biotechnology, St.Jude Children's Research Hospital, Memphis TN 38105



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Jun 15 01:45:00 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:32:36 EST