[R] gss package - predict.ssanova

Chong Gu chong at stat.purdue.edu
Tue May 8 16:55:58 CEST 2001


   Sender: ernesto at neptuno.ipimar.pt
   Date: Tue, 08 May 2001 09:51:08 +0100
   From: =?iso-8859-1?Q?Jos=E9?= Ernesto Jardim <ernesto at ipimar.pt>
   Organization: IPIMAR
   X-Mailer: Mozilla 4.72 [en] (X11; I; Linux 2.2.14 i686)
   X-Accept-Language: en
   MIME-Version: 1.0
   CC: R-help at stat.math.ethz.ch
   References: <200105071639.LAA88280 at odds.stat.purdue.edu>
   Content-Type: text/plain; charset=us-ascii
   Content-Transfer-Encoding: 7bit

   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



recn.loc should be a matrix, not a data frame.
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list