Displaying 1 result from an estimated 1 matches for "b_12_i".
Did you mean:
b_12
2006 Nov 07
1
gamm(): nested tensor product smooths
...where X represents the unpenalized part of the spline functions and Z the "wiggly" parts, this would be:
y=X%*%beta+ Z_1%*%b_1+ Z_2%*%b_2
vs
y=X%*%beta+ Z_1%*%b_1+ Z_2%*%b_2 + Z_12 %*% b_12
where b are random effect vectors and the hypothesis to be tested is
H_0: Var(b_12)=0 (<=> b_12_i == 0 for all i)
the problem:
gamm() does not seem to support the use of nested tensor product splines,
does anybody know how to work around this?
example code: (you'll need to source the P-spline constructor from ?p.spline beforehand)
###########
test1<-function(x,z,sx=0.3,sz=0.4)
{...