Displaying 2 results from an estimated 2 matches for "error_i".
Did you mean:
error_r
2010 Jul 21
0
Piecewise regression using lme()
...is:
y_t=
(beta_01+beta_11*t+error_1)*(group==1)+(beta_02+beta_12*t+error_2)*(group==2)+(beta_03+beta_13*t+error_3)*(group==3)
It looks like a mixed effects model were the fixed effect are the piecewise
linear regression parts (beta_0i+beta_1i*t) and the random effects are the
variance components error_i.
The problem is that I can't find the way the write this model correctly
using the lme() function of the package nlme. I can't remove the global
intercept and the global variance component. Here's what I tried:
#1 Using a prior piecewise linear regression
lm.list=lmList(y~t|group,Da...
2007 Dec 13
1
Probelms on using gam(mgcv)
Dear all,
Following the help from gam(mgcv) help page, i tried to analyze my
dataset with all the default arguments. Unfortunately, it can't be run
successfully. I list the errors below.
#m.gam<-gam(mark~s(x,y)+s(lstday2004)+s(slope)+s(ndvi2004)+s(elevation)+s(disbinary),family=binomial(logit),data=point)