Я пытаюсь создать модель для оценки центилей переменной угла по возрастной переменной (age). Однако я не могу запустить следующую модель, потому что я не знаю, как сгенерировать p [1], и я не знаю, что означает p [1, 2 или 3] в примере. Кто-нибудь может мне помочь?
Tage<-(Homens$idade)^(m0$power)
mOD <- gamlss(angulo~cs(Tage,df=op$par[1]),
sigma.formula=~cs(Tage,df=op$par[2]),
nu.formula=~cs(Tage,df=op$par[3]),c.spar=(-1.5,2.5), family=BCCGo,data=Homens)