# NOT RUN {
r1 <- gamlss(R~nn(~Fl+A+H+loc, size=10, decay=0.2), data=rent,
family=GA, gd.tol=1000, n.cyc=5)
getSmo(r1)
plot(getSmo(r1), y.lab=expression(eta[1]))
plot(getSmo(r1), y.lab=expression(g[1](mu)))
# }
# NOT RUN {
r2 <- gamlss(R~nn(~Fl+A+H+loc, size=10, decay=0.2),
sigma.fo=~nn(~Fl+A+H+loc, size=10, decay=0.2),data=rent,
family=GA, gd.tol=1000, n.cyc=5)
plot(getSmo(r2), y.lab=expression(g[1](mu)))
plot(getSmo(r2, what="sigma"), y.lab=expression(g[2](sigma)))
# }
Run the code above in your browser using DataLab