ProProcessCapability.FitJohnsonByPercentile

SVB Example

Estimate the parameters of the Johnson distribution by percentiles:

Sub Main

    Dim newanalysis As Analysis
    Set newanalysis = Analysis (scProcessAnalysis, _
        Application.Spreadsheets.OpenReadOnly(Path & "\Examples\Datasets\Pistons2.sta") )
    Dim oStaDocs As StaDocuments

    ' Process Analysis Procedures: Pistons2
    Dim oAD1 As STAProcessAnalysis.ProStartup
    Set oAD1 = newanalysis.Dialog
    oAD1.TypeOfAnalysis = scProRawDataProcessCapabilityAndToleranceIntervals

    newanalysis.Run

    ' Process Capability Analysis Setup--Raw Data: Pistons2
    Dim oAD2 As STAProcessAnalysis.ProProcessCapability
    Set oAD2 = newanalysis.Dialog
    oAD2.Variables = "1 "
    oAD2.ProcessSpecifications = """SAMPLE_1"" ( 74 0.05 MISSING MISSING ) "
    oAD2.ConstantSampleSize = True
    oAD2.ValueOfConstantSampleSize = 1
    oAD2.EstimateSigmaFromRanges = True
    oAD2.BetaOffset = 0
    oAD2.BetaScale = 1
    oAD2.ExponentialOffset = 0
    oAD2.GammaOffset = 0
    oAD2.LogNormalOffset = 0
    oAD2.NormalAndGeneralNonNormal = True
    'Estimate the parameters of the Johnson distribution via the method outlined in Chou, Polansky, and Mason (1998)
    oAD2.FitJohnsonByPercentile = True
    oAD2.RayleighOffset = 0
    oAD2.WeibullOffset = 0
    oAD2.FoldedNormalOffset = 0
    oAD2.MethodForModelA1 = scProCapabilityTimeDistributionModelISOMethodM141
    oAD2.MethodForModelA2 = scProCapabilityTimeDistributionModelISOMethodM1l5
    oAD2.MethodForModelB = scProCapabilityTimeDistributionModelISOMethodM1l5
    oAD2.MethodForModelC1 = scProCapabilityTimeDistributionModelISOMethodM144
    oAD2.MethodForModelC2 = scProCapabilityTimeDistributionModelISOMethodM1l5
    oAD2.MethodForModelC3 = scProCapabilityTimeDistributionModelISOMethodM1l5
    oAD2.MethodForModelC4 = scProCapabilityTimeDistributionModelISOMethodM1l5
    oAD2.MethodForModelD = scProCapabilityTimeDistributionModelISOMethodM1l5

    newanalysis.Run

    ' Process Capability Analysis--Normal and General Non-Normal Distribution: Pistons2
    Dim oAD3 As STAProcessAnalysis.ProProcessCapabilityResults
    Set oAD3 = newanalysis.Dialog
    oAD3.TypeOfProcessSpecifications = scProNominalDelta
    oAD3.NominalIsDefined = True
    oAD3.NominalValue = 74
    oAD3.DeltaValue = 0.05
    oAD3.UserDefinedMean = False
    oAD3.UserDefinedSigma = False
    oAD3.SigmaLimits = 6
    oAD3.AverageTiedRanks = False
    oAD3.QuantilesProbabilitiesOnHorizontalAxis = True
    oAD3.ValueForAdjustingRanks = -0.33333333
    oAD3.ValueForAdjustingN = 0.33333333
    oAD3.NumberOfCasesForToleranceLimits = 100
    oAD3.PercentIncludedForToleranceLimits = 95
    oAD3.ConfidenceLevelForToleranceLimits = 95
    oAD3.BaseOnPercentOfCasesIncluded = True
    oAD3.FrequencyTabulationLowerLimit = 73.964999978431
    oAD3.FrequencyTabulationUpperLimit = 74.024999977089
    oAD3.FrequencyTabulationNumberOfCategories = 13
    oAD3.ShowProcessCapabilityIndices = True
    oAD3.SigmaLimitsAtEquivalentPercentiles = True
    oAD3.NeatIntervals = True
    oAD3.MoveTitlesToAvoidOverlappingText = True

    Set oStaDocs = oAD3.Summary
    newanalysis.RouteOutput(oStaDocs).Visible = True
    Set oStaDocs = Nothing

End Sub