Displaying 2 results from an estimated 2 matches for "feedstock".
2013 Jan 28
2
Why are the number of coefficients varying? [mgcv][gam]
...ts = 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="cr",k=length(BCAR.knots),by=as.factor(pot.trial))
+
s(INAR.imp,bs="cr",k=length(INAR.knots),by=as.factor(crop.legume))
+ s(bcph.imp,bs="cr",k=length(bcph.knots),by=as.factor...
2013 Mar 21
1
[mgcv][gam] Odd error: Error in PredictMat(object$smooth[[k]], data) : , `by' variable must be same dimension as smooth arguments
...:length(varstochange)){
if (varstochange[i] == "pot.trial"){j=1}
if (varstochange[i] == "year"){j=2}
if (varstochange[i] == "crop.legume"){j=3}
if (varstochange[i] == "crop.fruit"){j=4}
if (varstochange[i] == "feedstock"){j=5}
if (varstochange[i] == "BCAR.imp"){j=8}
if (varstochange[i] == "INAR.imp"){j=9}
if (varstochange[i] == "bcph.imp"){j=10}
if (varstochange[i] == "phi.imp"){j=11}
if (varstochange[i] == "htt.imp"...