[R] gss package - predict.ssanova

[José Ernesto Jardim]

Chong Gu wrote:

>    Date: Mon, 07 May 2001 12:08:34 +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
>    Content-Type: text/plain; charset=us-ascii
>    Content-Transfer-Encoding: 7bit
>    Sender: owner-r-help at stat.math.ethz.ch
>    Precedence: SfS-bulk
>
>    Hi
>
>    I'm using the gss package to fit thin plate splines. I've fitted the tps
>    without problems and got a ssanova object. Then I wanted to do some
>    prediction using a new set of data (latlon data) and I got an error
>    message:
>
>    > predict.ssanova(recs.spliney,recn.grid,se.fit=FALSE)
>    Error in array(x, c(length(x), 1), if (!is.null(names(x)))
>    list(names(x),  :
>            attempt to set an attribute on NULL
>
>    Can someone help me with this ?
>
>    Thanks
>
>    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
>    _._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
>
> Is recn.grid a data frame with the same variable names as used in the
> fitting formula?

Yes.

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

-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._