cancel
Showing results for 
Show  only  | Search instead for 
Did you mean: 

Discussions

Solve problems, and share tips and tricks with other JMP users.
%3CLINGO-SUB%20id%3D%22lingo-sub-738591%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3ECrow-AMSAA%20Confidence%20Intervals%20-%20Delta%20Method%3C%2FLINGO-SUB%3E%3CLINGO-BODY%20id%3D%22lingo-body-738591%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3E%3CP%3EHello%2C%3C%2FP%3E%3CP%3E%26nbsp%3B%3C%2FP%3E%3CP%3EI%20am%20using%20JMP%2016.2.0%20to%20calculate%20Reliability%20Growth%20and%20am%20trying%20to%20understand%20the%20math%20behind%20the%20profiler%20confidence%20intervals%20on%20the%20Crow%20AMSAA%20models.%20The%20help%20article%20states%20that%3A%20%22%3CSPAN%3EThe%20delta%20method%20is%20applied%20to%20this%20expression%20to%20obtain%20an%20estimate%20of%20its%20variance.%20This%20estimate%20is%20used%20to%20construct%20a%2095%25%20Wald-based%20confidence%20interval.%20The%20resulting%20confidence%20limits%20are%20then%20transformed%20using%20the%20exponential%20function%20to%20give%20confidence%20limits%20for%20the%20estimated%20cumulative%20number%20of%20events%20at%20time%26nbsp%3B%3C%2FSPAN%3E%3CSPAN%20class%3D%22%22%3Et%3C%2FSPAN%3E%3CSPAN%3E.%22%20I%20have%20not%20been%20able%20to%20find%20an%20example%20of%20utilizing%20the%20delta%20method%20that%20I%20can%20follow.%20Is%20there%20a%20way%20to%20on%20JMP%20to%20look%20at%20the%20formula%2Fmath%20that%20goes%20into%20these%20calculations%3F%3C%2FSPAN%3E%3C%2FP%3E%3CP%3E%26nbsp%3B%3C%2FP%3E%3CP%3E%3CSPAN%3EThank%20you!%3C%2FSPAN%3E%3C%2FP%3E%3CP%3E%26nbsp%3B%3C%2FP%3E%3CP%3E%3CSPAN%3E-John%26nbsp%3B%3C%2FSPAN%3E%3C%2FP%3E%3C%2FLINGO-BODY%3E%3CLINGO-LABS%20id%3D%22lingo-labs-738591%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3E%3CLINGO-LABEL%3EReliability%20Analysis%3C%2FLINGO-LABEL%3E%3C%2FLINGO-LABS%3E%3CLINGO-SUB%20id%3D%22lingo-sub-743200%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3ERe%3A%20Crow-AMSAA%20Confidence%20Intervals%20-%20Delta%20Method%3C%2FLINGO-SUB%3E%3CLINGO-BODY%20id%3D%22lingo-body-743200%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3E%3CP%3EThank%20you%20for%20this%20information%20but%20I%20admit%20I%20am%20still%20lost.%20I%20am%20uncertain%20what%20the%20resulting%20formula%20of%20the%20delta%20method%20applied%20to%20the%20h(.)%20should%20give%20me.%20I%20am%20attempting%20to%20recreate%20confidence%20intervals%20for%20numerous%20time%20inputs%20simultaneously%2C%20so%20I%20would%20love%20to%20be%20able%20to%20input%20my%20parameters%20and%20times%20of%20interest%20in%20a%20table%20and%20have%20a%20JMP%20formula%20to%20create%20these%20confidence%20bands.%26nbsp%3B%3C%2FP%3E%3C%2FLINGO-BODY%3E%3CLINGO-SUB%20id%3D%22lingo-sub-742636%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3ERe%3A%20Crow-AMSAA%20Confidence%20Intervals%20-%20Delta%20Method%3C%2FLINGO-SUB%3E%3CLINGO-BODY%20id%3D%22lingo-body-742636%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3E%3CP%3EThe%20formula%2Fmath%20of%20the%20delta%20method%20can%20be%20found%20here%3A%20%3CA%20href%3D%22https%3A%2F%2Fen.wikipedia.org%2Fwiki%2FDelta_method%22%20target%3D%22_blank%22%20rel%3D%22nofollow%20noopener%20noreferrer%22%3Ehttps%3A%2F%2Fen.wikipedia.org%2Fwiki%2FDelta_method%3C%2FA%3E%3C%2FP%3E%0A%3CP%3EThe%20method%20can%20be%20applied%20well%20if%20the%20result%20of%20g(.)%20function%20is%20approximately%20normal.%3C%2FP%3E%0A%3CP%3EFor%20any%20g(.)%20function%20in%20hands%2C%20the%20above%20may%20not%20be%20true.%20Sometimes%20a%20transformation%20on%20g(.)%20function%20is%20necessary%2C%20e.g.%20h(.).%3C%2FP%3E%0A%3CP%3EFor%20such%20a%20situation%2C%20apply%20the%20delta%20method%20on%20h(g(.))%20first%2C%20then%20use%20the%20inverse%20of%20h(.)%20to%20get%20desired%20results%20of%20g(.).%3C%2FP%3E%0A%3CP%3EIn%20other%20words%2C%20get%20intervals%20of%20h(g(.))%2C%20say%20h_Low%20and%20h_High%2C%20using%20delta%20method%20on%20h(g(.))%2C%20then%20use%20h%5E(-1)(h_Low)%20and%20h%5E(-1)(h_High)%20to%20get%20the%20intervals%20for%20g(.)%2C%20where%20h%5E(-1)%20is%20the%20inverse%20of%20h(.).%3C%2FP%3E%0A%3CP%3EIn%20above%20JMP%20documentation%2C%20log%20function%20is%20the%20h(.).%3C%2FP%3E%3C%2FLINGO-BODY%3E%3CLINGO-SUB%20id%3D%22lingo-sub-743227%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3ERe%3A%20Crow-AMSAA%20Confidence%20Intervals%20-%20Delta%20Method%3C%2FLINGO-SUB%3E%3CLINGO-BODY%20id%3D%22lingo-body-743227%22%20slang%3D%22en-US%22%20mode%3D%22CREATE%22%3E%3CP%3EAt%20the%20moment%2C%20JMP%20does%20not%20produce%20a%20formula%20in%20JSL%20that%20a%20user%20can%20directly%20reproduces%20the%20results.%3C%2FP%3E%0A%3CP%3EThe%20mathematical%20details%20are%20available%20in%20statistical%20analysis%20books.%20My%20reference%20is%20%3CA%20href%3D%22https%3A%2F%2Fwww.amazon.com%2FStatistical-Reliability-Analysis-Probability-Statistics%2Fdp%2F1118115457%22%20target%3D%22_self%22%20rel%3D%22nofollow%20noopener%20noreferrer%22%3EStatistical%20Methods%20for%20Reliability%20Data%3C%2FA%3E%20.%20Section%20C.7.4%20and%20C.7.7%20in%20that%20book%20are%20what%20you%20need.%3C%2FP%3E%0A%3CP%3EThe%20needed%20results%2C%20however%2C%20are%20available%20in%20the%20report%2C%20especially%20the%20profilers%20in%20the%20report%3B%20see%20this%20screenshot%3A%3C%2FP%3E%0A%3CP%3E%3CSPAN%20class%3D%22lia-inline-image-display-wrapper%20lia-image-align-inline%22%20image-alt%3D%22peng_liu_0-1712061319408.png%22%20style%3D%22width%3A%20400px%3B%22%3E%3Cspan%20class%3D%22lia-inline-image-display-wrapper%22%20image-alt%3D%22peng_liu_0-1712061319408.png%22%20style%3D%22width%3A%20400px%3B%22%3E%3Cimg%20src%3D%22https%3A%2F%2Fcommunity.jmp.com%2Ft5%2Fimage%2Fserverpage%2Fimage-id%2F62968iD08E41F2E71FD7B5%2Fimage-size%2Fmedium%3Fv%3Dv2%26amp%3Bpx%3D400%22%20role%3D%22button%22%20title%3D%22peng_liu_0-1712061319408.png%22%20alt%3D%22peng_liu_0-1712061319408.png%22%20%2F%3E%3C%2Fspan%3E%3C%2FSPAN%3E%3C%2FP%3E%0A%3CP%3EBased%20on%20your%20description%2C%20seems%20that%20you%20would%20like%20to%20get%20a%20table%20which%20lists%20desired%20confidence%20intervals.%20I%20suggest%20looking%20into%20the%20%22Remember%20Settings%22%20feature%20in%20the%20profiler%3A%20%3CA%20href%3D%22https%3A%2F%2Fwww.jmp.com%2Fsupport%2Fhelp%2Fen%2F17.1%2F%23page%2Fjmp%2Fprediction-profiler-options.shtml%22%20target%3D%22_blank%22%20rel%3D%22noopener%20noreferrer%22%3Ehttps%3A%2F%2Fwww.jmp.com%2Fsupport%2Fhelp%2Fen%2F17.1%2F%23page%2Fjmp%2Fprediction-profiler-options.shtml%3C%2FA%3E%3C%2FP%3E%0A%3CP%3EFollowing%20is%20a%20screenshot%20of%20remembered%20settings%2C%20which%20record%20confidence%20intervals%20at%20three%20different%20time%20points.%3C%2FP%3E%0A%3CP%3E%3CSPAN%20class%3D%22lia-inline-image-display-wrapper%20lia-image-align-inline%22%20image-alt%3D%22peng_liu_1-1712061454183.png%22%20style%3D%22width%3A%20400px%3B%22%3E%3Cspan%20class%3D%22lia-inline-image-display-wrapper%22%20image-alt%3D%22peng_liu_1-1712061454183.png%22%20style%3D%22width%3A%20400px%3B%22%3E%3Cimg%20src%3D%22https%3A%2F%2Fcommunity.jmp.com%2Ft5%2Fimage%2Fserverpage%2Fimage-id%2F62969i571DD013F86AD42A%2Fimage-size%2Fmedium%3Fv%3Dv2%26amp%3Bpx%3D400%22%20role%3D%22button%22%20title%3D%22peng_liu_1-1712061454183.png%22%20alt%3D%22peng_liu_1-1712061454183.png%22%20%2F%3E%3C%2Fspan%3E%3C%2FSPAN%3E%3C%2FP%3E%0A%3CP%3E%26nbsp%3B%3C%2FP%3E%3C%2FLINGO-BODY%3E
Choose Language Hide Translation Bar

