R version 3.1.0 (2014-04-10) -- "Spring DanceĀ. Platform:
x86_64-apple-darwin13.1.0 (64-bit).
I would like to convert the following function from nlme to nlmer. I am
finding it difficult to apply the documentation I can find on this procedure to
this problem. The function predicts earthquake ground-shaking (lnIa) as a
function of M, Rrup, Ztor, Fevent, Vs30 and e1. I am converting to nlmer, as
once I have it running, I will explore crossed and nested mixed effects.
eq1 <- nlme( lnIa ~ c0 + m1*(M-5) + (r1a + r1b*M)*log(sqrt(Rrup^2 + Ztor^2)) +
f1*Fevent + v1*log(Vs30/760) +
e1,fixed=c0+m1+r1a+r1b+f1+v1~1,random=e1~1|Event,data = data, start =
c(c0=11,m1=-0.9,r1a=0.5,r1b=-0.52,f1=0.9,v1=-1.32))
Many thanks for your help,
Roxane.
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
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.
