GAMLSS を逆に使用することは可能ですか? つまり、成長チャートの LMS ポイントが与えられており、そのパーセンタイルを計算したいと考えています。
現在、データは次の形式になっています。
年齢、L、M、S
GAMLSS を逆に使用することは可能ですか? つまり、成長チャートの LMS ポイントが与えられており、そのパーセンタイルを計算したいと考えています。
現在、データは次の形式になっています。
年齢、L、M、S
これには GAMLSS は必要ありません。
私はあなたがこのようなことを意味していると思います
参照テーブル rt には、次の列が含まれています。
年齢、l、m、s
4 歳で身長 100 cm など、特定の年齢で何らかの価値があるよりも
次に、このような関数が必要です
sds <- function(age,value){
m <- approx(rt$age,rt$m,xout=age,rule=1)$y
l <- approx(rt$age,rt$l,xout=age,rule=1)$y
s <- approx(rt$age,rt$s,xout=age,rule=1)$y
((value/m)**l-1)/(l*s)
}
セックスを含める必要があるかもしれません - 関数をどのように使用するかによって異なります
興味があれば: 私はこれを行うためのパッケージを作成しました。これは、多くの異なる参照に対して簡単な方法です (いくつか追加するのも簡単です)。まだベータ版ですが、データ フレームを操作するためのより便利な方法を提供します。
または多分あなたはこれを意味しますか???:
テーブルしかなく、おそらく 90 パーセンタイルを描画したいとします。
次に、次のような関数が必要です(l、m、sはベクトルにすることもできるので、目的に便利だと思います):
percval <- function(l,m,s,perc){
(qnorm(perc)*l*s+1)**(1/l) * m
}
次に、次の方法で90パーセンタイルの値を取得します
percval(rt$l,rt$m,rt$s,0.9)