Hi, All I want to create a cross table to summaries favorite before and after test the tasty and put stage 1 and stage 2 together. I have created codes to create the result, but could not report it as follow in rtf file, I am wondering if anyone can help?   rtf ouput I want as follow:         The code as following: data raw; input (Subject Stage Before After) ($); cards; 1 I Like Like 2 II Abit Like 3 I Like Abit 4 II Abit Like 5 I Like Abit 6 II Like Like 7 II Abit Like 8 I Like Dislike 9 II Like Abit 10 I Dislike Like 11 II Abit Dislike 12 I Like ND ; run; proc sort data=raw; by Stage; run; data class; format Before $8. After $8.; do Before='Like', 'Abit', 'Dislike', 'ND'; do After='Like', 'Abit', 'Dislike', 'ND'; output; end; end; run; proc format; value $FBefore 'Like'='Like' 'Abit'='A bit like' 'Dislike'='Dislike' 'ND'='Not Done'; value $FAfter 'Like'='Like' 'Abit'='A bit like' 'Dislike'='Dislike' 'ND'='Not Done'; run; proc tabulate data=raw missing classdata=class; format Before $FBefore. After $FAfter.; class Before After/order=data; table After, Before*(N=' ')/misstext='0'; by Stage; run;   ... View more

Dear all   I am trying to create a cross table to see how food test change, I used following code:   data raw;   input person firsttest $ secondtest $;   datalines;   1 Good Good   2 Poor Good   3 Good Bad   4 Poor Bad   5 Poor Poor ; run;   proc freq data=raw;   table firsttest*secondtest/nocol norow nopercent; run;   due to results in firttest dose not have "Bad" result, therefore, it shows:   Table of firsttest by secondtest 　 secondtest firsttest Bad Good Poor Total Good 1 1 0 2 Poor 1 1 1 3 Total 2 2 1 5   I am wondering how to display a full result as following table:   Table of firsttest by secondtest 　 secondtest firsttest Bad Good Poor Total  Bad 0 0 0 0 Good 1 1 0 2 Poor 1 1 1 3 Total 2 2 1 5 I also tried to use prco tabulate to do it, but it failed.   Appreciate the help.     Charlie ... View more

I am trying to do probit analysis for 95% confidence Limit and it is my first time to use this procedure.   I've got an example in .xls file which is downloaded from https://probitanalysis.wordpress.com/2016/07/07/first-blog-post/.   I used codes as following (reference to LD50 calculation for mixture of dose concentrations in community😞 data A; input dose total dead; cards; 0.03 10 0 0.10 10 1 0.30 10 5 3.00 10 8 ; proc probit data=A log10 plot=predpplot; model Dead / Total = dose / d=logistic InverseCL(prob=0.1 to 0.9 by 0.1); run; However, the results of 95% Fiducial Limits are different from example and proc probit dose not provide information for curve fitting ( I need to use proc reg?). How can I reproduce the result from the example file?   Furthermore, I used the same code for my data; data B; input dose total dead; cards; 0.64 12 9 0.8 12 10 1 12 8 1.25 12 9 1.56 12 9 1.95 12 8 ; proc probit data=B log10 plot=predpplot; model dead / total = dose / d=logistic InverseCL(prob=0.1 to 0.9 by 0.1); run; The 95% Fiducial Limits didnot estimate in the result and there was no warring from SAS. I am wondering if I could find any suggestion or solution for my data.           ... View more

I use SAS program (Version 9.4) and Phoenix/Winnonlin to calculate  two one-sided T test. The results of point estimate and 90% of confidence interval of SAS and Phoenix/Winnonlin are same. But "TOST" have different value of upper and lower. Here are the results of TOST value of SAS and Phoenix/Winnonlin. SAS:   Result in TOST Level 0.1 Equivalence Analysis shows: SAS code as follow ods graphics off; proc ttest data=temp1_1 dist=lognormal tost(1.0114488825,1.2793655184) alpha=0.1; class Treatment; var lnCmax; run;   Phoenix/Winnonlin There are questions: 1. Why the results are different? 2. How to calculate t1_TOST, t2_TOST, Prob_80_00 and Prob_125_00 in Phoenix by using SAS?   The data as attachment. Please help me answer these questions. Thank you. ... View more