stat_rethinking_2023
stat_rethinking_2023 copied to clipboard
height and age as piecewise linear in 03_howell_new_weight_model.r not working
Let me start by expressing my gratitude for the great lecture and book. Thanks for making this available on youtube and GitHub.
I did notice, that the last few lines in https://github.com/rmcelreath/stat_rethinking_2023/blob/main/scripts/03_howell_new_weight_model.r are unfinished. There is a missing closing parenthesis as well as missing prior specifications.
#######
# height and age as piecewise linear
data(Howell1)
d <- Howell1
dat <- list(
H = d$height,
A = d$age )
m <- quap(
alist(
H ~ dnorm(mu,sigma),
mu <- a1*(1-exp(-b1*A)) + a2*(1-exp(-b2*(A))
)
)