Displaying 2 results from an estimated 2 matches for "abslat".
Did you mean:
ablat
2013 Jan 28
2
Why are the number of coefficients varying? [mgcv][gam]
...BCAR.knots = c(2,15,60,120)
INAR.knots = c(50,100,200,300)
bcph.knots = c(7.5,8.5,9.5,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.fac...
2013 Mar 21
1
[mgcv][gam] Odd error: Error in PredictMat(object$smooth[[k]], data) : , `by' variable must be same dimension as smooth arguments
....prc.C.imp"){j=13}
if (varstochange[i] == "CEC.imp"){j=14}
if (varstochange[i] == "soc.imp"){j=15}
if (varstochange[i] == "sand.imp"){j=16}
if (varstochange[i] == "clay.imp"){j=17}
if (varstochange[i] == "abslat.imp"){j=18}
preddat[j] = val[i]
}
predict.gam(MI,newdata=preddat,se.fit=TRUE)
}
I then make predictions that look like this:
a = makepred(c("phi.imp","bcph.imp","year"),c(4.5,7.25,1))
b = makepred(c("phi.imp","bcph.imp&...