Re: Question on SAS Code!

--On 10/19/07 1:25 AM +0000 yahoo.com">mrfurley15@yahoo.com wrote:

> I am lost. I am getting all kinds of errors.

[...]
> 27 GOptions Device=gif Transparency NoBorder
> --------
> 180
> ERROR 180-322: Statement is not valid or it is used out of proper
> order.
[...]
> 113 Proc Gplot Data=Boards;
> ERROR: Procedure GPLOT not found.
> 114 Plot (X1 X2 X3 X4)*(X1 X2 X3 X4)=1 / Name="Boards";
> 115 Symbol1 C=Black V=Dot H=1;
> 116 Run;
>
> NOTE: The SAS System stopped processing this step because of errors.
[...]
> 132 Proc IML;
> ERROR: Procedure IML not found.
> 133 Reset NoLog; /* send output to the listing file */
[...]

To me, these errors look like what I would expect if SAS/GRAPH and SAS/IML were
not licensed or not installed. Are you sure these products are present and
licensed?


Dave Scocca

On Oct 18, 9:38 pm, d...@SCOCCA.ORG (Dave Scocca) wrote:
> --On 10/19/07 1:25 AM +0000 mrfurle...@yahoo.com wrote:
>
> > I am lost. I am getting all kinds of errors.
>
>
> To me, these errors look like what I would expect if SAS/GRAPH and SAS/IML were
> not licensed or not installed. Are you sure these products are present and
> licensed?
>
> Dave Scocca

Dave is correct. You do not have the proper SAS procedures licensed.

--
Paige Miller
paige\dot\miller \at\ kodak\dot\com

On Oct 19, 8:07 am, Paige Miller <paige.mil...@kodak.com> wrote:
> On Oct 18, 9:38 pm, d...@SCOCCA.ORG (Dave Scocca) wrote:
>
> > --On 10/19/07 1:25 AM +0000 mrfurle...@yahoo.com wrote:
>
> > > I am lost. I am getting all kinds of errors.
>
> > To me, these errors look like what I would expect if SAS/GRAPH and SAS/IML were
> > not licensed or not installed. Are you sure these products are present and
> > licensed?
>
> > Dave Scocca
>
> Dave is correct. You do not have the proper SAS procedures licensed.
>
> --
> Paige Miller
> paige\dot\miller \at\ kodak\dot\com

I see, thanks. I'll check it out on a computer that I know has that.
But what do I need to change here in the code to make the oliverdata
work? I know it isn't too much, but I am stumped! Please can you
help?

On Oct 21, 11:41 pm, "mrfurle...@yahoo.com" <mrfurle...@yahoo.com>
wrote:
> On Oct 19, 8:07 am, Paige Miller <paige.mil...@kodak.com> wrote:

> I see, thanks. I'll check it out on a computer that I know has that.
> But what do I need to change here in the code to make the oliverdata
> work? I know it isn't too much, but I am stumped! Please can you
> help?

You will need to attempt to modify the code yourself. I don't think it
is too much to ask for someone to take an example with four variables
and modify it so that it now runs with 8 variables. If you run into
errors, post appropriate portions of the SASLOG here which contain the
errors.

--
Paige Miller
paige\dot\miller \at\ kodak\dot\com

On Oct 22, 8:03 am, Paige Miller <paige.mil...@kodak.com> wrote:
> On Oct 21, 11:41 pm, "mrfurle...@yahoo.com" <mrfurle...@yahoo.com>
> wrote:
>
> > On Oct 19, 8:07 am, Paige Miller <paige.mil...@kodak.com> wrote:
> > I see, thanks. I'll check it out on a computer that I know has that.
> > But what do I need to change here in the code to make the oliverdata
> > work? I know it isn't too much, but I am stumped! Please can you
> > help?
>
> You will need to attempt to modify the code yourself. I don't think it
> is too much to ask for someone to take an example with four variables
> and modify it so that it now runs with 8 variables. If you run into
> errors, post appropriate portions of the SASLOG here which contain the
> errors.
>
> --
> Paige Miller
> paige\dot\miller \at\ kodak\dot\com

OK thanks...can you give me a hint on what should be changed besides
the variables.?

On Oct 22, 10:13 am, "mrfurle...@yahoo.com" <mrfurle...@yahoo.com>
wrote:
> On Oct 22, 8:03 am, Paige Miller <paige.mil...@kodak.com> wrote:

> OK thanks...can you give me a hint on what should be changed besides
> the variables.?

I believe you need to incorporate your data into the program
somewhere.

--
Paige Miller
paige\dot\miller \at\ kodak\dot\com

Ok Paige, here is my SAS code:


************************************************** **;
* olivetst.sas
************************************************** **;

Options PS=55 LS=80 PageNo=1 NoDate MPrint
FORMCHAR='|----|+|---+=|-/\<>*';

************************************************** **;
* The following section generates HTML pages for *;
* web presentation. Delete the ODS and GOptions *;
* statements if you do not want the web pages. *;
************************************************** **;

ODS HTML body= "olivetst-body.html"
contents="olivetst-contents.html"
frame= "olivetst-frame.html"
page= "olivetst-page.html"
headtext="<title>Multivariate Normality</title>"
anchor="olivetst";

GOptions Device=gif Transparency NoBorder
VSize=6 in HSize=6 in
HText=1 FText=Swiss HTitle=1 FTitle=Swiss;

