# Approximate Calculation Method for Noncentral t-Distribution Quantile

Authors: Jun Gao; Jitao Yao
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

In the process of structural design and structural performance evaluation, the inference of the reliable life of the structure and the representative value of the material strength is a necessary work. The determination of material strength is the presumption of the quantile of the normal distribution, and the determination of the confidence level of the quantile of the normal distribution involves the noncentral t-distribution function. However, the calculation of the quantile is very complicated and is often provided in the form of a numerical table, which often involves multiparameter interpolation calculation, so it is not convenient to apply. The existing approximate calculation methods for noncentral t-distribution quantiles have strict application conditions, and the calculation process is relatively cumbersome. It is still difficult to meet actual needs in terms of fitting accuracy, application range, and convenience. In this paper, a new calculation method for noncentral t-distribution quantiles is proposed by introducing new probability expressions and related approximate distributions, based on theoretical derivation and numerical fitting. The comparative analysis results show that the method not only is convenient for calculation but also has the advantages of higher accuracy and wider application range, and it is more in line with the actual needs of engineering.


---

## Body

## 1. Introduction

The structural reliability design covers the entire civil engineering and is the basic guarantee for the normal operation of the building structures [1, 2]. There are various reliability modeling methods. In practical engineering, uncertain failure modes can result in the system failure [3–5]. The noncentral t-distribution is mainly used for small sample inference of normal distribution quantile, and it is widely used in the estimation of structural reliability life, representative value of material strength, and so on [6, 7]. The probability density function of the noncentral t-distribution is complex, and its quantile is usually calculated by numerical methods and provided in the form of a numerical table. However, in the current numerical table, only some typical values are provided, and in other cases, interpolation calculation is required, and multiple parameters are involved, which is inconvenient to apply [8–11]. Therefore, a variety of approximate calculation methods have been proposed, but they are still difficult to meet the actual needs in terms of fitting accuracy, scope of application, and convenience [12, 13].In this paper, a new calculation method for noncentral t-distribution quantiles is proposed by introducing new probability expressions and related approximate distributions, based on theoretical derivation and numerical fitting. The comparative analysis results show that the method not only is convenient for calculation, but also has the advantages of higher accuracy and wider application range, and it is more in line with the actual needs of engineering.


## 2. Current Approximate Calculation Method

The random variableX obeys the normal distribution N(μ, σ2), and the distribution parameters μ and σ2 are unknown. Then, the upper p and 1 − p quantile values of X can be expressed as(1)xp=μ−z1−pσ,x1−p=μ+z1−pσ,where z1−p is the upper 1 − p quantile of the standard normal distribution; p is the assurance rate of X ≥ xp or X ≤ x1−p, and its value is generally more than 0.5, and xp and x1−p are small and large quantile values, respectively.LetX1, X2, …, Xn be the n samples of the random variable X. According to the statistical method, after obtaining the sample realization value x1, x2, …, xn, the lower limit estimate of xp and the upper limit estimate of x1−p are, respectively,(2)x^p=x¯−ks,x^1−p=x¯+ks,where(3)x¯=1n∑i=1nxi,s=1n−1∑i=1nxi−x¯2,k=tn−1,z1−pn,αn,where tn−1,z1−pn,α is the upper α quantile value of the noncentral t-distribution with a degree of freedom of X and a parameter of z1−pn, and α is the significance level.The coefficientk is related to the sample size n, the distribution parameter z1−pn and the significance level α, etc., and is mainly determined by the quantile value tn−1,z1−pn,α of the noncentral t-distribution. According to the approximate statistics of the F distribution [12], it is recommended to approximate as(4)tn−1,λ,α=λCnFγ1,γ1,α,λ>0,λCnFγ1,γ1,1−α,λ<0,where(5)Cn=2n−1Γn/2Γn−1/2,γ1=2λ2,γ2=2Cn21−Cn2,where tn−1,γ,α is the upper α quantile value of the noncentral t-distribution with the degree of freedom n − 1 and the parameter γ. Fγ1,γ1,α and Fγ1,γ1,1−α are the upper α and 1 − α quantile values of the F distribution with parameters γ1 and γ2, respectively. When γ1 and γ2 are nonpositive integers, the values of Fγ1,γ1,α and Fγ1,γ1,1−α can be determined by linear interpolation method. This method is mainly applicable to the case of n > 9, and, in general, it is still inconvenient to check the table.When(6)T1=X¯−xpS/n,X¯=1n∑i=1nXi,S=1n−1∑i=1nXi−X¯2.Then,T1 obeys a noncentral t-distribution with a degree of freedom of n − 1 and a parameter of z1−pn. Jitao [14] approximates that the probability distribution of S is a normal distribution Nσ,σ2/2n−1. By numerical fitting, it is recommended to approximate the noncentral t-distribution of the quantile value as(7)tn−1,γ,α=γ+zα1−zα2−γ2/2n−m1−zα2/2n−mzα2≤2n−m,(8)m=1.51−α2,where zα is the upper α quantile value of the standard normal distribution.This method is mainly applicable to the case ofn ≥ 5, has a wider scope of application, and is simpler than the method proposed by Huirong [12]. But it can be known from the comparison analysis of the accuracy of the following, the approximation methods proposed by Huirong [12] and Jitao and Yaokui [13] are difficult to meet the actual needs in terms of accuracy, and the scope of application is relatively limited.