Crow-AMSAA Confidence Intervals - Delta Method

Hello,

 

I am using JMP 16.2.0 to calculate Reliability Growth and am trying to understand the math behind the profiler confidence intervals on the Crow AMSAA models. The help article states that: "The delta method is applied to this expression to obtain an estimate of its variance. This estimate is used to construct a 95% Wald-based confidence interval. The resulting confidence limits are then transformed using the exponential function to give confidence limits for the estimated cumulative number of events at time t." I have not been able to find an example of utilizing the delta method that I can follow. Is there a way to on JMP to look at the formula/math that goes into these calculations?

 

Thank you!

 

-John 

1 ACCEPTED SOLUTION

Accepted Solutions
peng_liu
Level VII

Re: Crow-AMSAA Confidence Intervals - Delta Method

At the moment, JMP does not produce a formula in JSL that a user can directly reproduces the results.

The mathematical details are available in statistical analysis books. My reference is Statistical Methods for Reliability Data . Section C.7.4 and C.7.7 in that book are what you need.

The needed results, however, are available in the report, especially the profilers in the report; see this screenshot:

peng_liu_0-1712061319408.png

Based on your description, seems that you would like to get a table which lists desired confidence intervals. I suggest looking into the "Remember Settings" feature in the profiler: https://www.jmp.com/support/help/en/17.1/#page/jmp/prediction-profiler-options.shtml