************************************************** **;
* The Data Analysis Begins Here.
************************************************** **;
Title1 "Classification Analysis of Olive Data";

Data olivetst;
Infile palmitic palmitoleic stearic oleic
linoleic eicosanoic linolenic eicosenoic;
cards;

1.454 0.183 0.196 7.057 1.014 0.0149 0.046 0.0067
1.347 0.194 0.197 7.277 0.895 0.0122 0.046 0.0045
1.364 0.204 0.225 6.929 1.084 0.0082 0.05 0.0041
1.41 0.199 0.216 7.13 0.955 0.0082 0.048 0.0067
1.384 0.178 0.208 7.105 0.999 0.0182 0.067 0.0135
1.412 0.185 0.217 6.842 1.203 0.03 0.072 0.0245
1.41 0.232 0.28 6.715 1.233 0.0245 0.06 0.011
1.509 0.209 0.257 6.647 1.24 0.0667 0.062 0.0201
1.317 0.197 0.256 7.036 1.067 0.0546 0.06 0.009
1.286 0.192 0.203 7.132 1.053 0.0447 0.065 0.0164
1.273 0.191 0.202 6.862 1.303 0.0737 0.07 0.0164
1.463 0.183 0.183 6.747 1.307 0.0366 0.06 0.011
1.399 0.187 0.191 6.861 1.233 0.0447 0.06 0.0055
1.413 0.193 0.208 6.875 1.202 0.0201 0.06 0.006
1.369 0.206 0.203 6.953 1.168 0.0331 0.05 0.005
1.488 0.172 0.17 6.92 1.144 0.0404 0.054 0.0041
1.323 0.16 0.205 6.911 1.298 0.011 0.05 0.0055
1.311 0.166 0.17 6.902 1.312 0.0603 0.069 0.0164
1.286 0.163 0.183 7.04 1.23 0.0182 0.057 0.0033
1.38 0.173 0.188 7.038 1.139 0.0222 0.044 0.0041
1.394 0.164 0.223 7.086 1.042 0.011 0.043 0.01
1.324 0.174 0.198 6.863 1.289 0.0366 0.07 0.0082
1.29 0.157 0.192 7 1.263 0.0135 0.051 0.0067
1.361 0.163 0.196 6.888 1.273 0.0404 0.058 0.011
1.387 0.182 0.242 6.913 1.101 0.0815 0.068 0.0201
1.369 0.18 0.181 7 1.13 0.0494 0.045 0.011
1.303 0.165 0.175 7.025 1.243 0.0222 0.041 0.005
1.346 0.16 0.169 7.072 1.151 0.0494 0.048 0.0045
1.369 0.171 0.184 6.937 1.246 0.0201 0.048 0.0045
1.305 0.172 0.169 7.004 1.26 0.0164 0.05 0.003
1.351 0.179 0.186 6.935 1.243 0.0366 0.05 0.0067
1.283 0.151 0.182 7 1.271 0.0546 0.052 0.0082
1.449 0.175 0.198 6.883 1.162 0.0546 0.07 0.009
1.31 0.18 0.183 7.054 1.202 0.0135 0.032 0.0033
1.36 0.163 0.176 6.901 1.28 0.0164 0.065 0.0149
1.3 0.187 0.196 6.92 1.253 0.0603 0.076 0.0122
1.368 0.171 0.218 7.01 1.057 0.0603 0.054 0.0135
1.207 0.151 0.156 7.159 1.234 0.0149 0.051 0.0041
1.348 0.154 0.183 6.917 1.277 0.1215 0.056 0.005
1.632 0.186 0.229 7.261 0.827 0.03 0.056 0.0074
1.301 0.156 0.207 7.003 1.229 0.0603 0.048 0.0041
1.226 0.181 0.213 6.961 1.23 0.1099 0.074 0.0135
1.201 0.168 0.19 7.1 1.216 0.0737 0.064 0.005
1.297 0.153 0.177 7.004 1.26 0.0331 0.06 0.005
1.248 0.163 0.158 7.103 1.222 0.0222 0.06 0.0041
1.335 0.159 0.197 6.974 1.22 0.0366 0.06 0.0055
1.219 0.167 0.171 7.087 1.254 0.0331 0.05 0.005
1.318 0.179 0.177 7.03 1.194 0.0331 0.042 0.0122
1.264 0.167 0.166 7.13 1.187 0.009 0.052 0.0033
1.201 0.175 0.201 7.129 1.193 0.0366 0.049 0.0045
1.252 0.18 0.181 7.055 1.214 0.0222 0.059 0.0447
1.273 0.182 0.209 6.965 1.191 0.0737 0.074 0.01
1.351 0.179 0.17 7.034 1.154 0.0331 0.066 0.0027
1.336 0.155 0.212 7.103 1.086 0.0271 0.055 0.0074
1.499 0.201 0.182 6.803 1.204 0.0201 0.056 0.011
1.425 0.198 0.193 7.032 1.041 0.0222 0.052 0.0055
1.358 0.204 0.227 6.962 1.109 0.0603 0.065 0.03
1.346 0.181 0.257 7.147 0.933 0.0546 0.06 0.0366
1.392 0.186 0.256 6.732 1.278 0.2003 0.064 0.0182
1.311 0.166 0.222 7.006 1.147 0.0603 0.08 0.0149
1.314 0.171 0.229 6.923 1.198 0.1099 0.076 0.0667
1.409 0.2 0.207 6.842 1.224 0.0222 0.06 0.0149
1.342 0.174 0.221 6.993 1.147 0.0366 0.064 0.01
1.387 0.182 0.206 7.1 1.02 0.03 0.054 0.0055
1.413 0.202 0.205 6.92 1.165 0.0366 0.046 0.0037
1.43 0.209 0.225 6.8 1.2 0.0245 0.059 0.0149
1.336 0.185 0.223 6.956 1.155 0.2704 0.073 0.005
1.372 0.2 0.2 6.916 1.189 0.0271 0.05 0.009
1.33 0.157 0.228 7.055 1.108 0.0667 0.055 0.0122
1.412 0.207 0.208 6.822 1.239 0.0366 0.051 0.0164
1.321 0.209 0.217 6.948 1.178 0.0667 0.062 0.01
1.401 0.2 0.217 6.98 1.073 0.0546 0.068 0.0082
1.401 0.214 0.217 6.734 1.293 0.0815 0.069 0.0149
1.457 0.168 0.242 6.724 1.266 0.2214 0.059 0.0201
1.451 0.199 0.221 6.835 1.177 0.0404 0.051 0.0182
1.438 0.206 0.248 6.806 1.183 0.03 0.057 0.0164
1.462 0.204 0.237 6.644 1.309 0.0667 0.054 0.0164
1.529 0.215 0.203 6.602 1.31 0.09 0.069 0.0149
1.51 0.189 0.245 6.752 1.188 0.0366 0.052 0.0164
1.437 0.222 0.184 6.803 1.24 0.0737 0.056 0.005
1.327 0.129 0.247 7.024 1.157 0.0447 0.056 0.009
1.438 0.172 0.252 6.63 1.38 0.0546 0.064 0.011
1.447 0.176 0.189 6.849 1.18 0.0667 0.064 0.0135
1.355 0.144 0.214 6.972 1.198 0.0271 0.06 0.011
1.369 0.156 0.241 6.89 1.209 0.0667 0.063 0.0201
1.471 0.188 0.276 6.697 1.269 0.03 0.051 0.005
1.456 0.179 0.24 6.738 1.267 0.0603 0.065 0.0041
1.314 0.14 0.207 7.02 1.22 0.0164 0.059 0.0033
1.408 0.176 0.192 6.909 1.195 0.09 0.05 0.0122
1.397 0.172 0.191 7.107 1.018 0.0366 0.05 0.0182
1.413 0.191 0.186 6.937 1.18 0.0222 0.046 0.0037
1.539 0.194 0.213 6.764 1.178 0.0447 0.058 0.005
1.304 0.159 0.234 7.019 1.174 0.0447 0.053 0.0067
1.341 0.16 0.231 7.033 1.069 0.0546 0.067 0.0271
1.508 0.208 0.249 6.641 1.311 0.0122 0.043 0.0074
1.515 0.226 0.257 6.595 1.287 0.0603 0.063 0.005
1.262 0.165 0.235 7.12 1.113 0.0245 0.051 0.0082
1.307 0.197 0.238 7.003 1.144 0.0404 0.05 0.011
1.294 0.159 0.253 7.009 1.19 0.0201 0.052 0.0037
1.46 0.187 0.215 6.843 1.172 0.0331 0.056 0.0245
1.476 0.187 0.203 6.837 1.197 0.0366 0.048 0.009
1.482 0.178 0.197 6.814 1.201 0.0546 0.064 0.011
1.388 0.176 0.185 7.008 1.111 0.1215 0.053 0.0222
1.367 0.172 0.235 7.066 1.054 0.0331 0.045 0.0135
1.272 0.207 0.205 7.152 1.098 0.0404 0.052 0.009
1.323 0.157 0.234 7.132 1.022 0.0447 0.058 0.0222
1.206 0.218 0.242 7.193 1.002 0.0404 0.054 0.0122
1.383 0.157 0.217 7.018 1.09 0.0546 0.06 0.0404
1.521 0.19 0.238 6.956 0.986 0.0366 0.05 0.01
1.35 0.168 0.227 6.986 1.165 0.0182 0.058 0.0055
1.422 0.181 0.218 6.813 1.23 0.0201 0.059 0.0082
1.298 0.166 0.224 6.986 1.162 0.03 0.065 0.0222
1.447 0.236 0.245 6.607 1.336 0.0271 0.051 0.0082
1.347 0.197 0.211 6.795 1.3 0.0245 0.059 0.03
1.339 0.17 0.253 6.989 1.11 0.0182 0.063 0.01
1.388 0.183 0.216 6.867 1.208 0.0164 0.061 0.0082
1.527 0.26 0.232 6.488 1.37 0.0222 0.045 0.0074
1.495 0.237 0.236 6.571 1.318 0.0245 0.058 0.0135
1.487 0.246 0.251 6.504 1.39 0.0182 0.053 0.0067
1.399 0.18 0.232 6.855 1.19 0.0245 0.066 0.009
1.489 0.215 0.242 6.777 1.145 0.0201 0.06 0.009
1.339 0.166 0.226 6.928 1.198 0.0201 0.06 0.01
1.482 0.246 0.238 6.444 1.462 0.0149 0.05 0.0074
1.434 0.172 0.255 6.646 1.354 0.0149 0.059 0.0122
1.347 0.156 0.214 6.85 1.313 0.0122 0.048 0.0067
1.34 0.158 0.233 6.848 1.272 0.0245 0.063 0.0122
1.453 0.18 0.244 6.752 1.238 0.03 0.054 0.01
1.306 0.149 0.226 7.082 1.097 0.0271 0.061 0.011
1.349 0.161 0.217 6.997 1.138 0.0222 0.062 0.01
1.254 0.151 0.205 7.319 0.947 0.0164 0.054 0.01
1.168 0.144 0.22 7.23 1.109 0.0222 0.052 0.0164
1.346 0.167 0.224 6.959 1.111 0.0201 0.049 0.01
1.39 0.184 0.212 6.898 1.189 0.0182 0.044 0.0067
1.283 0.149 0.224 7.077 1.104 0.0201 0.057 0.0245
1.214 0.137 0.232 7.269 1.005 0.0245 0.055 0.01
1.491 0.227 0.205 6.941 0.988 0.0271 0.068 0.03
1.479 0.218 0.207 7.039 0.887 0.0366 0.065 0.0366
1.445 0.174 0.228 6.875 1.123 0.0182 0.069 0.0222
1.439 0.183 0.218 6.775 1.226 0.0245 0.066 0.0182
1.387 0.154 0.204 6.991 1.09 0.03 0.074 0.0245
1.426 0.169 0.192 7.025 1.043 0.0222 0.064 0.0149
1.451 0.2 0.208 6.98 1.006 0.0201 0.062 0.0222
1.493 0.204 0.188 6.913 1.044 0.0245 0.061 0.0331
1.419 0.192 0.207 6.996 1.014 0.0366 0.07 0.0366
1.342 0.177 0.199 7.172 0.952 0.03 0.065 0.0271
1.349 0.152 0.236 7.145 0.949 0.0331 0.075 0.0182
1.44 0.196 0.208 6.938 1.07 0.0245 0.061 0.0135
1.46 0.215 0.197 6.918 1.081 0.0164 0.055 0.01
1.249 0.133 0.205 7.417 0.827 0.0271 0.072 0.0271
1.348 0.159 0.238 7.017 1.081 0.0222 0.067 0.0122
1.341 0.155 0.244 6.958 1.144 0.0245 0.068 0.0135
1.398 0.149 0.204 7.182 0.907 0.0182 0.076 0.0201
1.454 0.2 0.199 6.91 1.09 0.0201 0.062 0.0122
1.334 0.153 0.219 6.928 1.214 0.0271 0.066 0.011
1.438 0.204 0.189 7.107 0.91 0.0271 0.063 0.0149
1.303 0.138 0.212 7.17 1.016 0.03 0.069 0.0122
1.323 0.147 0.21 7.108 1.07 0.0271 0.061 0.0074
1.417 0.169 0.207 6.875 1.184 0.03 0.057 0.0149
1.36 0.167 0.225 6.883 1.22 0.0222 0.055 0.0149
1.42 0.179 0.214 6.923 1.121 0.0271 0.056 0.0149
1.472 0.218 0.214 6.724 1.238 0.0182 0.053 0.01
1.368 0.174 0.205 7.042 1.066 0.0222 0.057 0.0135
1.367 0.173 0.228 6.948 1.141 0.0245 0.053 0.011
1.403 0.173 0.209 6.843 1.21 0.0271 0.063 0.0271
1.413 0.197 0.206 6.737 1.387 0.03 0.06 0.0222
1.201 0.138 0.207 7.011 1.269 0.0404 0.064 0.0331
1.359 0.18 0.207 6.895 1.203 0.0271 0.061 0.0201
1.518 0.198 0.225 6.681 1.243 0.0182 0.057 0.011
1.434 0.185 0.189 6.771 1.269 0.0201 0.062 0.0122
1.367 0.162 0.179 6.772 1.368 0.0271 0.064 0.0149
1.461 0.181 0.197 6.783 1.246 0.0135 0.057 0.01
1.368 0.161 0.198 7.03 1.095 0.0271 0.059 0.0222
1.419 0.159 0.215 6.862 1.193 0.0331 0.06 0.0222
1.514 0.162 0.298 6.725 1.119 0.09 0.093 0.0201
1.328 0.171 0.253 6.987 1.03 0.0447 0.083 0.0494
1.469 0.16 0.337 6.675 1.127 0.0815 0.094 0.0366
1.283 0.153 0.196 7.107 1.115 0.0404 0.06 0.0164
1.263 0.155 0.199 7.14 1.148 0.0222 0.042 0.006
1.369 0.158 0.215 7.16 0.958 0.0447 0.069 0.0245
1.353 0.172 0.175 6.965 1.212 0.0164 0.075 0.0067
1.187 0.139 0.185 7.427 0.952 0.0182 0.056 0.009
1.732 0.231 0.156 6.437 1.313 0.09 0.062 0.01
1.62 0.255 0.166 6.628 1.212 0.0182 0.062 0.0149
1.543 0.172 0.193 6.74 1.157 0.1813 0.087 0.03
1.498 0.17 0.195 6.804 1.206 0.0331 0.066 0.01
1.399 0.169 0.171 7.011 1.1 0.0366 0.072 0.005
1.293 0.156 0.191 7.101 1.111 0.0245 0.06 0.0222
1.42 0.175 0.152 7.004 1.149 0.0149 0.05 0.0074
1.721 0.238 0.255 6.3 1.35 0.0331 0.07 0.0164
1.742 0.221 0.156 6.415 1.315 0.0737 0.082 0.01
1.391 0.187 0.189 6.975 1.062 0.1813 0.07 0.09
1.517 0.206 0.249 6.68 1.205 0.0271 0.08 0.0149
1.269 0.157 0.193 7.14 1.148 0.0222 0.04 0.006
1.577 0.204 0.208 6.732 1.183 0.0074 0.052 0.0074
1.59 0.241 0.195 6.705 1.149 0.0149 0.068 0.0082
1.621 0.28 0.197 6.608 1.179 0.0164 0.058 0.0149
1.753 0.275 0.236 6.367 1.214 0.01 0.061 0.0149
1.679 0.26 0.177 6.568 1.191 0.0201 0.059 0.0271
1.419 0.203 0.176 6.973 1.083 0.0447 0.078 0.0149
1.693 0.236 0.174 6.499 1.204 0.164 0.102 0.0404
1.692 0.27 0.234 6.499 1.196 0.0222 0.059 0.0045
1.638 0.252 0.215 6.57 1.199 0.0494 0.053 0.0182
1.497 0.247 0.219 6.621 1.27 0.0366 0.073 0.0245
1.442 0.222 0.194 6.677 1.314 0.0366 0.072 0.0447
1.68 0.27 0.17 6.44 1.31 0.0222 0.062 0.0164
1.463 0.164 0.185 6.909 1.154 0.1343 0.058 0.0055
;


