Я создаю центильные кривые, используя функцию lms в R. Код:
Value <- data[data$Patient.Age>50,i]
age <- data[data$Patient.Age>50,]$Patient.Age
par(col.main='white')
ct = lms(Value, age, families = "BCCGo",
k=7, cent=c(5,15,25,35,45,50,55,65,75,85,95),
mu.df = 1, sigma.df = 1)
par(col.main = 'black')
title(nl_names[i])
График с кривыми от 5 до 95 процентилей:
Я хочу извлечь эти кривые центилей в файл .mat, включая все точки данных из кривых сантили. Я не мог найти способ извлечь эту информацию.