[R] Illegal Instruction

2006-02-28 Thread David Pleydell
I am running Debian Etch on a dual Xeon 64 bit Dell Precision with the latest versions of R and RandomFields. I am unable to run the first example in ?GaussRF (in RandomFields). The bug occurs when GaussRF runs the line error - InitSimulateRF(x = x, y = y, z = z, T = T, grid = grid, model =

[R] Re Illegal Instrunction

2006-02-28 Thread David Pleydell
OK I found the problem. locate libblas.so.3 showed there were two such files, one in /usr/lib/, the other in /usr/lib/atlas/ removing the atlas installation cured the bug. Not sure what I should do if I later find that I need atlas for some reason. cheers David

[R] corCompSymm in nlme package

2005-02-23 Thread David Pleydell
We are trying to use the corCompSymm function in nlme The example from the help pages for the corAR1 function gives the following corAR1(0.2, form = ~ 1 | Mare) Correlation structure of class corAR1 representing Phi 0.2 We are expecting a somewhat similar correlation specification with the

[R] strange results with dmvnorm

2004-10-31 Thread David Pleydell
I am experiencing strange results using dmvnorm. I define a scaled distance matrix from the coordinates bellow and then calculate a covariance matrix using a spherical correlation function. Then with certain combinations of range and sill parameters dmvnorm is returning values greater than 1.

Re: [R] non-linear trends in kriging model

2003-10-06 Thread David Pleydell
) fit you model by alternating between using nls() and passing the residuals to the functions in geoR or 2) try trying to fit using the functions in the nlme() package Regards P.J. On Tue, 30 Sep 2003, David Pleydell wrote: Hi I am struggling to fit a non-linear trend using the likfit function

[R] non-linear trends in kriging model

2003-09-30 Thread David Pleydell
Hi I am struggling to fit a non-linear trend using the likfit function in geoR. Specifically I want a sigmoidal function, something like SSfpl in the nls package to fit the trend. But it seems trend.spatial in geoR only works with lm or glm type models. Any ideas how I can specify the model to