************************************************** **;
* Check each variable separately for normality. *;
************************************************** **;
Title2 "Assessment of Marginal Normality";
Proc Univariate Data=olivetst;
Var palmitic palmitoleic stearic oleic
linoleic eicosanoic linolenic eicosenoic;
Symbol1 V=Dot;

Histogram / Normal(Mu=Est Sigma=Est Fill)
CFill=Yellow Name="olivetst";
Inset Mean Std / Header='Normal Parameters'
Position=(95,95) RefPoint=TR;

QQPlot / Normal(Mu=Est Sigma=Est L=1 NoPrint)
PctlMinor PCTLSCALE Name="olivetst";
Inset Mean Std / Header='Normal Parameters'
Position=(95,5) RefPoint=BR;

ProbPlot / Normal(Mu=Est Sigma=Est L=1 NoPrint)
PctlMinor Name="olivetst";
Inset Mean Std / Header='Normal Parameters'
Position=(95,5) RefPoint=BR;
Run;


************************************************** **;
* Examine bivariate normality through scatter *;
* plots of pairs of variables. *;
************************************************** **;
Title2 "Scatter Plot Matrix";
Proc Gplot Data=olivetst;
Plot (palmitic palmitoleic stearic oleic linoleic eicosanoic
linolenic eicosenoic)*(palmitic palmitoleic stearic oleic linoleic
eicosanoic linolenic eicosenoic)=1 / Name="olivetst";
Symbol1 C=Black V=Dot H=1;
Run;
Quit;

