[R] Problem with "different" proj4string that are actually identical

MacQueen, Don macqueen1 at llnl.gov
Sat Mar 21 00:18:03 CET 2015


Hmm, that does complicate things. Since the error message references
"join(pts, ta)" I think you're going to have to find out more about this
join() function and how niche.test() constructs the arguments that it
passes to join().

Starting with
  traceback() 
immediately after getting the error might help.

If you move this question over to r-sig-geo your chances of getting help
from someone familiar with niche.test() and/or join() is better.

-- 
Don MacQueen

Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062





On 3/20/15, 11:24 AM, "Hemery, Lenaig" <Lenaig.Hemery at oregonstate.edu>
wrote:

>Dear Don,
>Thanks for your answer.
>proj4string(spCro) == proj4string(paramdf), as well as
>identical(proj4string(spCro), proj4string(paramdf)), returns a TRUE so I
>forced them to be identical as you suggested, and repeated the
>niche.test(), but I got the same error message.
>Would you have any other ideas?
>Thanks a lot,
>Lenaig
>________________________________________
>From: MacQueen, Don [macqueen1 at llnl.gov]
>Sent: Friday, March 20, 2015 8:23 AM
>To: Hemery, Lenaig; r-help at R-project.org
>Subject: Re: [R] Problem with "different" proj4string that are actually
>identical
>
>This question would be better asked on r-sig-geo.
>
>Here are some things to try (type at the R prompt):
>
>  proj4string(spCro) == proj4string(paramdf)
>
>Does it return TRUE or FALSE?
>
>Or, instead of using summary(), directly inspect the proj4strings:
>
>  proj4string(spCro)
>  proj4string(paramdf)
>
>If there still appears to be no difference between them,
>try forcing them to be identical
>  proj4string(spCro) <- proj4string(paramdf)
>
>then repeat your niche.test()
>
>
>--
>Don MacQueen
>
>Lawrence Livermore National Laboratory
>7000 East Ave., L-627
>Livermore, CA 94550
>925-423-1062
>
>
>
>
>
>On 3/19/15, 5:29 PM, "Hemery, Lenaig" <Lenaig.Hemery at oregonstate.edu>
>wrote:
>
>>Dear everybody
>>I am trying to run an Ecological Niche Factor Analysis (enfa) with the
>>package "adehabitatHS" but before I do so, I want to test the
>>significance of its parameters marginality and tolerance, by using the
>>function niche.test() that runs a Monte Carlo test. I used to do that
>>with the old package "adehabitat" that is now deprecated so I updated all
>>my codes. Everything I want to do seems to work fine except the
>>niche.test() function and I get the error message here below.
>>
>>> testCro<-niche.test(paramdf, spCro, nrep = 999, o.include = TRUE)
>>Error in join(pts, ta) : different proj4string in x and xy
>>
>>I checked the proj4string for x (paramdf) and xy (spCro) and there are
>>exactly identical (see below).
>>
>>> summary(spCro)
>>Object of class SpatialPoints
>>Coordinates:
>>              min     max
>>longitude -124.91 -124.13
>>latitude    39.86   46.95
>>Is projected: FALSE
>>proj4string :
>>[+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0]
>>Number of points: 84
>>
>>> summary(paramdf)
>>Object of class SpatialPixelsDataFrame
>>Coordinates:
>>         min       max
>>x -125.50000 -123.7998
>>y   39.50384   47.0000
>>Is projected: FALSE
>>proj4string :
>>[+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0]
>>Number of points: 108237
>>Grid attributes:
>>  cellcentre.offset cellsize cells.dim
>>x        -125.49517  0.00966       176
>>y          39.50867  0.00966       776
>>Data attributes:
>>      bath              slope              rugo               sumT
>> Min.   :-3611.83   Min.   : 0.0000   Min.   :0.000000   Min.   : 1.184
>> 1st Qu.:-1885.11   1st Qu.: 0.5271   1st Qu.:0.008044   1st Qu.: 2.150
>> Median : -894.61   Median : 1.4974   Median :0.036921   Median : 3.978
>> Mean   :-1157.85   Mean   : 3.4237   Mean   :0.147369   Mean   : 4.794
>> 3rd Qu.: -150.85   3rd Qu.: 4.1152   3rd Qu.:0.197974   3rd Qu.: 7.388
>> Max.   :   -0.05   Max.   :61.3592   Max.   :0.992856   Max.   :14.436
>>[...]
>>
>>I think I looked for all the help I could find on the internet and didn't
>>find any solution so here is my last chance.
>>Does anybody have an idea on how to solve this problem?
>>By the way, I use the R version 3.1.3 and up-to-date "adehabitatHS",
>>"ade4", "sp" and other packages.
>>Thanks a lot,
>>Lenaig
>>
>>‰‰‰‰ ><)))°> <°(((>< ‰‰‰‰
>>Lenaïg Hemery, Ph. D.
>>Research Associate - post-doc
>>Hatfield Marine Science Center
>>Oregon State University
>>2030 S.E. Marine Science Drive, Newport, OR 97365, USA
>>Cell:+1-541-272-7196
>>______________________________________________
>>R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>>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.
>



More information about the R-help mailing list