我有以下分配:如何計算R中曲線下面積的95%可信極限?
x<-c(22.5,28.14285714,33.78571429,39.42857143,45.07142857,50.71428571,56.35714286,62,67.64285714,73.28571429,78.92857143,84.57142857,90.21428571,95.85714286,101.5,107.1428571,112.7857143,118.4285714,124.0714286,129.7142857,135.3571429,141,146.6428571,152.2857143,157.9285714,163.5714286,169.2142857,174.8571429,180.5,186.1428571,191.7857143,197.4285714,203.0714286,208.7142857,214.3571429,220,225.6428571,231.2857143,236.9285714,242.5714286,248.2142857,253.8571429,259.5,265.1428571,270.7857143,276.4285714,282.0714286,287.7142857,293.3571429,299)
y<-c(0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.003541755,0.00328839614285714,0.00296425985714286,0.002655899,0.00236187857142857,0.002080895,0.00181184271428571,0.00155376085714286,0.00130578928571429,0.001074706,0.000877193,0.000709397142857142,0.000567189714285714,0.000447254,0.000346858571428571,0.000263689142857143,0.000195768428571429,0.000141427,9.92657142857141e-05,6.77857142857142e-05,4.48571428571428e-05,2.86428571428571e-05,1.75142857142857e-05,1.01357142857143e-05,5.52e-06,2.78857142857142e-06,1.27285714285713e-06,5.00714285714284e-07,1.5742857142857e-07,3.29857142857142e-08,2.78857142857137e-09,1.74e-12)
plot(x,y)
我想找到的x
,分配到左和麪積的0.05到右側下分離的0.95的區域中的值(單尾95%信譽間隔)。
我想我必須將我的經驗曲線擬合到一個函數中,然後整合函數以便我可以獲得所需的值,但是我不知道從哪裏開始。
這怎麼可能在R中完成?
GSee的答案是要走的路。但我想指出,源數據的數值積分不僅比創建擬合函數和積分更容易,而且計算誤差也更小(一般而言)。 –
@CarlWitthoft,我不太確定我的答案(這是'quantile(x,0.95)')。它將'x'分成95%和5%,但根本不考慮區域('y')。 – GSee
@Gsee - 我想你應該生成辛普森積分的分位數。我仍然不會產生合適的功能,壽。 –