/*Using R to Analyze Data in a SAS Data Set*/
proc iml;

    
/* Transfer the data to R */
    
call ExportDataSetToR("Sashelp.Class""Class");

    
/* Run a linear model using R */
    YVar = 
"Weight";
    XVar = 
"Height";
    
submit XVar YVar / R; /*Passing Parameters to R*/
        # lm=R linear model function
        Model <- lm(&YVar ~ &XVar, data=Class, 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"}];

quit;