[R-sig-Geo] Error In proj4string

[alannie]

Hi list!

I have been trying to run some code that I have run tons of times before,
only to get a warning message:

##add in ASCII raster maps bio1-11
> bio1 <- readGDAL("bio1.asc") 
bio1.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio2 <- readGDAL("bio2.asc") 
bio2.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio3 <- readGDAL("bio3.asc") 
bio3.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio4 <- readGDAL("bio4.asc") 
bio4.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio5 <- readGDAL("bio5.asc") 
bio5.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio6 <- readGDAL("bio6.asc")
bio6.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio7 <- readGDAL("bio7.asc")  
bio7.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio8 <- readGDAL("bio8.asc")
bio8.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio9 <- readGDAL("bio9.asc")
bio9.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio10 <- readGDAL("bio10.asc")
bio10.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
> bio11 <- readGDAL("bio11.asc")
bio11.asc has GDAL driver AAIGrid 
and has 1486 rows and 1347 columns
###combine raster maps
> temp <- cbind("bio1", "bio2", "bio3",
> "bio4","bio5","bio6","bio7","bio8","bio9","bio10","bio11")
> names(temp) <- c("bio1", "bio2", "bio3",
> "bio4","bio5","bio6","bio7","bio8","bio9","bio10","bio11") 
###read in point file lin.shp
> lin <- readShapePoints("lin.shp")
###set proj4strings to the same proj4string: *this is where the issue is*
> proj4string(lin)= proj4string(temp)
*Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘proj4string’ for
signature ‘"matrix"’*
###what is in "temp" ?
> temp
     [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]   [,10]  
[,11]  
[1,] "bio1" "bio2" "bio3" "bio4" "bio5" "bio6" "bio7" "bio8" "bio9" "bio10"
"bio11"
attr(,"names")
 [1] "bio1"  "bio2"  "bio3"  "bio4"  "bio5"  "bio6"  "bio7"  "bio8"  "bio9" 
"bio10" "bio11"
> class(temp)
[1] "character"
###end

I am not sure why this issue is occurring. 
I have checked all the projections - all files are in the same projection.
I have restarted R - nothing gets better.
I have tried data for a different organism (it is the same data, but
belonging to a different species) - this does not work.
I have the following packages loaded:
library(maptools)
library(adehabitatHS)
library(SDMTools)
library(spatial)
library(sp)
library(rgdal)
library(dismo)
library(shapefiles)


Thanks, List!

 







--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/Error-In-proj4string-tp7587305.html
Sent from the R-sig-geo mailing list archive at Nabble.com.