Skip to content
Prev 378955 / 398502 Next
Original Raw

Obtaining individual parameters from a sitar growth model.

Rayt Chiruka 
I want to obtain the individual biological parameters (ie age at takeoff,
final height velocity at takeoff) after a sitar model. (I do understand a
sitar model fits a mean curve for the population).



I have managed to obtain the parameters for one individual using the code
below,





library(sitar)

data(heights)

library(tidyverse)

m1 <- sitar(x=age, y=height, id=id, data=heights, df=5)



#velocity at takeoff for individual 1



(vel_attoff<-getPeakTrough(plot_v(m1, subset=age < 14 & id==1) ,
peak=FALSE))





What I want is to produce a dataframe with the id and the velocity at
takeoff for each individual in the dataset.



I have tried to do this



for (i in seq_along(heights[1])){

                print(getPeakTrough(plot_v(m1, subset=age < 14 & id==i) ,
peak=FALSE))

}



But it prints only one x and y  value



Any assistance will be greatly appreciated.