Dear Ivan
Dear community
Quite nice book recommendation.
Yes indeed my raster "s" (the shape file for the boxplot classes) has
several layers. That's way I tried to select a layer by "
s<-sf$Unterregio".
> sf <- read_sf("C:/Users/....._BiogeoRegion.shp")
> names(sf)
> names(sf)
[1] "RegionNumm" "RegionName" "Unterregio"
"Unterreg_1" "ObjNummer" "Version"
"Shape_Leng" "Shape_Area" "DERegionNa"
"FRRegionNa" "ITRegionNa" "DEBioBedeu"
"FRBioBedeu" "ITBioBedeu"
[15] "geometry"
> s<-sf$Unterregio
> r<-allrasters_pres[[1]]
Kind regards
Sibylle
-----Original Message-----
From: Ivan Krylov <ikrylov at disroot.org>
Sent: Friday, August 23, 2024 5:30 PM
To: sibylle.stoeckli at gmx.ch
Cc: 'SIBYLLE ST?CKLI via R-help' <r-help at r-project.org>
Subject: Re: [R] paired raster boxplots
? Fri, 23 Aug 2024 10:15:55 +0200
<sibylle.stoeckli at gmx.ch> ?????:
> > s<-sf$Unterregio
> > r<-allrasters_pres[[1]]
> >
> >
> > rs <- stack(r, s)
> > names(rs) <- c('r', 's')
> Error in `names<-`(`*tmp*`, value = c("r", "s")) :
> incorrect number of layer names
It looks like at least one of the rasters 'r' and 's' has
multiple layers. What does names(rs) return? I would offer more detailed advice,
but I don't know 'raster' that well.
The "R Inferno" book [1] offers a lot of generic-R troubleshooting
advice, which should help you progress past errors like this one without waiting
for someone on R-help to reply.
--
Best regards,
Ivan
[1] https://www.burns-stat.com/documents/books/the-r-inferno/