search for: trialid

Displaying 1 result from an estimated 1 matches for "trialid".

2013 Jan 28
2
Why are the number of coefficients varying? [mgcv][gam]
...10.5) htt.knots = c(350,450,550,650) bc.prc.C.knots = c(.3,.45,.6,.8) phi.knots = c(4.5,5.5,6.5,7.5) CEC.knots = c(5,12,19,26) soc.knots = c(10,20,30,40) sand.knots = c(.2,.4,.6,.8) clay.knots = c(.15,.3,.45,.6) abslat.knots = c(10,20,30,45) lon.knots = c(-50,0,50,125) dum = as.vector(rep(1,length(trialid))) doyee = NA r.ints = tryCatch(gam(RR ~ as.factor(pot.trial) + as.factor(year) + as.factor(crop.legume) + as.factor(crop.fruit) + as.factor(feedstock) + s(trialid,bs="re",by=dum) + s(BCAR.imp,bs...