1

GAMLSS を逆に使用することは可能ですか? つまり、成長チャートの LMS ポイントが与えられており、そのパーセンタイルを計算したいと考えています。

現在、データは次の形式になっています。

年齢、L、M、S

4

1 に答える 1

3

これには 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)
于 2013-10-01T13:42:45.320 に答える