Goptions Reset=Symbol;

************************************************** **;
* Compute squared distances from each observation *;
* to the sample mean (centroid) and see if their *;
* distribution is approximately chi-squared. *;
* Here, we'll use the Interactive Matrix Language *;
* to compute the squared distances. *;
************************************************** **;
Title2 "Assessing Multivariate Normality";
Title3 'Generalized Distance Information';
Proc IML;
Reset NoLog; /* send output to the listing file */

/* Read the data into the matrix X */
use olivetst var{palmitic palmitoleic stearic oleic linoleic
eicosanoic linolenic eicosenoic};
read all var{palmitic palmitoleic stearic oleic linoleic eicosanoic
linolenic eicosenoic} into x;
close olivetst;

/* Number of Observations and Variables */
n=nrow(x);
p=ncol(x);

/* Compute sample mean, covariance, and inverse */
one=J(n,1,1);
xbar=(X`*one)/n;
print "Sample Means", xbar;
xstar=X-one*xbar`;
s=xstar`*xstar/(n-1);
print "Sample Covariance Matrix", s;
sinv=inv(s);
print "Inverse of Sample Covariance Matrix", sinv;
free xstar one;

/* Now compute the squared distances */
d2=J(n,1,1);
do i=1 to n;
dstar=(x(|i,|)-xbar`)`;
d2(|i|)=dstar`*sinv*dstar;
end;

