search for: sibylle

Dear community Dear Ivan Thanks a lot. The code works now. Solution: direct and full path to the .tif files. I confused back and forward slash #first import all files in a single folder as a list rastlist <- list.files(path = "C:/Users/Sibylle St?ckli/Desktop/NCCS_Impacts_Lot2_2022/InVEST/Species_Input/valpar_bee_presence", pattern='.tif$', all.files= T, full.names= T) At the end I got a warning #to run a function on an individual raster e.g., plot > boxplot(allrasters[[1]]) > > p<-boxplot(allrasters[[1]]) War...

...7,16,3,1,2237,611,403,240,280,0,0), 2,7)) barplot(F1, beside = TRUE, col = c("blue", "grey")) axis(2, at=c(0,10,20,30,40,50,60, labels=c(0,10,20,30,40,50,60))) axis(4, at = c(0,500,1000,1500,2000,2500), labels = c(0,500,1000,1500,2000,2500)) -----Message d'origine----- De?: sibylle.stoeckli at gmx.ch <sibylle.stoeckli at gmx.ch> Envoy??: 24 novembre 2023 11:27 ??: 'Charles-?douard Gigu?re' <ce.giguere at gmail.com>; r-help at r-project.org Objet?: RE: [R] ggplot adjust two y-axis Dear Charles-Edouard Thanks a lot. Yes indeed barplot sounds excellent. U...

...le(matrix(c(50,11,6,17,16,3,1,2237,611,403,240,280,0,0), 2,7)) barplot(F1, beside = TRUE, col = c("blue", "grey")) axis(2, at=c(0,10,20,30,40,50,60, labels=c(0,10,20,30,40,50,60))) axis(4, at = c(0,500,1000,1500,2000,2500), labels = c(0,500,1000,1500,2000,2500)) Kind regards Sibylle -----Original Message----- From: Charles-?douard Gigu?re <ce.giguere at gmail.com> Sent: Friday, November 24, 2023 3:57 PM To: sibylle.stoeckli at gmx.ch; r-help at r-project.org Subject: RE: [R] ggplot adjust two y-axis Hi, I don't know the axis mecanism well enough in ggplot but u...

...lot function you can add an axis on the right using the axis function. Here is an example: test <- as.table(matrix(c(2,10,3,11), 2,2)) barplot(test, beside = TRUE, col = scales::brewer_pal(palette = 1)(2)) axis(4, at = c(0, 5, 10), labels = c(0,50,100)) -----Message d'origine----- De?: sibylle.stoeckli at gmx.ch <sibylle.stoeckli at gmx.ch> Envoy??: 24 novembre 2023 09:27 ??: 'Charles-?douard Gigu?re' <ce.giguere at gmail.com>; r-help at r-project.org Objet?: RE: [R] ggplot adjust two y-axis Dear Charles-Edouard Thanks a lot. So no way in R to just simply have one...

Dear Ivan Thanks a lot. I tried now to provide the full path. However probably the "?" in the path produces the error, would that be possible? > #first import all files in a single folder as a list > rastlist <- list.files(path = "C:\Users\Sibylle St?ckli\Desktop\NCCS_Impacts_Lot2_2022\InVEST\Species_Input\valpar_bee_presence", pattern='.tif$', all.files= T, full.names= T) Error: '\U' used without hex digits in character string (<input>:1:35) Kind regards Sibylle -----Original Message----- From: Ivan Krylov <i...

...non-missing values: all.equal(getValues(r1)[!mask], getValues(r2)[!mask]) --> what do you mean by use is.na() in getValues(). So I need to call getValues a second time? I suppose you mean to first prepare a mask using is.na without getValues and then in the second step your code? Kind regards Sibylle -----Original Message----- From: Ivan Krylov <ikrylov at disroot.org> Sent: Friday, August 16, 2024 9:28 AM To: SIBYLLE ST?CKLI via R-help <r-help at r-project.org> Cc: sibylle.stoeckli at gmx.ch Subject: Re: [R] allequal diff ? Fri, 16 Aug 2024 07:19:38 +0200 SIBYLLE ST?CKLI via R-h...

...ght? is not working. mylist<-list(p1, p2) dev.new(width=28, height=18) fig1<- ggarrange(plotlist=mylist, common.legend = TRUE, legend="topright", labels = c("(A)", "(B)"), font.label = list(size = 18, color = "black"), ncol=2) fig1 Kind regards Sibylle From: John Kane <jrkrideau at gmail.com> Sent: Monday, February 5, 2024 1:59 PM To: sibylle.stoeckli at gmx.ch Cc: r-help at r-project.org Subject: Re: [R] ggarrange & legend Could you supply us with a MWE (minimal working example)of what you have so far? Thanks. On Mon, 5 F...

...testing? https://www.rdocumentation.org/packages/base/versions/3.6.2/topics/all.equal I am actually interested in the x values (x-y coordinates). Test if x-y coordinates of both 25-m-pixel rasters are the same. Ther may be a small shift or differences in the number of decimal places. Kind regards Sibylle -----Original Message----- From: Ivan Krylov <ikrylov at disroot.org> Sent: Friday, August 16, 2024 11:45 AM To: sibylle.stoeckli at gmx.ch Cc: 'SIBYLLE ST?CKLI via R-help' <r-help at r-project.org> Subject: Re: [R] allequal diff ? Fri, 16 Aug 2024 11:32:58 +0200 <sibyll...

