topic Create SAS data set from procedure reg in SAS Procedures

Create SAS data set from procedure reg

Pankaj_23 — Thu, 08 Sep 2011 05:40:30 GMT

Dear All,

Create SAS data set contain output of standard errror of slope, standard errror of intercept. output generated by proc glm.

Following is SAS code:

Data inputdata;

input x y;

datalines;

1 0.9450 0.060958

2 1.8580 0.143329

4 3.7480 0.279329

8 7.4240 0.597336

10 9.1880 0.820075

20 11.7960 8.247110

30 12.6310 17.404420

40 12.7820 27.297152

;

Run;

proc reg data=inputdata ;

model y = x;

run;

After execution of above sas code, we get sas output contain slope, intercept, se of slope, se of intercept.

how to save sas output (create SAS data set) from proc reg?

Please anyone help me?

Thanks in advance.

Best Reagrads,

Pankaj Bhangale

Create SAS data set from procedure reg

Cynthia_sas — Thu, 08 Sep 2011 12:02:42 GMT

Hi,

You can either look in the documentation for PROC REG to see how it might use the OUTPUT statement or the OUT= option to create a SAS dataset from your procedure results OR you can read about the

ODS OUTPUT statement to create a SAS dataset from PROC REG output objects.

Cynthia

Create SAS data set from procedure reg

SteveDenham — Thu, 08 Sep 2011 12:03:11 GMT

Add an ods output statement.

proc reg data=inputdata ;

model y = x;

ods output parameterestimates=parms;

run;

The dataset parms will have what you request.

Steve Denham

Create SAS data set from procedure reg

Pankaj_23 — Fri, 09 Sep 2011 06:34:38 GMT

Thank you very much for this quick reply.

It is helpful for me in another application.

Best Regards,

Pankaj Bhangale

Create SAS data set from procedure reg

Pankaj_23 — Fri, 09 Sep 2011 06:37:14 GMT

Thank you very much for this quick reply.Your answer is correct.

It is helpful for me in another application.

Best Regards,

Pankaj Bhangale

Create SAS data set from procedure reg

sunilpusarla — Tue, 13 Sep 2011 23:23:54 GMT

Hi,

If you want to know the dataset names of ods output, use ods trace on/listing before the stat proc. Then you can pick what you want, from output window.

Thanks,

Sunil

Re: Create SAS data set from procedure reg

Pankaj_23 — Mon, 19 Sep 2011 11:55:50 GMT

Dear Steve Denham,

Also I want required sum of squares.

Let Y is response variable & X is independent variable.

How to create SAS data set contain output of sum of squares like Sxx, Sxy, Syy & SST (Total sum of square) through proc reg procedure?

Please anyone help me.

Thanks in advance.

Best Regards,

Pankaj Bhangale

Re: Create SAS data set from procedure reg

art297 — Mon, 19 Sep 2011 12:30:51 GMT

I believe that one way to get the calculations you are looking for is to specify the outsscp option. Take a look at: http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/viewer.htm#statug_reg_sect007.htm#statug.reg.procoutsscp

Re: Create SAS data set from procedure reg

Ksharp — Wed, 21 Sep 2011 03:34:08 GMT

ods output anova=inputanova;

proc reg data=sashelp.class ;

model weight=height;

run;

The dataset inputanova will contain what you required.

Ksharp

Re: Create SAS data set from procedure reg

Pankaj_23 — Wed, 21 Sep 2011 05:09:25 GMT

Dear Ksharp,

But the dataset inputanova does not contain Sxx & Sxy.

Slope_Hat =(Sxy/Sxx)

Slope ~ N(Slope_Hat, Sigma^2/Sxx)

I want Sxy & Sxx for above application.

Best Regards

Pankaj Bhangale

Re: Create SAS data set from procedure reg

Ksharp — Wed, 21 Sep 2011 09:22:25 GMT

Oh.....

For corrected sum of squares for x or y, You can use proc means with css option to get.

But for corrected cross product,I am afraid that you need to calculated by hand.

proc means noprint data=sashelp.class;
 var weight height;
 output out=Sxx(drop=_:) css=weight_css height_css mean=weight_mean height_mean;
run;



data class;
 set sashelp.class;
 if _n_ eq 1 then set Sxx;
run;
proc sql ;
 create table Sxy_Sxx as
  select weight_css,height_css ,
       sum((weight-weight_mean)*(height-height_mean))as Sxy
    from class;
quit;

Ksharp

Re: Create SAS data set from procedure reg

Pankaj_23 — Wed, 21 Sep 2011 11:49:10 GMT

Thank Ksharp for quick reply. you are great.

Best Regards,

Pankaj Bhangale

Re: Create SAS data set from procedure reg

art297 — Wed, 21 Sep 2011 13:05:42 GMT

KSharp and Pankaj, take a look at the reference I provided earlier in this thread, particularly regarding what is included with the outsscp option. The documentation states:

OUTSSCP=

outputs a data set that contains sums of squares and crossproducts

Re: Create SAS data set from procedure reg

Ksharp — Thu, 22 Sep 2011 02:56:55 GMT

Yes.

Art is right. I should check documentation firstly. Sorry! Pankaj Bhangale ,You should use outsscp= option to get that.

proc reg data=fitness outsscp=sscp;

var Oxygen RunTime Age Weight RestPulse RunPulse MaxPulse;

proc print data=sscp;

run;

Ksharp

Re: Create SAS data set from procedure reg

Ksharp — Thu, 22 Sep 2011 03:49:21 GMT

HaHa.

Art297.

I just found that outsscp= option gave the Uncorrected Sum of Square, not Corrected Sum of Square and

The cross product is also uncorrected .So Maybe My code is what Op want.

I am so curious that why sas can not give us Corrected Sum of Square in proc reg.

Ksharp

Re: Create SAS data set from procedure reg

Ksharp — Thu, 22 Sep 2011 03:55:40 GMT

Yes. Art.

I am right. I revise it again.

Ksharp

Re: Create SAS data set from procedure reg

Pankaj_23 — Thu, 22 Sep 2011 05:04:02 GMT

Ksharp you are right. I got correct solution by using your SAS code and also matched with manual calculation.

I Happy....

Ok but your SAS Code based on Proc Means.

Again question is arise, I thought Slope & Intercept obtained from proc reg.

Slope=(Sxy/Sxx) where Sxx=Corrected sum of Square.

Why sas can not give us Corrected Sum of Square in proc reg?

Thanks,

Best Regards,

Pankaj Bhangale

Re: Create SAS data set from procedure reg

Ksharp — Thu, 22 Sep 2011 08:06:48 GMT

Hi.

Maybe SAS calculated estimated coefficient is based on Uncorrected Sum of Square not Corrected.

It is also correct .But you use Corrected Sum of Square , you are also right. They are actually the same thing,but formula is different.

So I also hope SAS can offer CSS in proc reg in the future.

Ksharp