/*Using R to Analyze Data in a SAS Data Set*/
proc iml;
    
start lm(InData, XVar, YVar);

        
/* Transfer the data to R */
        
call ExportDataSetToR(InData, "InDF");

        
/* Run a linear model using R */
        
submit  XVar YVar / R; /*Passing Parameters to R*/

            # lm=R linear model function
            Model <- lm(&YVar ~ &XVar, data=InDF, na.action=
"na.exclude")
            ParamEst <- coef(Model)                           
            Pred     <- fitted(Model)
            Resid    <- residuals(Model)
            print (Model$call)
        endsubmit;

        
/*Transfer the data from R*/
        
call ImportMatrixFromR(pe, "ParamEst");

        
/* Print the results */
        
print pe[r={"Intercept" "Height"}];
    finish;

    
store module=lm;

quit;

proc iml;

    
load module=lm;
    call lm("sashelp.class","Weight""Height");

quit;