1

是否可以反过来使用 GAMLSS？也就是说，我已经获得了增长图表的 LMS 点，我想计算它的百分位数。

数据目前采用以下格式：

年龄，L，M，S

4

1 回答 1

3

为此，您不需要 GAMLSS：

我想要么你的意思是这样的

参考表 rt 包含以下列：

年龄，l，m，s

并且比你在给定年龄有一些价值，例如 4 岁时身高 100 厘米

那么你需要一个像这样的功能

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 回答