Date:         Wed, 22 Oct 2003 10:41:39 -0700
Reply-To:     Dale McLerran <stringplayer_2@YAHOO.COM>
Sender:       "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU>
From:         Dale McLerran <stringplayer_2@YAHOO.COM>
Subject:      Re: SAS/IML
Comments: To: sonia afroz <afroz_sonia@hotmail.com>
In-Reply-To:  <Sea1-F52IIzQtO9vXLs000237f2@hotmail.com>
Content-Type: text/plain; charset=us-ascii

Sonia,

It would be most helpful if you presented the log file, rather
than just the SAS code.  The log should indicate any problems
opening the data (file not found type messages) as well as
problems with variable identification (variables named on the
READ statement are not on the input dataset).

I might mention that I took a quick look at your regression
code.  I would note that you have not included an intercept
in your regression.  Is that what you mean to do, or is it a
coding problem?

Dale


--- sonia afroz <afroz_sonia@hotmail.com> wrote:
> Hi,
>
>
> I am running this proc IML program and not getting any outputs. Seems
> like,
> I am not being able to read the data. Could you please check and see
> what I
> am doing wrong here?
> Thanks,
> sonia
>
> options nocenter linesize=256;
> libname ssd 'U:\Users\sa81\My Documents\ORSdata';
>
> proc iml;
> use ssd.inmoodata;
>
> read all var {rp12} into y ;
> read all var { x1 x2 x3 rm12 rm36} into x;
>
>    start reg;
> reset print;
>       n=nrow(x);                      /* number of observations */
>       k=ncol(x);                         /* number of variables */
>       xpx=x`*x;                               /* cross-products */
>       xpy=x`*y;
>       xpxi=inv(xpx);                   /* inverse crossproducts */
>       b=xpxi*xpy;                        /* parameter estimates */
>       yhat=x*b;                             /* predicted values */
>       resid=y-yhat;                                /* residuals */
>       sse=resid`*resid;                /* sum of squared errors */
>       dfe=n-k;                      /* degrees of freedom error */
>       mse=sse/dfe;                        /* mean squared error */
>       rmse=sqrt(mse);                /* root mean squared error */
>       covb=xpxi#mse;                 /* covariance of estimates */
>       stdb=sqrt(vecdiag(covb));              /* standard errors */
>       t=b/stdb;                        /* ttest for estimates=0 */
>       probt=1-probf(t#t,1,dfe);     /* significance probability */
>       print name b stdb t probt;
>
> _________________________________________________________________
> Concerned that messages may bounce because your Hotmail account has
> exceeded
> its 2MB storage limit? Get Hotmail Extra Storage!
> http://join.msn.com/?PAGE=features/es
>


=====
---------------------------------------
Dale McLerran
Fred Hutchinson Cancer Research Center
mailto: dmclerra@fhcrc.org
Ph:  (206) 667-2926
Fax: (206) 667-5977
---------------------------------------

__________________________________
Do you Yahoo!?
The New Yahoo! Shopping - with improved product search
http://shopping.yahoo.com