Hi Sibylle, For that kind of data with two different scales, I generally use two graphs that I name gg1 and gg2 and join them using gridExtra::grid.arrange(gg1, gg2). This way, the red part of your graph is easier to interpret. Have a nice day, Charles-?douard -----Message d'origine----- De?: R-help...

...otlist=mylist, common.legend = TRUE, legend="top", labels = c("(A)", "(B)"), font.label = list(size = 18, color = "black"), ncol=2) fig1 #===================================================================================== On Mon, 5 Feb 2024 at 08:44, <sibylle.stoeckli at gmx.ch> wrote: > Dear John Kane > > Dear R community > > > > Here my working example > > 1. Example that is working with legend=?top?. However, as mentioned, > the legend is in the middle of the top axis. > > mylist<-list(p1, p2) > &gt...

-------- Original-Nachricht -------- Datum: Wed, 27 Oct 2010 15:03:48 +0200 Von: "Sibylle St?ckli" <sibylle.stoeckli at gmx.ch> An: Rhelp <r-help at r-project.org> Betreff: as.list Dear R-users sorry, here Rcode included (attachment has been removed) I would like to read a txt file as list, to select rows and columns, and to create barplots. (1) selection of paramet...

...female trap catches barplot(m$trapcatch, ylab="Y", space=0.5, col=c("grey0"), ylim=c(0,450), las=2, cex.lab=0.9, cex.axis=0.9, cex.names=0.9) barplot(w$trapcatch,space=0.5, add=TRUE, beside=FALSE, col=c("grey50"), xaxt="n", yaxt="n") Thanks a lot Sibylle [[alternative HTML version deleted]]

...rm of multiple dispatch) and spend some time with it. There are many good ones out there. This S3 system is widely used within R and many packages, so doing this sort of homework now should serve you well in your future R journey. All IMO obviously. Cheers, Bert On Sun, Aug 18, 2024 at 11:00?AM SIBYLLE ST?CKLI via R-help <r-help at r-project.org> wrote: > > Dear Ivan > > Thanks a lot for this very nice example. > > Is it true that all.equal just compares y values? > Based on this help here I think so and the value I got is the difference for the y-values. > https://w...

...= W[[1]], time.scale = year, precip = Precip,...) where precip (small "p") is a formal argument of the dd() function and Precip is a column in the data frame W[[1]]. If this is wrong, my apologies, and feel free to ignore without responding. Best, Bert On Wed, Nov 6, 2024 at 10:59?AM Sibylle St?ckli <sibylle.stoeckli at gmx.ch> wrote: > Dear Rui > Dear Bert > > Many thanks > > Solution > filelist <- list.files(path = > "O:/Data-Work/2.../Daten_RA-MeteoCH_1990-2021", pattern='*.dat', all.files= > T, full.names= T) > AAR<-read....

? Thu, 18 Apr 2024 11:08:33 +0200 SIBYLLE ST?CKLI via R-help <r-help at r-project.org> ?????: > > #to check the index numbers of all imported raster list elements > > allrasters > list() > > > > #call single raster element > > allrasters[[1]] > Error in allrasters[[1]] : subscript out of boun...

...ar Gerrit Thanks a lot. "rbind" seems to be the right function. Unfortunately there is a shift in the x-axis (see pdf). There are 52 trapcatch values each, m and w, but m$trapcatch and w$trapcatch are shifted up to x-value 60. The follow-up lines for temp and humidity are fine. Thanks Sibylle setwd("~/Desktop/DatenLogger2017") # am Mac sks trap = read.delim("SWD_Trap_week-new.txt", na.strings="*", header=TRUE) climate = read.delim("Agrometeo_week-2017.txt", na.strings="*", header=TRUE) names(trap) names(climate) Hecke<...

Could you supply us with a MWE (minimal working example)of what you have so far? Thanks. On Mon, 5 Feb 2024 at 05:00, SIBYLLE ST?CKLI via R-help < r-help at r-project.org> wrote: > Dear R community > > It is possible to adjust the legend in combined ggplots using ggarrange > with > be positions top, bottom, left and right. > My question: Is there a function to change the position of the legend to...

...atrices. Adj. R square for the unique fraction [a] is 0.25. Does anyone know why the decomposition by hand using rda gives me a different result for [a] (constrained proportion is 0.32)? I used cbind() for the conditional fractions, but it should be similar to condition()? Thanks very much Sibylle > AZ_var<-varpart(AZ, OeAF, Acker, Gruen) > plot(AZ_var) > > AZ_var Partition of variation in RDA Call: varpart(Y = AZ, X = OeAF, Acker, Gruen) Explanatory tables: X1: OeAF X2: Acker X3: Gruen No. of explanatory tables: 3 Total variation (SS): 239101 Variance...

.../sd to complete a histogram (skewed). xs:range of y-values, ys: dnorm function Problem: I expected to multiply the ys function with the sample size (n=250-300). I was wondering about a factor between 12'000 and 30'000 to match the size of the dnorm line with the specific histogram. Thanks Sibylle hist(Biotree[Ld,]$Height2008, main="Larix decidua", ylim=c(0,50), xlab="Tree Height 2008 (cm)",col="aquamarine", font.main=3, cex.axis=0.8) xs<-0:650 ys<-dnorm(xs, mean=397.8, sd=97.6) lines(xs,ys*12000) -- GMX DSL SOMMER-SPECIAL: Surf & Phone Flat 16.000...

Hi, Sibylle, since you write '"mathematically" add', does barplot(rbind(m$trapcatch, w$trapcatch)) do what you want (modulo layout details)? Hth -- Gerrit --------------------------------------------------------------------- Dr. Gerrit Eichner Mathematical Institute...