Could anyone explain what's the difference when calculate the mean and CI of "class" variable in proc means and proc glm. For example I want to calculate CI for 2 levels Y and N, (Note: not the difference between Y and N) I have 2 codes give me 2 different CI Code 1 : Proc means data=KK clm alpha=0.05 ; Var UOSMOLVN; class lotresp; Output out=textn n=n mean=mean lclm=lclm uclm=uclm std=std ; run;                                                       Output 1: Analysis Variable : UOSMOLVN Osmolality(mOsm/kg) Numerical ValueLot AEOccurrence(Yes/No) N Obs Lower 95%CL for Mean Upper 95%CL for MeanN 116Y 18 344.1497139 350.8675275 328.3620766 344.0823679 Code 2 : proc glm data=KK; class lotresp; model UOSMOLVN=lotresp; means lotresp/clm t ; quit;                                                        output 2 :  lotresp N Mean 95% Confidence LimitsNY 116 347.509 344.209 350.808 18 336.222 327.846 344.599   Both codes generate the same means, but why CI are different?     Thanks in advance. H     ... View more

 I use a simple example to illustrate my question.   I want to check the fdr for a raw_palue, say 0.61, using the code below:  data a; input Test$ Raw_P @@; datalines; test01  0.61 ; proc multtest inpvalues=a bon fdr ; run; p-ValuesTest Raw Bonferroni False DiscoveryRate1 0.6100 0.6100 0.6100 But if I use another dataset:    data a; input Test$ Raw_P @@; datalines; test01 0.001 test02 0.001 test03 0.001 test04 0.02 test05 0.22 test06 0.59 test07 0.87 test08 0.61 ; proc multtest inpvalues=a bon fdr ; run;   the output is  p-ValuesTest Raw Bonferroni False DiscoveryRate12345678 0.0010 0.0080 0.0027 0.0010 0.0080 0.0027 0.0010 0.0080 0.0027 0.0200 0.1600 0.0400 0.2200 1.0000 0.3520 0.5900 1.0000 0.6971 0.8700 1.0000 0.8700 0.6100 1.0000 0.6971  So the fdr for 0.61 become 0.6971?     Can anyone tell me why different dataset, with the same value, give me the different fdr?  Did I do something wrong?    Best,    H ... View more

Hi,  I try to compare convergence rate using different model procedure: logistic, genmod, and glimmix.  However each procedure has difference convergence criterion.  I listed some example below.   Are these criterion similar to each other? Can anyone give me advice which one I should use ?  Thanks so much in advance !    Proc genmod: CICON, CONVH, CONVER, default is 1e-4, INITI Proc logistic: GCONV : default is 1e-8 proc glimmix:  pconv, default is 1e-4   ... View more

I try to use proc logistic regression to estimate the probability difference between placebo and treatment group. It gives me odds ratio difference, how can I change to risk difference?    Also any difference between "risk difference" or "proportional difference"  or "probability difference" between proc freq and proc logistic?  Thanks,  ... View more

Hi  I need to do simulation analysis for survival study. I have one real clinical dataset with only 40 patients. need to create 1000 simulation samples.    Thanks ... View more