## 3. New Approximate Calculation Method

If(9)T2=U+γV/n−1,where U obeys the standard normal distribution and V obeys the chi-square distribution with degrees of freedom of n – 1, then, T2 obeys a noncentral t-distribution with a degree of freedom of n − 1 and a parameter of γ. When(10)PT2≤tn−1,γ,α=PU+γV/n−1≤tn−1,γ,α=1−α,W=Vn−1.Then,(11)PU+γ−tn−1,γ,αW≤0=1−αThe probability density function and mean value ofV are(12)fVv=12n−1/2Γn−1/2e−v/2vn−1/2−1,v≥0,μV=n−1.So, the probability density function ofW is(13)fWw=2n−1wfVn−1w2,w≥0.Therefore,(14)μW=∫0∞wfWwdw=∫0∞vn−1fVvdv=2n−1Γn/2Γn−1/2,EW2=∫0∞w2fWwdw=∫0∞vn−1fVvdv=1,σW=EW2−μW2=1−2n−1Γ2n/2Γ2n−1/2.Referring to the method proposed by Shisong et al. [15] and Jitao [14], approximating that W obeys the normal distribution NμW,σW2, then(15)PU+γ−tn−1,γ,αW≤0≈Φ−γ−tn−1,γ,αμW1+tn−1,γ,α2σW2=1−α,tn−1,γ,αμW−γ1+tn−1,γ,α2σW2≈zα.It can be obtained that(16)tn−1,γ,α=γ+zα1−zα2−γ2δW2μW1−zα2δW2,(17)δW=σWμW=n−12Γ2n−1/2Γ2n/2−1≈1−0.25n−1.After the overall numerical fit, it can be found in equation (16):(18)μW=1−0.25n−1p+0.6−0.1α,p=Φγn.If(19)γ=z1−pn,then(20)k=z1−pn+zα1−zα2−z1−p2nδW2μWn1−zα2δW2.


## 4. Comparative Analysis

