Here is the information about this objects. > recs.spliney Call: ssanova(formula = recs.mds[[1]][, 2] ~ recs.loc$x * recs.loc$y, type = "tp") Terms: [1] "1" "recs.loc$x" "recs.loc$y" [4] "recs.loc$x:recs.loc$y" Number of fixed and random effects: Fixed Random 1 1 0 recs.loc$x 1 1 recs.loc$y 1 1 recs.loc$x:recs.loc$y 1 3 total 4 5 Smoothing parameters are selected by GCV. > is.data.frame(recn.grid) [1] TRUE > names(recn.grid) [1] "x" "y" Regards EJ The variables in the fit are "recs.loc$x" and "recs.loc$y", but the variables in recn.grid are "x" and "y". BTW, are recs.loc$x and recs.loc$y the x-y coordinate of geographic locations? If the answer is yes, then you probably want to enter it as a single variable with two columns. What you are now doing appears to fit a tensor product spline. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
gss package - predict.ssanova
3 messages · Chong Gu, Ernesto Jardim
Chong Gu wrote:
Here is the information about this objects.
> recs.spliney
Call:
ssanova(formula = recs.mds[[1]][, 2] ~ recs.loc$x * recs.loc$y, type = "tp")
Terms:
[1] "1" "recs.loc$x" "recs.loc$y"
[4] "recs.loc$x:recs.loc$y"
Number of fixed and random effects:
Fixed Random
1 1 0
recs.loc$x 1 1
recs.loc$y 1 1
recs.loc$x:recs.loc$y 1 3
total 4 5
Smoothing parameters are selected by GCV.
> is.data.frame(recn.grid)
[1] TRUE
> names(recn.grid)
[1] "x" "y" Regards EJ The variables in the fit are "recs.loc$x" and "recs.loc$y", but the variables in recn.grid are "x" and "y". BTW, are recs.loc$x and recs.loc$y the x-y coordinate of geographic locations? If the answer is yes, then you probably want to enter it as a single variable with two columns. What you are now doing appears to fit a tensor product spline.
Thanks I'll take a look at what you're saying. Tomorrow there will be news :-) Regards EJ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Chong Gu wrote:
Here is the information about this objects.
> recs.spliney
Call:
ssanova(formula = recs.mds[[1]][, 2] ~ recs.loc$x * recs.loc$y, type = "tp")
Terms:
[1] "1" "recs.loc$x" "recs.loc$y"
[4] "recs.loc$x:recs.loc$y"
Number of fixed and random effects:
Fixed Random
1 1 0
recs.loc$x 1 1
recs.loc$y 1 1
recs.loc$x:recs.loc$y 1 3
total 4 5
Smoothing parameters are selected by GCV.
> is.data.frame(recn.grid)
[1] TRUE
> names(recn.grid)
[1] "x" "y" Regards EJ The variables in the fit are "recs.loc$x" and "recs.loc$y", but the variables in recn.grid are "x" and "y". BTW, are recs.loc$x and recs.loc$y the x-y coordinate of geographic locations? If the answer is yes, then you probably want to enter it as a single variable with two columns. What you are now doing appears to fit a tensor product spline. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
Hi I've tried your proposal and I got an error message
ssanova(recn.mds[[1]][,1]~recn.loc,"tp")
Error in model.frame(formula, rownames, variables, varnames, extras, extranames, :
invalid variable type
The objects are in aggreement with your description
is.data.frame(recn.loc)
[1] TRUE
is.vector(recn.mds[[1]][,1])
[1] TRUE
names(recn.loc)
[1] "x" "y" Regards EJ -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._