ksプロットのデータの入手方法

Question

Gumbelフィットを使用して気象データのKolmogorov-Smirnovテストを行いました。

ks.gumbel(data$mm, 63.2899, 30.06080, alternative = "two.sided", plot = TRUE) 

そして、y = 0.996にプロットされた曲線のx値を取得したいとします。

アイデア？ありがとう！

Answer 1

ks.gumbelのソースコードを見ると、プロットはplot(ecdf(x))で生成されますが、プロットからの値は関数から返されません。従って、あなたが直接求める情報をks.gumbelから得ることは不可能ですが、私たちは自分自身の中に入ることができるほど十分に知っています。

https://stackoverflow.com/a/16819127/1017276の答えを使用し、unirootでこの知識を組み合わせることで、我々はyの所望の値にxの値を返すようにecdfによって返される関数を解くことができます。

uniroot(f = function(x) ecdf(data$mm)(v = x) - .996, 
     interval = range(data$mm, na.rm = TRUE)) 

unirootによって返されたリストのroot要素は、あなたが探している値が得られます。