/* Order the squared distances from smallest to largest */
temp=d2;
r=rank(d2);
d2(|r|)=temp;

Print "Ordered Squared Distances", d2;

/* Save this information in a data set */
create Distance var{d2};
append var{d2};
close Distance;

* Reset log;
Quit;

************************************************** **;
* Construct histogram and q-q and p-p plots of the *;
* squared distances to compare with the chi-squared*;
* distribution with 4 degrees of freedom (4 X's). *;
* NOTE: The chi-squared distribution is a special *;
* case of the gamma distribution. As parameterized *;
* by SAS, take ALPHA=df/2, SIGMA=2, and THETA=0 to *;
* get the chi-squared distribution with df degrees *;
* of freedom. *;
************************************************** **;

Title3 "Goodness of Fit of Chi-squared Distribution To Squared
Distances";
Proc Univariate Data=Distance;
Var D2;
Label d2="Squared Distance";
Symbol1 V=Dot;

Histogram / Gamma(Alpha=2 Sigma=2 Theta=0 Fill)
CFill=Yellow Name="olivetst";
Inset Gamma(KSD KSDPval) / Header='Goodness of Fit Test'
Position=(95,95) RefPoint=TR;

QQPlot / Gamma(Alpha=2 Sigma=2 Theta=0 L=1 NoPrint)
PctlMinor PCTLSCALE Name="olivetst";

ProbPlot / Gamma(Alpha=2 Sigma=2 Theta=0 L=1 NoPrint)
PctlMinor Name="olivetst";
Run;







Is that right? I'm still get errors.

Here is my log, if needed:

NOTE: This session is executing on the XP_PRO platform.



NOTE: SAS initialization used:
real time 1.46 seconds
cpu time 0.45 seconds

1 ************************************************** **;
2 * olivetst.sas
3 ************************************************** **;
4
5 Options PS=55 LS=80 PageNo=1 NoDate MPrint
6 FORMCHAR='|----|+|---+=|-/\<>*';
7
8 ************************************************** **;
9 * The following section generates HTML pages for *;
10 * web presentation. Delete the ODS and GOptions *;
11 * statements if you do not want the web pages. *;
12 ************************************************** **;
13
14 ODS HTML body= "olivetst-body.html"
15 contents="olivetst-contents.html"
16 frame= "olivetst-frame.html"
17 page= "olivetst-page.html"
18 headtext="<title>Multivariate Normality</title>"
19 anchor="olivetst";
NOTE: Writing HTML Body file: olivetst-body.html
ERROR: Insufficient authorization to access
C:\WINDOWS\system32\olivetst-body.html.
WARNING: No body file. HTML output will not be created.
20
21 GOptions Device=gif Transparency NoBorder
22 VSize=6 in HSize=6 in
23 HText=1 FText=Swiss HTitle=1 FTitle=Swiss;
24
25 ************************************************** **;
26 * The Data Analysis Begins Here.
27 ************************************************** **;
28 Title1 "Classification Analysis of Olive Data";
29
30 Data olivetst;
31 Infile palmitic palmitoleic stearic oleic
-----------
23
ERROR 23-2: Invalid option name PALMITOLEIC.

31 ! Infile palmitic palmitoleic stearic oleic
-------
23
ERROR 23-2: Invalid option name STEARIC.

31 ! Infile palmitic palmitoleic stearic oleic
-----
23
ERROR 23-2: Invalid option name OLEIC.

32 linoleic eicosanoic linolenic eicosenoic;
--------
23
ERROR 23-2: Invalid option name LINOLEIC.

32 ! linoleic eicosanoic linolenic eicosenoic;
----------
23
ERROR 23-2: Invalid option name EICOSANOIC.

32 ! linoleic eicosanoic linolenic eicosenoic;
---------
23
ERROR 23-2: Invalid option name LINOLENIC.

32 ! linoleic eicosanoic linolenic eicosenoic;
----------
23
ERROR 23-2: Invalid option name EICOSENOIC.

33 cards;

NOTE: The SAS System stopped processing this step because of errors.
WARNING: The data set WORK.OLIVETST may be incomplete. When this step
was
stopped there were 0 observations and 0 variables.
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds

241 ;


242
243
244 ************************************************** **;
245 * Check each variable separately for normality. *;
246 ************************************************** **;
247 Title2 "Assessment of Marginal Normality";
248 Proc Univariate Data=olivetst;
249 Var palmitic palmitoleic stearic oleic
ERROR: Variable PALMITIC not found.
ERROR: Variable PALMITOLEIC not found.
ERROR: Variable STEARIC not found.
250 linoleic eicosanoic linolenic eicosenoic;
ERROR: Variable OLEIC not found.
ERROR: Variable LINOLEIC not found.
ERROR: Variable EICOSANOIC not found.
ERROR: Variable LINOLENIC not found.
ERROR: Variable EICOSENOIC not found.
251 Symbol1 V=Dot;
252
253 Histogram / Normal(Mu=Est Sigma=Est Fill)
254 CFill=Yellow Name="olivetst";
255 Inset Mean Std / Header='Normal Parameters'
256 Position=(95,95) RefPoint=TR;
257
258 QQPlot / Normal(Mu=Est Sigma=Est L=1 NoPrint)
259 PctlMinor PCTLSCALE Name="olivetst";
260 Inset Mean Std / Header='Normal Parameters'
261 Position=(95,5) RefPoint=BR;
262
263 ProbPlot / Normal(Mu=Est Sigma=Est L=1 NoPrint)
264 PctlMinor Name="olivetst";
265 Inset Mean Std / Header='Normal Parameters'
266 Position=(95,5) RefPoint=BR;
267 Run;

NOTE: The SAS System stopped processing this step because of errors.
NOTE: PROCEDURE UNIVARIATE used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds

268
269
270 ************************************************** **;
271 * Examine bivariate normality through scatter *;
272 * plots of pairs of variables. *;
273 ************************************************** **;
274 Title2 "Scatter Plot Matrix";


275 Proc Gplot Data=olivetst;
276 Plot (palmitic palmitoleic stearic oleic linoleic eicosanoic
linolenic
276! eicosenoic)*(palmitic palmitoleic stearic oleic linoleic
eicosanoic
276! linolenic eicosenoic)=1 / Name="olivetst";
ERROR: Variable PALMITIC not found.
ERROR: Variable PALMITOLEIC not found.
ERROR: Variable STEARIC not found.
ERROR: Variable OLEIC not found.
ERROR: Variable LINOLEIC not found.
ERROR: Variable EICOSANOIC not found.
ERROR: Variable LINOLENIC not found.
ERROR: Variable EICOSENOIC not found.
ERROR: Variable PALMITIC not found.
ERROR: Variable PALMITOLEIC not found.
ERROR: Variable STEARIC not found.
ERROR: Variable OLEIC not found.
ERROR: Variable LINOLEIC not found.
ERROR: Variable EICOSANOIC not found.
ERROR: Variable LINOLENIC not found.
ERROR: Variable EICOSENOIC not found.
NOTE: The previous statement has been deleted.
277 Symbol1 C=Black V=Dot H=1;
278 Run;

NOTE: No variables in data set WORK.OLIVETST.
279 Quit;

NOTE: PROCEDURE GPLOT used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds


280
281 Goptions Reset=Symbol;
282
283 ************************************************** **;
284 * Compute squared distances from each observation *;
285 * to the sample mean (centroid) and see if their *;
286 * distribution is approximately chi-squared. *;
287 * Here, we'll use the Interactive Matrix Language *;
288 * to compute the squared distances. *;
289 ************************************************** **;
290 Title2 "Assessing Multivariate Normality";
291 Title3 'Generalized Distance Information';
292 Proc IML;
NOTE: IML Ready
293 Reset NoLog;
293! /* send output to the listing file */
294
295 /* Read the data into the matrix X */
296 use olivetst var{palmitic palmitoleic stearic oleic linoleic
eicosanoic
296! linolenic eicosenoic};
WARNING: Data set WORK.OLIVETST is empty.
ERROR: Variable PALMITIC is not on file WORK.OLIVETST.

statement : USE at line 296 column 1
297 read all var{palmitic palmitoleic stearic oleic linoleic
eicosanoic
297! linolenic eicosenoic} into x;
ERROR: No data set is currently open for input.

statement : READ at line 297 column 1
298 close olivetst;
NOTE: Cannot close WORK.OLIVETST; it is not open.
299
300 /* Number of Observations and Variables */
301 n=nrow(x);
302 p=ncol(x);
303
304 /* Compute sample mean, covariance, and inverse */
305 one=J(n,1,1);
ERROR: (execution) Invalid operand to operation.

operation : J at line 305 column 6
operands : N, *LIT1003, *LIT1004

N 1 row 1 col (numeric)

0

*LIT1003 1 row 1 col (numeric)

1

*LIT1004 1 row 1 col (numeric)

1

statement : ASSIGN at line 305 column 1
306 xbar=(X`*one)/n;
ERROR: (execution) Matrix has not been set to a value.

operation : `* at line 306 column 9
operands : X, ONE

X 0 row 0 col (type ?, size 0)


ONE 0 row 0 col (type ?, size 0)


statement : ASSIGN at line 306 column 1
307 print "Sample Means", xbar;
ERROR: Matrix XBAR has not been set to a value.

statement : PRINT at line 307 column 1
308 xstar=X-one*xbar`;
ERROR: (execution) Matrix has not been set to a value.

operation : *` at line 308 column 12
operands : ONE, XBAR

ONE 0 row 0 col (type ?, size 0)


XBAR 0 row 0 col (type ?, size 0)


statement : ASSIGN at line 308 column 1
309 s=xstar`*xstar/(n-1);
ERROR: (execution) Matrix has not been set to a value.

operation : `* at line 309 column 9
operands : XSTAR, XSTAR

XSTAR 0 row 0 col (type ?, size 0)


XSTAR 0 row 0 col (type ?, size 0)


statement : ASSIGN at line 309 column 1
310 print "Sample Covariance Matrix", s;
ERROR: Matrix S has not been set to a value.

statement : PRINT at line 310 column 1
311 sinv=inv(s);
ERROR: (execution) Matrix has not been set to a value.

operation : INV at line 311 column 9
operands : S

S 0 row 0 col (type ?, size 0)


statement : ASSIGN at line 311 column 1
312 print "Inverse of Sample Covariance Matrix", sinv;
ERROR: Matrix SINV has not been set to a value.

statement : PRINT at line 312 column 1
313 free xstar one;
314
315 /* Now compute the squared distances */
316 d2=J(n,1,1);
ERROR: (execution) Invalid operand to operation.

operation : J at line 316 column 5
operands : N, *LIT1009, *LIT1010

N 1 row 1 col (numeric)

0

*LIT1009 1 row 1 col (numeric)

1

*LIT1010 1 row 1 col (numeric)

1

statement : ASSIGN at line 316 column 1
317 do i=1 to n;
318 dstar=(x(|i,|)-xbar`)`;
319 d2(|i|)=dstar`*sinv*dstar;
320 end;
321
322 /* Order the squared distances from smallest to largest */
323 temp=d2;
ERROR: (execution) Matrix has not been set to a value.

operation : MOVE at line 323 column 1
operands : D2

D2 0 row 0 col (type ?, size 0)


statement : ASSIGN at line 323 column 1
324 r=rank(d2);
ERROR: (execution) Matrix has not been set to a value.

operation : RANK at line 324 column 7
operands : D2

D2 0 row 0 col (type ?, size 0)


statement : ASSIGN at line 324 column 1
325 d2(|r|)=temp;
ERROR: (execution) Matrix has not been set to a value.

operation : (| at line 325 column 3
operands : D2, R, TEMP

D2 0 row 0 col (type ?, size 0)


R 0 row 0 col (type ?, size 0)


TEMP 0 row 0 col (type ?, size 0)


statement : ASSIGN at line 325 column 1
326
327 Print "Ordered Squared Distances", d2;
ERROR: Matrix D2 has not been set to a value.

statement : PRINT at line 327 column 1
328
329 /* Save this information in a data set */
330 create Distance var{d2};
331 append var{d2};
WARNING: All data set variables are unvalued. No APPEND done.

statement : APPEND at line 331 column 1
332 close Distance;
NOTE: The data set WORK.DISTANCE has 0 observations and 1 variables.
333
334 * Reset log;
335 Quit;
NOTE: Exiting IML.
NOTE: The SAS System stopped processing this step because of errors.
NOTE: PROCEDURE IML used (Total process time):
real time 0.07 seconds
cpu time 0.07 seconds

336
337 ************************************************** **;
338 * Construct histogram and q-q and p-p plots of the *;
339 * squared distances to compare with the chi-squared*;
340 * distribution with 4 degrees of freedom (4 X's). *;
341 * NOTE: The chi-squared distribution is a special *;
342 * case of the gamma distribution. As parameterized *;
343 * by SAS, take ALPHA=df/2, SIGMA=2, and THETA=0 to *;
344 * get the chi-squared distribution with df degrees *;
345 * of freedom. *;
346 ************************************************** **;
347
348 Title3 "Goodness of Fit of Chi-squared Distribution To Squared
Distances";


349 Proc Univariate Data=Distance;
350 Var D2;
351 Label d2="Squared Distance";
352 Symbol1 V=Dot;
353
354 Histogram / Gamma(Alpha=2 Sigma=2 Theta=0 Fill)
355 CFill=Yellow Name="olivetst";
356 Inset Gamma(KSD KSDPval) / Header='Goodness of Fit Test'
357 Position=(95,95) RefPoint=TR;
358
359 QQPlot / Gamma(Alpha=2 Sigma=2 Theta=0 L=1 NoPrint)
360 PctlMinor PCTLSCALE Name="olivetst";
361
362 ProbPlot / Gamma(Alpha=2 Sigma=2 Theta=0 L=1 NoPrint)
363 PctlMinor Name="olivetst";
364 Run;

NOTE: No observations in data set WORK.DISTANCE.
NOTE: PROCEDURE UNIVARIATE used (Total process time):
real time 0.04 seconds
cpu time 0.01 seconds

On Oct 22, 3:03 pm, "mrfurle...@yahoo.com" <mrfurle...@yahoo.com>
wrote:
> Here is my log, if needed:
>
> NOTE: This session is executing on the XP_PRO platform.
>
> NOTE: SAS initialization used:
> real time 1.46 seconds
> cpu time 0.45 seconds
>
> 1 ************************************************** **;
> 2 * olivetst.sas
> 3 ************************************************** **;
> 4
> 5 Options PS=55 LS=80 PageNo=1 NoDate MPrint
> 6 FORMCHAR='|----|+|---+=|-/\<>*';
> 7
> 8 ************************************************** **;
> 9 * The following section generates HTML pages for *;
> 10 * web presentation. Delete the ODS and GOptions *;
> 11 * statements if you do not want the web pages. *;
> 12 ************************************************** **;
> 13
> 14 ODS HTML body= "olivetst-body.html"
> 15 contents="olivetst-contents.html"
> 16 frame= "olivetst-frame.html"
> 17 page= "olivetst-page.html"
> 18 headtext="<title>Multivariate Normality</title>"
> 19 anchor="olivetst";
> NOTE: Writing HTML Body file: olivetst-body.html
> ERROR: Insufficient authorization to access
> C:\WINDOWS\system32\olivetst-body.html.
> WARNING: No body file. HTML output will not be created.
> 20
> 21 GOptions Device=gif Transparency NoBorder
> 22 VSize=6 in HSize=6 in
> 23 HText=1 FText=Swiss HTitle=1 FTitle=Swiss;
> 24
> 25 ************************************************** **;
> 26 * The Data Analysis Begins Here.
> 27 ************************************************** **;
> 28 Title1 "Classification Analysis of Olive Data";
> 29
> 30 Data olivetst;
> 31 Infile palmitic palmitoleic stearic oleic
> -----------
> 23
> ERROR 23-2: Invalid option name PALMITOLEIC.
>
> 31 ! Infile palmitic palmitoleic stearic oleic
> -------
> 23
> ERROR 23-2: Invalid option name STEARIC.
>


What version of SAS are you using?

It looks to me like you are using an ancient version of SAS, in which
case you will need to truncate all variable names to a maximum of
eight characters.

In the ODS HTML statement, you would want to add a PATH= and GPATH= to
point to folders where your output should go. Something like:

ODS HTML path='e:\output' gpath='e:\output' body="olivetst-
body.html" ...

You pick the proper folder location.

--
Paige Miller
paige\dot\miller \at\ kodak\dot\com