Displaying 3 results from an estimated 3 matches for "fitted_model".
2010 Feb 16
2
HELP on Non-Linera Mixed-Effect model
Hi,
I'm trying to fit nonlinear mixed effects model using nlme function but getting an error message. Here is what I have:
fitted_model = nlme(scores~spline(b1,b2,b3,kt,time),
fixed = list(b1~1, b2~1, b3~1, kt~1),
random = b1+b2+b3~1,
groups= ~id,
data = sdat,
start = c(b1=3.5,b2=2,b3=.60,kt=3.5),verbose=T)
Error:
Error in spline(b1, b2, b3, kt, time) : invalid interpolati...
2007 Aug 02
1
Xyplot - adding model lines to plotted points
...possible.
I have trawled R help, and it seems that panel.superpose may be one
method of attempting this, however, I am unable to produce any working
code.
My dataset is too large to put in as an example, but it is in the basic
form as below (where sex is 1 for male and 5 for female).
len age fitted_model cohort sex
1 24 2 30.05771 1977 1
2 31 3 36.64122 1977 1
3 22 2 27.73938 1978 1
4 34 3 36.64122 1977 1
5 22 2 27.73938 1978 1
6 31 3 36.64122 1977 1
7 34 3 36.64122 1977 1
8 28 2 27.73938 1978 1
9 23 2 27.73938 1978 1
10 24...
2008 Jun 05
1
R-code embedded in VBE -- Type mismatch errors
...= summary(lm_output)[9]
# For a complete output:
# type summary(lm_output)
param_vec = as.numeric(c(a,se_a,b,se_slope,R2))
# Fitted curve
xx = c(1:length(lr))
yy = intercept+slope*xx
fitted_model = exp(yy)+1
# Add a tail factor if max_settle_time > no_of_years
if (max_settle_time > no_of_years){
xx_ext = c(no_of_years:max_settle_time)
yy_ext = intercept + slope*xx_ext
}...