The smooth.sspline function in sspline v0.1-5 is giving different results depending whether I run it in 32-bit or 64-bit platform. Correct results (TRUE result below) occur on the 32-bit platform. I have run it on Windows 7 x32 +R2.15.0 (works) and R2.14.2 (works), Windows 7 x64 +R2.15.0 (does not work) and R2.14.2 (does not work), and Ubuntu 12.04 x64 +R2.15.0 (does not work). "Does not work" means it gives results that are sometimes off by dozens of orders of magnitude (and sometimes the sign). A test case is given below.

The library uses Fortran subroutines, and I am not familiar with interfacing Fortran with R, so I hope there is someone who is familiar and has an x64 platform to debug with can take a look at it.

library(sspline)
dlon <- 125.4637
dlat <- 6.93252
grd <- structure(list(lonr = c(125.25, 125.25, 125.25, 125.25, 125.25,
125.333333333333, 125.333333333333, 125.333333333333, 125.333333333333,
125.333333333333, 125.416666666667, 125.416666666667, 125.416666666667,
125.416666666667, 125.416666666667, 125.5, 125.5, 125.5, 125.5,
125.5, 125.583333333333, 125.583333333333, 125.583333333333,
125.583333333333, 125.583333333333), latr = c(6.75, 6.83333333333333,
6.91666666666667, 7, 7.08333333333333, 6.75, 6.83333333333333,
6.91666666666667, 7, 7.08333333333333, 6.75, 6.83333333333333,
6.91666666666667, 7, 7.08333333333333, 6.75, 6.83333333333333,
6.91666666666667, 7, 7.08333333333333, 6.75, 6.83333333333333,
6.91666666666667, 7, 7.08333333333333), value = c(3.75, 3.65,
3.7, 3.7, 3.75, 3.75, 3.45, 3.6, 3.55, 3.65, 3.95, 3.7, 3.75,
3.75, 3.8, 3.95, 4.1, 3.8, 3.85, 3.8, 3.95, 4, 4.05, 3.95, 3.85
)), .Names = c("lonr", "latr", "value"), row.names = c(271L,
272L, 273L, 274L, 275L, 297L, 298L, 299L, 300L, 301L, 323L, 324L,
325L, 326L, 327L, 349L, 350L, 351L, 352L, 353L, 375L, 376L, 377L,
378L, 379L), class = "data.frame")

linkespline <- smooth.sspline( grd$lonr, grd$latr, grd$value )
all.equal( 3.800515, predict.smooth.sspline( linkespline, dlon, dlat ), tolerance=0.001 )

---------------------------------------------------------------------------
Jeff Newmiller                        The     .....       .....  Go Live...
DCN:<jdnew...@dcn.davis.ca.us>        Basics: ##.#.       ##.#.  Live Go...
                                      Live:   OO#.. Dead: OO#..  Playing
Research Engineer (Solar/Batteries            O.O#.       #.O#.  with
/Software/Embedded Controllers)               .OO#.       .OO#.  rocks...1k

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to