Compare and analyze the calculation accuracy of the four methods: method 1: the method proposed by Huirong [12] calculates k according to (4); method 2: the method proposed by Jitao [13] calculates k according to (7); method 3: the method in this paper calculates k according to (20). Tables 1−5 list the corollary results of the coefficient k according to the different values of C.Table 1
Corollary results of coefficientk (C = 0.5).
p0.7nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)30.5979860.415034−30.590.524401−12.310.6080011.6740.5724870.435142−23.990.524401−8.400.5773220.8450.5600640.449987−19.650.524401−6.370.5631150.5460.5526640.460839−16.61——0.5549210.4170.5478270.469046−14.380.524401−4.280.549590.3280.5444070.475423−12.670.524401−3.680.5458440.2690.5418680.480519−11.320.524401−3.220.5430680.22100.5398770.484634−10.230.524401−2.870.5851110.19200.531660.503979−5.210.524401−1.370.5626830.08300.5291470.510691−3.490.524401−0.900.5542240.05400.5279070.514103−2.610.524401−0.660.5495340.04Error range−30.59∼−2.61−12.31∼−0.660.04∼1.67p0.8nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)30.9707150.88509−8.820.841621−13.30.990143240.9245470.859236−7.060.841621−9.00.9351351.1550.9026060.849782−5.85——0.9098610.860.8898440.845484−4.99——0.8953420.6270.8814790.843315−4.330.841621−4.520.8859170.580.8755640.842031−3.830.841621−3.880.8793060.4390.8712020.841351−3.430.841621−3.40.8744120.37100.8677980.840878−3.10.841621−3.020.8706430.33200.8538530.840228−1.60.841621−1.430.8551230.15300.8495850.840502−1.070.841621−0.940.8504190.1400.8475340.84068−0.810.841621−0.70.8481450.07Error range−8.82∼−0.81−13.3∼−0.70.07∼2p0.9nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)31.4985451.481939−1.111.281552−14.51.5302112.1141.4188831.403378−1.09——1.4372541.2951.3818841.368009−1——1.3948860.9461.3604191.348175−0.9——1.3706430.7571.3465661.335561−0.821.281552−4.831.3549440.6281.3367291.326928−0.731.281552−4.131.3439490.5491.329651.320672−0.681.281552−3.621.3358190.46101.3241261.315862−0.621.281552−3.221.3295640.41201.3012721.29683−0.341.281552−1.521.3038540.2301.2944491.291352−0.241.281552−11.2960770.13401.2910371.28876−0.181.281552−0.731.2923210.1Error range−1.11∼−0.18−14.5∼−0.730.1∼2.11p0.95nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)31.9384451.9549450.85——1.9787712.0841.8295251.838530.49——1.8533561.351.779261.7846940.31——1.7964270.9661.750431.7538740.20——1.7639180.7771.7318161.7342340.141.644854−5.021.7428910.6481.7187221.720310.091.644854−4.31.7281760.5591.7090371.7103370.081.644854−3.761.7173020.48101.7016621.7025890.051.644854−3.341.7089390.43201.6711631.6711301.644854−1.571.67460.21301.6619821.66191301.644854−1.031.6642220.13401.6575361.65741701.644854−0.771.6592120.1Error range0∼0.85−5.02∼−0.770.1∼2.08Table 2
Corollary results of coefficientk (C = 0.6).
p0.7nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)30.7998950.58474−26.90.627209−21.60.8018970.2540.7332650.578577−21.1——0.7336470.0550.6980540.577128−17.320.7028670.690.698110.0160.6756720.576407−14.690.659712−2.360.675599−0.0170.6597580.575744−12.730.642013−2.690.659767080.6478480.574975−11.250.630945−2.610.6478810.0190.638520.574221−10.070.622911−2.440.638550100.6309960.57343−9.120.616635−2.280.6309810200.5941810.566305−4.690.586898−1.230.5941980300.5796870.561364−3.160.574889−0.830.579710400.5714850.557856−2.380.567906−0.630.5714980Error range−26.9∼−2.38−21.6∼0.690∼0.25p0.8nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)31.2112671.111396−8.250.867739−28.361.2161480.441.1109311.038738−6.5——1.1137180.2551.0602451.003604−5.341.0822992.081.0624970.2161.029140.982315−4.551.010636−1.81.0309330.1771.0075690.967799−3.950.983455−2.391.0091840.1680.9916410.957063−3.490.967535−2.430.9931130.1590.979350.948756−3.120.95652−2.330.9806580.13100.9695020.942029−2.830.948208−2.20.9706620.12200.9228910.909472−1.450.911562−1.230.9235350.07300.9052380.89636−0.980.897651−0.840.9056740.05400.8954060.888792−0.740.889712−0.640.8957350.04Error range−8.25∼−0.74−28.36∼2.080.04∼0.4p0.9nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)31.8053991.77636−1.61——1.8122080.3841.64971.628998−1.25——1.6553570.3451.5742051.557654−1.051.6187572.831.5792050.3261.5289631.514778−0.931.506009−1.51.5332320.2871.4982671.485975−0.821.464439−2.261.5020450.2581.4759271.465039−0.741.440912−2.371.4792850.2391.4586421.449051−0.661.42511−2.31.4618250.22101.4451121.436289−0.611.413474−2.191.4479340.2201.3824941.377987−0.331.365121−1.261.3840650.11301.3595981.356631−0.221.347844−0.861.3606580.08401.3470751.344737−0.171.338184−0.661.3478530.06Error range−1.61∼−0.17−2.37∼2.830.06∼0.38p0.95nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)32.306762.3083310.07——2.3128660.2642.1021182.1049010.13——2.1087690.3252.0047262.0065780.092.0656853.042.0109190.3161.9465731.9480.071.918951−1.421.9523690.371.907811.90860.041.86521−2.231.9129190.2781.8796921.8801880.031.835126−2.371.8842890.2491.8581261.858480.021.815139−2.311.8624260.23101.8412391.84125801.800564−2.211.84510.21201.7642251.764128−0.011.741567−1.281.7663580.12301.7364881.73640601.721125−0.881.7379630.08401.7214371.721317−0.011.709825−0.671.722560.07Error range−0.01∼0.07−2.37∼3.040.07∼0.26Table 3
Corollary results of coefficientk (C = 0.7).
p0.7nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)31.0568560.823382−22.09——1.041203−1.4840.928040.770563−16.970.9829415.920.919138−0.9650.8610970.742317−13.790.850421−1.240.855142−0.6960.8186490.723444−11.630.799625−2.320.814274−0.5370.7886320.709445−10.040.769229−2.460.785322−0.4280.7661560.698382−8.850.747888−2.380.763447−0.3590.7484130.689293−7.90.731641−2.240.74618−0.3100.734010.681582−7.140.718652−2.090.732107−0.26200.6634380.639206−3.650.655701−1.170.662707−0.11300.6352610.619653−2.460.6302−0.80.634827−0.07400.6191440.607677−1.850.615421−0.60.61886−0.05Error range−22.09∼−1.85−2.46∼5.92−1.48∼−0.05p0.8nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)31.5283691.412346−7.59——1.507871−1.3441.3432641.267934−5.611.4582188.561.333216−0.7551.2513791.1953−4.481.248493−0.231.245424−0.4861.194891.150102−3.751.173717−1.771.190894−0.3371.155921.118601−3.231.131398−2.121.153027−0.2581.12711.09512−2.841.102841−2.151.12485−0.291.1045951.076739−2.521.081717−2.071.102872−0.16101.08651.061783−2.271.06519−1.961.085136−0.13201.000390.988779−1.160.988889−1.150.999985−0.04300.9670350.959519−0.780.959368−0.790.966858−0.02400.9482610.942748−0.580.942525−0.60.948161−0.01Error range−7.59∼−0.58−2.15∼8.56−1.34∼−0.01p0.9nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)32.2199992.162173−2.6——2.190339−1.3441.9439071.912762−1.62.1423510.211.931962−0.6151.812241.790519−1.21.8199840.431.806181−0.3361.7331631.716398−0.971.708671−1.411.729754−0.271.6795321.665911−0.811.647691−1.91.67754−0.1281.6403741.628914−0.71.60763−21.63918−0.0791.6103641.600302−0.621.578625−1.971.609567−0.05101.5863511.577401−0.561.556317−1.891.585873−0.03201.4745041.470441−0.281.457608−1.151.4748110.02301.4326151.429924−0.191.421102−0.81.4328840.02401.4093061.407292−0.141.400619−0.621.4095560.02Error range−2.6∼−0.14−2∼10.21−1.34∼0.02p0.95nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)32.80898212.769737−1.4——2.769095−1.4242.45302262.438936−0.572.71734710.82.437778−0.6252.28583262.278508−0.322.3009630.662.278716−0.3162.1867872.181884−0.222.158486−1.292.183009−0.1772.1200932.116571−0.172.081324−1.832.118108−0.0982.0716572.069042−0.132.031161−1.952.070709−0.0592.0347522.032585−0.111.995165−1.952.034297−0.02102.00523482.003434−0.091.967689−1.872.005280201.87020061.869459−0.041.848538−1.161.8708550.03301.8202881.819858−0.021.805492−0.811.8208760.03401.79274411.7925−0.011.781563−0.621.7932720.03Error range−1.4∼−0.01−1.95∼10.8−1.42∼0.03Table 4
Corollary results of coefficientk (C = 0.8).
p0.7nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)31.4414591.198475−16.862.24055555.441.392952−3.3741.1990031.056408−11.891.2330372.841.171917−2.2651.0796990.980532−9.181.065791−1.291.06182−1.6661.0062620.931047−7.470.984826−2.130.993137−1.370.9553580.895292−6.290.933649−2.270.945115−1.0780.9174790.867706−5.420.897132−2.220.909135−0.9190.8878130.845589−4.760.869215−2.090.880892−0.78100.8639640.827366−4.240.846899−1.980.857968−0.69200.7481680.733359−1.980.739788−1.120.745878−0.31300.702540.693546−1.280.697122−0.770.70116−0.2400.6765690.670218−0.940.672596−0.590.675594−0.14Error range−16.86∼−0.94−2.27∼55.44−3.37∼−0.14p0.8nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)32.0159631.864198−7.533.35025366.191.959245−2.8141.674761.595514−4.731.7562944.871.644344−1.8251.5138341.461578−3.451.508846−0.331.494208−1.361.4172331.378831−2.711.395374−1.541.403193−0.9971.351641.321513−2.231.326381−1.871.340847−0.881.3035231.278875−1.891.278522−1.921.294852−0.6791.2663891.245617−1.641.242699−1.871.259186−0.57101.2365721.218831−1.431.214529−1.781.230522−0.49201.0963321.089234−0.651.084534−1.081.094046−0.21301.0426141.038382−0.411.03481−0.751.041283−0.13401.0124751.009494−0.291.00666−0.571.011538−0.09Error range−7.53∼−0.29−1.92∼66.19−2.81∼−0.09p0.9nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)32.8704112.734096−4.754.9506872.472.800163−2.4542.3724812.313091−2.52.5208456.252.337578−1.4752.1452112.11023−1.632.1532160.372.124036−0.9962.0118881.988188−1.181.989734−1.11.997447−0.7271.9229261.905291−0.921.892856−1.561.912164−0.5681.8584961.844579−0.751.827017−1.691.850058−0.4591.8091821.797828−0.631.778539−1.691.802401−0.37101.7701821.760544−0.541.740922−1.651.764428−0.33201.589781.586399−0.211.573196−1.041.587908−0.12301.5226431.520797−0.121.511473−0.731.521655−0.06401.4855041.484164−0.091.477065−0.571.484813−0.05Error range−4.75∼−0.09−1.69∼72.47−2.45∼−0.05p0.95nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)33.6038363.451824−4.226.29378674.643.519022−2.3542.9680972.908834−23.1698196.82.928243−1.3452.6826752.65053−1.22.7001910.652.659239−0.8762.51697892.496323−0.822.493566−0.932.501342−0.6272.40694642.392608−0.62.372342−1.442.395769−0.4682.32787842.31704−0.472.29066−1.62.319351−0.3792.26777792.259295−0.372.230948−1.622.261001−0.3102.2202372.213298−0.312.184893−1.592.214701−0.25202.0035162.001622−0.091.982907−1.032.001987−0.08301.9241181.92322−0.051.910046−0.731.923348−0.04401.8804211.879741−0.041.869768−0.571.879932−0.03Error range−4.22∼−0.04−1.62∼74.64−2.35∼−0.03Table 5
Corollary results of coefficientk (C = 0.9).
p0.7nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)32.2277271.95304−12.333.92734576.292.172989−2.4641.6926641.592182−5.941.8186627.441.627362−3.8651.4554461.409432−3.161.4682630.881.406016−3.461.3172331.294781−1.71.30651−0.811.279076−2.971.2250711.214403−0.871.208123−1.381.19427−2.5181.158051.154134−0.341.139942−1.561.132505−2.2191.1064821.1067850.031.088974−1.581.084946−1.95101.0655161.0683620.271.048945−1.561.046877−1.75200.8736740.881020.840.86496−10.866203−0.86300.8005920.8069510.790.794931−0.710.79604−0.57400.7594730.7649010.710.755368−0.540.756278−0.42Error range−12.33∼0.84−1.58∼76.29−3.86∼−0.42p0.8nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)33.0391172.738367−9.95.81481491.333.03893−0.0142.2945402.186456−4.712.53491410.482.236538−2.5351.9756911.92394−2.622.0233452.411.927019−2.4661.7944521.766061−1.581.7972660.161.755281−2.1871.6753291.658734−0.991.664035−0.671.643227−1.9281.5898811.580105−0.611.573872−1.011.563058−1.6991.5252381.519469−0.381.507702−1.151.502188−1.51101.4740961.470974−0.211.456488−1.191.454012−1.36201.2405911.2432430.211.229623−0.881.232126−0.68301.1541521.156890.241.146775−0.641.148946−0.45401.1062531.1086810.221.100748−0.51.102528−0.34Error range−9.9∼0.24−1.19∼91.33−2.53∼−0.01p0.9nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)34.2580373.827763−10.118.552616101.864.3405671.9443.1876423.026395−5.063.59473912.773.144705−1.3552.7422352.658306−3.062.8423163.652.698114−1.6162.4934912.442339−2.052.5181360.992.456109−1.572.3323152.298158−1.462.331059−0.052.301023−1.3482.2183582.193841−1.112.206595−0.532.191628−1.292.1327052.114345−0.862.116515−0.762.109516−1.09102.0655742.051333−0.692.047597−0.872.045142−0.99201.7651611.762608−0.141.751624−0.771.756306−0.5301.6569911.656266−0.041.647536−0.571.651498−0.33401.5978421.597666−0.011.590631−0.451.593871−0.25Error range−10.11∼−0.01−0.87∼101.86−1.61∼1.94p0.95nExact solutionMethod 1Relative error (%)Method 2Relative error (%)Method 3Relative error (%)35.3111584.75345−10.510.85606104.45.4574182.7543.9558873.743654−5.364.500213.763.924381−0.853.3996313.286259−3.333.5430874.223.35901−1.1963.09160043.020273−2.313.1344261.393.055638−1.1672.89344862.844041−1.712.900530.242.862743−1.0682.75407452.717418−1.332.746018−0.292.727541−0.9692.64980592.62142−1.072.634866−0.562.626593−0.88102.5680892.545597−0.882.550266−0.692.547805−0.79202.2076162.201785−0.262.19221−0.72.19873−0.4302.0797232.077018−0.132.06865−0.532.074132−0.27402.0101912.008534−0.082.00167−0.422.006152−0.2Error range−10.5∼−0.08−0.7∼104.4−1.19∼2.75In order to better research the trend of the coefficientk corollary results, the relative error scatter plots of the coefficient k corollary results can be drawn when C = 0.5, C = 0.6, C = 0.7, C = 0.8, and C = 0.9, respectively, as shown in Figures 1–5. Since the relative error of Method 2 when n < 5 is too large, this part of the data is deleted when the figure is made.Figure 1
Relative error (C = 0.5).Figure 2
Relative error (C = 0.6).Figure 3
Relative error (C = 0.7).Figure 4
Relative error (C = 0.8).Figure 5
Relative error (C = 0.9).From the results of comparative analysis, the error range of method 1 is −30.59%∼0.84%, the error range of method 2 is −28.36%∼104.4%, and the error range of method 3 is −3.86%∼2.75%. Separately, whenp = 0.7, the relative error ranges of these three methods are −30.59%∼0.84%, −21.6%∼76.29%, and −3.86%∼1.67%; when p = 0.8, the relative error ranges of these three methods are −9.9%∼0.24%, −28.36%∼91.33%, and −2.81%∼2.0%; when p = 0.9, the relative error ranges of these three methods are −10.11%∼−0.01%, −14.5∼101.86%, and −2.45%∼2.11%; and when p = 0.95, the relative error ranges of these three methods are −10.5%∼0.85%, −5.02%∼104.4%, and −2.35%∼2.75%. It can be obtained that the error range calculated by method 2 is the largest, method 1 is in the middle, and method 3 is the smallest, indicating that the relative error calculated by method 3 is the smallest with the highest precision.WhenC = 0.9, p = 0.7∼0.95, and n = 3∼40, the absolute values of the average and maximum relative errors of method 1 are 2.28%, 12.33%, the absolute values of the average and maximum relative errors of method 2 are 10.33%, 104.4%, and the absolute values of the average and maximum relative errors of method 3 are 1.38%, 3.9%. It can be seen that the precision of method 2 is the lowest, method 1 is in the middle, and method 3 has the highest precision.With an absolute value of relative error not exceeding 5%, the applicable range of the method 1 isn > 10, and the applicable range of method 2 is n > 7, and method 3 satisfies all the values of n in the above study. When n < 5, the difference of the relative error calculated by method 2 is the largest, indicating that method 2 is more suitable for the approximate solution of the coefficient k when n ≥ 5. It can be seen from the figure that when p is increased from 0.7 to 0.95, the absolute value of relative error calculated by these three methods has a decreasing trend, indicating that the precision of these three methods is improved with the increase of p. With the increase of n, the absolute value of the relative error calculated by these three methods has a decreasing trend, indicating that the precision of these three methods is improved with the increase of n. When n > 10, the relative error trend curve of the coefficient k calculated by these three methods tends to be gentle, indicating that the accuracy of these three methods is improved with the increase of n. And the relative error calculated by method 3 fluctuates around 0, and there is no large variability. In summary, method 3 has a wider range of applications than the first two methods.


## 5. Concluding Remarks

Through the study of the inference method of the noncentral t-distribution value, the following main conclusions can be obtained:(1)
By comparing and analyzing the error range of coefficientk under different conditions, it is shown that the quantile value corollary method proposed in this paper is more accurate than the current approximation method.(2)
By comparing the influence of different number of samplesn on the coefficient k, it indicates that the proposed method in this paper has a wider scope of application.(3)
The method of inferring the quantile value proposed in this paper does not need to repeatedly check the value table and does not need to perform linear interpolation calculation, which simplifies the calculation process and has better convenience.



---

*Source: 9917575-2021-05-25.xml*