Following is a screenshot of remembered settings, which record confidence intervals at three different time points.

peng_liu_1-1712061454183.png

 

View solution in original post

3 REPLIES 3
peng_liu
Level VII

Re: Crow-AMSAA Confidence Intervals - Delta Method

The formula/math of the delta method can be found here: https://en.wikipedia.org/wiki/Delta_method

The method can be applied well if the result of g(.) function is approximately normal.

For any g(.) function in hands, the above may not be true. Sometimes a transformation on g(.) function is necessary, e.g. h(.).

For such a situation, apply the delta method on h(g(.)) first, then use the inverse of h(.) to get desired results of g(.).

In other words, get intervals of h(g(.)), say h_Low and h_High, using delta method on h(g(.)), then use h^(-1)(h_Low) and h^(-1)(h_High) to get the intervals for g(.), where h^(-1) is the inverse of h(.).

In above JMP documentation, log function is the h(.).

Re: Crow-AMSAA Confidence Intervals - Delta Method

Thank you for this information but I admit I am still lost. I am uncertain what the resulting formula of the delta method applied to the h(.) should give me. I am attempting to recreate confidence intervals for numerous time inputs simultaneously, so I would love to be able to input my parameters and times of interest in a table and have a JMP formula to create these confidence bands. 

peng_liu
Level VII

Re: Crow-AMSAA Confidence Intervals - Delta Method

At the moment, JMP does not produce a formula in JSL that a user can directly reproduces the results.

The mathematical details are available in statistical analysis books. My reference is Statistical Methods for Reliability Data . Section C.7.4 and C.7.7 in that book are what you need.

The needed results, however, are available in the report, especially the profilers in the report; see this screenshot:

peng_liu_0-1712061319408.png

Based on your description, seems that you would like to get a table which lists desired confidence intervals. I suggest looking into the "Remember Settings" feature in the profiler: https://www.jmp.com/support/help/en/17.1/#page/jmp/prediction-profiler-options.shtml

Following is a screenshot of remembered settings, which record confidence intervals at three different time points.

peng_liu_1-1712061454183.png

 

Recommended Articles