Good Afternoon Everyone, I am looking for advice fitting a linear mixed model where the random components do not seem to fit within the model formulae for lmer. The columns of Z are not stratified and have the notional random formula (z1 | 1) + ... + (zk | 1).
Context I am fitting a penalised thin plate spline with knots k1 to kn. The basis functions Zk are |x-ki|^3 and the penalty matrix has elements |ki-kj|^3. Rather than fit basis functions Zk I have transformed the basis to Z = Zk %*% penalty^(-1/2). Full details are in section 2: Bayesian analysis for penalised spline regression using winBUGS. https://www.jstatsoft.org/article/view/v014i14/v14i14.pdf I am specifically looking for an alternative to the function spm() in package semipar. Thank you for any help or advice. Dave. David Mcnulty Medical Statistician Tel: +44 (0) 121 371 2448 Internal: 12448 Email: david.mcnu...@uhb.nhs.uk Web: http://www.uhb.nhs.uk Health Informatics - University Hospitals Birmingham NHS Foundation Trust Yardley Court, 11-13 Frederick Road, Edgbaston Birmingham, B15 1JD The Health Informatics Department is a team of specialist NHS information analysts who work together to provide analysis solutions. The Team filters, analyses and presents comprehensive information that empowers NHS professionals to make informed decisions. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
