這個問題是關於在R中基於逐點分位數估計值繪製同時空頻帶的簡便方法。假設我獲得了一個p值向量,我想通過使用uniform QQ情節,如果它的分佈接近均勻分佈(...如果不是,那麼哪個點「脫穎而出」,即落在該頻帶之外)。下面是一個例子(其中u_obs
表示p值從均勻分佈不一定 - 僅用於說明的目的):R中統一的QQ空間統一空頻帶
library(gap)
u_obs = qqrunif(1000)
qqunif(u_obs,pch=21,bg="blue",bty="n")
根據第5頁中討論的非參數的方法,圖1(左下小圖)本文:http://stat.wharton.upenn.edu/~buja/PAPERS/paper-sim.pdf,如果期望90%的同時界限,我們可以進行10,000次模擬並獲得5%的下限和上限。
我的問題是:
是否存在的R包/代碼段,能夠有效地執行這個任務?我認爲這屬於「功能數據分析」領域,但我不太確定。
我們可以從統一的QQ空間的同時空界限做出什麼推論?檢查那些不屬於樂隊的點?如果p值的真實分佈是均勻的(即來自零值),我預計有多少點可以合理地超出頻帶?
我不知道它是否是同時或不是,但你應該看看'汽車:: qqPlot()'和'mgcv :: qq.gam()' –