Turn on suggestions

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

Showing results for

- JMP User Community
- :
- Discussions
- :
- Discussions
- :
- Re: Accessing Control Chart constants/unbiasing constant

Topic Options

- Start Article
- 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
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Apr 15, 2019 4:29 AM
(1294 views)

I'm doing a variance analysis and need to look at the unbiased variance estimates for a number of subgroups I'm looking at. I can manually calculate them using a d2 table I have, but I was hoping there was a way to put that into a formula. I know that it can be derived from the Chi-squared distribution, but I don't have the math for that.

Any help on either a direct way to extract d2, or the math behind manually building it would be helpful.

Thanks

Evan Morris

1 ACCEPTED SOLUTION

Accepted Solutions

Highlighted

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Haven't seen that d2 thing in a while!

```
// https://v8doc.sas.com/sashtml/qc/chapc/sect9.htm
NamesDefaultToHere(1);
calculateD2 =
Function({n}, {Default Local},
d2 = Round(Integrate(1 - Power(1 - NormalDistribution(x), n) - Power(NormalDistribution(x), n), x, ., .), 5);
);
for(n=2, n<=25, n++, Print(n, calculateD2(n)));
```

6 REPLIES 6

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Re: Accessing Control Chart constants/unbiasing constant

Are you looking for the variance of a subgroup?

https://en.wikipedia.org/wiki/Variance#Sample_variance

That will be the square of the coumns standard deviation the formulas is (Col Std Dev( col, by1, by2, ...))^2

Run the script below. You can also find this using Table> Summary.

```
Names Default To Here( 1 );
dt = Open( "$SAMPLE_DATA/Big Class.jmp" );
dt << New Column( "Standard Variance for Each Age and Sex Group",
Formula( (Col Std Dev( :height, :age, :sex ))^2 )
);
```

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Re: Accessing Control Chart constants/unbiasing constant

No, I'm looking for the unbiased estimate of the population standard deviation (sigma-hat), that can be calculated from the subgroup sample standard deviation (s) and an unbiasing constant (either d2 or c4 depending on sample sizes and how I'm using them).

So lets say I have 30 subgroups of the same size (n=4). I can calculate S-bar by averaging the within subgroup standard deviation. I can then divide the s-bar by c4 and determine the unbiased estimate of the population standard deviation.

My problem is that I have varying subgroup sizes, so I can't just look it up once and put it in the table. I'm also a bit squirelly on whether or not I need d2 or c4 (but that's my problem, I'll sort it out).

Either way, I was just hoping those parameters were easily extractable by a function in JSL, they are pretty common in control charting/SQC work.

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Re: Accessing Control Chart constants/unbiasing constant

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Re: Accessing Control Chart constants/unbiasing constant

Have you taken a look at the JMP documentation Statistical Details section for the specific control charts you are using? Generally, JMP documentation provides background detail for specific questions such as yours in the Statistical Details section. Here's the section for X bar and S charts in the JMP online documentation:

https://www.jmp.com/support/help/14-2/control-limits-for-xbar-and-s-charts.shtml#1149056

Highlighted

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Haven't seen that d2 thing in a while!

```
// https://v8doc.sas.com/sashtml/qc/chapc/sect9.htm
NamesDefaultToHere(1);
calculateD2 =
Function({n}, {Default Local},
d2 = Round(Integrate(1 - Power(1 - NormalDistribution(x), n) - Power(NormalDistribution(x), n), x, ., .), 5);
);
for(n=2, n<=25, n++, Print(n, calculateD2(n)));
```

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Get Direct Link
- Email to a Friend
- Report Inappropriate Content

Re: Accessing Control Chart constants/unbiasing constant

This is exactly what I need, thanks.