turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Variance Analysis

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

11-10-2015 09:06 AM - edited 11-11-2015 08:09 AM

Hello all,

I am editing my question since I had progress with it.

I have lab results, the measurement is Y, and I have some factors X1, X2, X3,...

In addition, for each combination, a duplicate of measurements.

I am interested in the variance components, meaning, I want the % of variation that can be addressed to each factor.

My code is:

```
My code is:
proc varcomp method=reml data=RR;
class X1 X2 X3 X4;
model Y=X1 X2 X3 X4 / fixed=0;
run;
or
proc mixed data = RR;
class X1 X2 X3 X4;
model Y = ;
random X1 X2 X3 X4;
run;
```

I tried running this in JMP (just to validate, I must use SAS on this one), and JMP said that the REML estimators are negative and therefore it used Bayesian estimators instead. By doing so, the estimate of the total variability was not far from the observed variaility across all samples. When I forced JMP to stick with REML, it was identical to SAS, but I think wrong. In addition, I only got the right variability when I asked for nested factors.

My questions are:

1. How do I ask SAS for Bayesian estimators ?

2. How do I specify nested factors ?

3. If in addition to the variability between each factor (between lots, between days,...) I want also the CV, how do I obtain the mean for diviging the SD with ?

thank you