topic Graphing the probit procedure in Graphics Programming

Graphing the probit procedure

Pepperoni — Mon, 01 Sep 2014 18:54:48 GMT

Greetings! As a toxicologist I run the probit procedure in order to estimate the LC50 (the chemical concentration which is lethal to 50% of my test microscopic organisms). I get the output in both concentration units and log-transformed. I am interested in generating a graph showing the lethality curve for the log-transformed data. I don't know what code I can add into the proc to get such a graph. Could anyone help? I am copying my code with datalines below. Thank you very much!

data MAB111;

input form $ dead total conc @@;

prop = dead/total;

datalines;

A 1 30 0

A 15 35 1000

A 21 32 3000

A 34 39 5000

A 34 34 7000

A 36 36 9000

;

ods rtf;

proc sort;

by form;

run;

proc probit log10 inversecl lackfit;

model dead/total = conc/d = normal itprint;

output out=p p=prob std=std xbeta=xbeta;

by form;

run;

ods rtf close;

Re: Graphing the probit procedure

Reeza — Mon, 01 Sep 2014 19:19:39 GMT

Have you looked at the plots=all option and the default plots?

proc probit log10 inversecl lackfit plots=all;

model dead/total = conc/d = normal itprint;

output out=p p=prob std=std xbeta=xbeta;

by form;

run;

Re: Graphing the probit procedure

Pepperoni — Thu, 04 Sep 2014 20:22:35 GMT

No plots came up when I tried that. Does that require Java? There was no Java on the machine I used. Thanks.

Re: Graphing the probit procedure

Reeza — Fri, 05 Sep 2014 16:58:48 GMT

No idea. One thing to try is an ODS destination instead of listing,

Also, what version of SAS are you using?

ods html;

your code goes here;

ods html close;

Re: Graphing the probit procedure

Pepperoni — Sat, 06 Sep 2014 16:54:54 GMT

I am using SAS 9.3

No plots came up when using the html option.

Thanks so much for trying! I'll chart the data generated by SAS in Excel.

Best regards.

Re: Graphing the probit procedure

PGStats — Sat, 06 Sep 2014 18:32:51 GMT

Normally, the code

proc probit log10 inversecl lackfit plots=predpplot;

model dead/total = conc/ d = normal ;

output out=p p=prob std=std xbeta=xbeta;

by form;

run;

should produce the requested graph (it does for me)

ODS graphics definitely require Java. Look at the diagnostics from:

proc javainfo; run;

Also look seriously at the fit diagnostics from the procedure.