BookmarkSubscribeRSS Feed
ST1
Calcite | Level 5 ST1
Calcite | Level 5

Hi Everyone,

I am trying to get the 95% confidence interval for a disease severity progression data set in SAS 9.4. The data have been log transformed and are not normal (certain weeks have a lot of 0 values).  Sometimes the confidence limit output contains "." rather than numbers. In the year 2 output it seems that this only happens when the disease severy mean in 0 (see year 2 below), but in the year 1 data it happens for weeks when the disease severity mean is 0.3010 or 0 (see year 1 below).  My question is: what does the "." mean? Why is a "." returned even when the mean is not 0?

Thank you

Sherri 🙂

 

This is the coding for year 2:

Title2 "Confidence Interval";
proc means data= first maxdec=4 alpha=0.05 clm mean;
class trt;
var jul23 jul31 aug6 aug13 aug20 aug27 sep3 sep10 sep17;
run;

 

and this is part of the output: 

year2.png

 

Coding for year 1

Title2 "Confidence Interval";
proc means data= first maxdec=4 alpha=0.05 clm mean;
class trt;
var aug21 aug30 sep14 sep24 oct11;
run;

 

Output for year 1:

year1.png

8 REPLIES 8
Astounding
PROC Star

Is it possible you have missing values in your data?  If only one observation has a nonmissing value, you can get a mean.  But you can't get confidence limits.  Perhaps you could check your data for a few of these cases.

ST1
Calcite | Level 5 ST1
Calcite | Level 5

There are 3 missing values in the year 1 data set (Oct 11 trt 2, and Sept 24 and Oct 11 trt 3), but they don't seem to correspond to the missing values in the confidence limit output. If it helps provide context, both data sets have 7 trt and four reps

Astounding
PROC Star

Perhaps you could share some of the data.  There are other conditions where there would be a mean but no confidence limit, such as all values being identical.  Are you able to share all the data points for one set of values that has no confidence limits?

ST1
Calcite | Level 5 ST1
Calcite | Level 5

Here is data for trt 7 of year 1 

block trt  aug21 aug30 sep14 sep24 oct11;
1	7	0	0	0.2	0.5	1
2	7	0	0	0	0.6	1
3	7	0	0.1	0	0.9	1
4	7	0	0	0	0.5	1
;
Reeza
Super User

Add Standard Deviation to the chart please and repost. 

ST1
Calcite | Level 5 ST1
Calcite | Level 5

Please see orginal message for update chart.

Reeza
Super User

Do you see the issue? 🙂

 

All results with missing CI have a standard deviation of 0 - implying that the values are all that single value.

Ksharp
Super User
As Reeza pointed out . You can not calculated CI when it is at the highest level (nway), there is only one obs for that. Specify WAYS statement.
proc means data= first maxdec=4 alpha=0.05 clm mean;
class trt;
ways 1 ;
var jul23 jul31 aug6 aug13 aug20 aug27 sep3 sep10 sep17;
run;

SAS Innovate 2025: Register Now

Registration is now open for SAS Innovate 2025 , our biggest and most exciting global event of the year! Join us in Orlando, FL, May 6-9.
Sign up by Dec. 31 to get the 2024 rate of just $495.
Register now!

What is Bayesian Analysis?

Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.

Find more tutorials on the SAS Users YouTube channel.

SAS Training: Just a Click Away

 Ready to level-up your skills? Choose your own adventure.

Browse our catalog!

Discussion stats
  • 8 replies
  • 2358 views
  • 0 likes
  • 4 in conversation