Вычислить процентили для распределения LMS R GAMLSS

Question

Можно ли использовать GAMLSS наоборот? То есть мне дали точки LMS диаграмм роста, и я хочу рассчитать процентили для него.

Данные в настоящее время в следующем формате:

AGE, L, M, S

Answer 1

вам не нужно GAMLSS для этого:

Я думаю, что либо вы имеете в виду что-то вроде этого

справочной таблицы RT содержит следующие colums:

возраст, L, M, S

и чем у вас есть определенное значение в данном возрасте, например. 100 см в высоту в возрасте 4

тогда вам нужна функция, как этот один

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-й процентиль

Тогда вам нужна функция, как это (л, м, с также может быть векторами, поэтому я думаю, что это удобно для вас цели):

percval <- function(l,m,s,perc){ (qnorm(perc)*l*s+1)**(1/l) * m }

Тогда вы получите значение 90-го процентиля через