Here is the code I am using (which I know is not quite right).  It  
produces the output below.  The data frame summary seems to print  
correctly, but the call to 'lm' says that I do not have a 'terms  
component'.   I think I am confused about whether objects reside in R  
or Perl although my code assumes I can create a Perl object (in this  
case df) that can be passed back to R as a data frame.  Can someone  
clarify this for me?

--------------
R::initR("--silent");
&R::library("RSPerl");

my @df = &R::callWithNames( "data.frame", { 'date', $r_date[0], 'r1',  
$e_ret[0], 'r2', $e_ret[1] } );

print "Printing summary of df\n";
&R::call( "print", @df);

print "Computing linear model parameters\n";
&R::callWithNames( "lm", { 'r2 ~ r1' , 'data', [EMAIL PROTECTED] }  );

--------------

          date                    r2                    r1
1  2006-05-01   -0.0573445954863402   -0.0112231851115216
2  2006-05-02    0.0150507525376269    0.0290229885057473
3  2006-05-03 -0.000275862068965462  -0.00670203853672164
4  2006-05-04    0.0117273730684326 -0.000140567894293015
5  2006-05-05    0.0278876312559662    0.0106846618866864
6  2006-05-08    -0.054792703150912                     0
7  2006-05-09  -0.00357919854024857   -0.0119627208234803
8  2006-05-10  0.000986054373855577  -0.00605378009291858
9  2006-05-11  -0.00830284266816779   -0.0347025495750707
10 2006-05-12  -0.00744997871434661  -0.00660308143800444
11 2006-05-15  -0.00986489384516403   0.00132939438700153
12 2006-05-16   0.00678651360912568   -0.0414515415252987
13 2006-05-17   -0.0254571531014699   0.00430901815943369
14 2006-05-18   -0.0300956585724798   -0.0318725099601594
15 2006-05-19   0.00485547378802812      0.02105096549541
16 2006-05-22   -0.0130615326538316   -0.0175166640830879
17 2006-05-23    0.0155293757649939   -0.0036289050173557
18 2006-05-24    0.0207909604519773   0.00300870942201116
19 2006-05-25   0.00715814331045678    0.0156299336911903
20 2006-05-26   0.00600820633059805   -0.0121249805689414
Printing summary of df
          date                    r2                    r1
1  2006-05-01   -0.0573445954863402   -0.0112231851115216
2  2006-05-02    0.0150507525376269    0.0290229885057473
3  2006-05-03 -0.000275862068965462  -0.00670203853672164
4  2006-05-04    0.0117273730684326 -0.000140567894293015
5  2006-05-05    0.0278876312559662    0.0106846618866864
6  2006-05-08    -0.054792703150912                     0
7  2006-05-09  -0.00357919854024857   -0.0119627208234803
8  2006-05-10  0.000986054373855577  -0.00605378009291858
9  2006-05-11  -0.00830284266816779   -0.0347025495750707
10 2006-05-12  -0.00744997871434661  -0.00660308143800444
11 2006-05-15  -0.00986489384516403   0.00132939438700153
12 2006-05-16   0.00678651360912568   -0.0414515415252987
13 2006-05-17   -0.0254571531014699   0.00430901815943369
14 2006-05-18   -0.0300956585724798   -0.0318725099601594
15 2006-05-19   0.00485547378802812      0.02105096549541
16 2006-05-22   -0.0130615326538316   -0.0175166640830879
17 2006-05-23    0.0155293757649939   -0.0036289050173557
18 2006-05-24    0.0207909604519773   0.00300870942201116
19 2006-05-25   0.00715814331045678    0.0156299336911903
20 2006-05-26   0.00600820633059805   -0.0121249805689414
Computing linear model parameters
Error in terms.formula(formula, data = data) :
         argument is not a valid model


Begin forwarded message:

> From: Adam Wizon <[EMAIL PROTECTED]>
> Date: June 28, 2006 10:58:35 AM EDT
> To: [email protected]
> Subject: [R-SIG-Mac] Passing data to lm using RSPerl
>
> I have successfully installed RSPerl.  I am retrieving data from a
> database in Perl and trying to run a linear model.  I have read the
> documentation, but I do not quite understand how to set up the call
> to 'lm' inside Perl.  Can someone provide an example I can work
> from?  Thanks.
>
> Adam Wizon
>
> _______________________________________________
> R-SIG-Mac mailing list
> [email protected]
> https://stat.ethz.ch/mailman/listinfo/r-sig-mac


        [[alternative HTML version deleted]]

_______________________________________________
R-SIG-Mac mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-mac

Reply via email to