Re: [R] apply

2024-10-04 Thread Jeff Newmiller via R-help
ctor: >>> >>> se<-apply(dd,1,sd) >>> >>> How can I calculate the covariance matrix using apply? Thanks. >>> >>> __ >>> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >

Re: [R] apply

2024-10-04 Thread Rui Barradas via R-help
273267965 0.0151896601 > On 10/4/2024 4:57 PM, Uwe Ligges wrote: Homework questions are not answered on this list. Best, Uwe Ligges On 04.10.2024 10:32, Steven Yen wrote: The following line calculates standard deviations of a column vector: se<-apply(dd,1,sd) How can I calculat

Re: [R] apply

2024-10-04 Thread Ivan Krylov via R-help
iance matrix for them. -- Best regards, Ivan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] apply

2024-10-04 Thread Ivan Krylov via R-help
one number. -- Best regards, Ivan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide comme

Re: [R] How to install this package

2024-10-03 Thread David Winsemius via R-help
- > *Roslinazairimah Zakaria* > *Tel: +609-5492370; Fax. No.+609-5492766* > > *Email: roslinazairi...@ump.edu.my ; > roslina...@gmail.com * > Faculty of Industrial Sciences & Technology > University Malaysia Pahang > Lebuhraya Tun Razak, 26300 Gambang, Pahang, Malaysia >

[R] (no subject)

2024-10-01 Thread Izmirlian, Grant (NIH/NCI) [E] via R-help
into submission, sort of). My overlooking unsplit() is somewhat impressive in view of "svn diff -c 18591" -pd > On 27 Sep 2024, at 17:08 , Martin Maechler wrote: > >>>>>> Chris Evans via R-help >>>>>>on Fri, 27 Sep 2024 12:20:47 +0200 w

Re: [R] Problem with converting grib file to excel

2024-09-30 Thread Roy Mendelssohn - NOAA Federal via R-help
df$d2m) > > # Convert valid_time from numeric to POSIXct assuming it's in seconds since > the epoch > df$valid_time <- as.POSIXct(df$valid_time, origin = "1970-01-01") > > # Extract year, month, day, and hour from valid_time > df$Year <- year(df$valid_time) > df$Month

Re: [R] when installing packages for R on Linux, is it better to use my distro's package manager, or install.packages()?

2024-09-29 Thread Ivan Krylov via R-help
no longer available. -- Best regards, Ivan [1] https://eddelbuettel.github.io/r2u/ __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] when installing packages for R on Linux, is it better to use my distro's package manager, or install.packages()?

2024-09-28 Thread Jeff Newmiller via R-help
install.packages("some_new_package")? > >I've done both over the years (perhaps a mistake that I will some day >regret.) > >Pros and cons of these two methods? > >Thanks. > >--Chris Ryan > -- Sent from my phone. Please excuse my brevity. ___

Re: [R] when installing packages for R on Linux, is it better to use my distro's package manager, or install.packages()?

2024-09-28 Thread Stephen H. Dawson, DSL via R-help
both over the years (perhaps a mistake that I will some day regret.) Pros and cons of these two methods? Thanks. --Chris Ryan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Is there a sexy way ...? Fortune nomination

2024-09-28 Thread Achim Zeileis via R-help
__ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code

Re: [R] Is there a sexy way ...?

2024-09-27 Thread Chris Evans via R-help
uador & Honorary Professor, University of Roehampton, London, UK. CORE site: http://www.coresystemtrust.org.uk Other work web site: https://www.psyctc.org/psyctc/ Personal site: https://www.psyctc.org/pelerinage2016/ ______ R-help@r-project.org mailing list -

Re: [R] Is there a sexy way ...?

2024-09-27 Thread Stephen Berman via R-help
of v, >> > >>> corresponding to level l of f are the entries of x[[l]]. I.e. I want >> > >>> v to equal >> > >>> >> > >>> c(7, 2, 6, 13, 5, 9, 1, 14, 15, 4, 8, 12, 10, 11, 3) >> > >>> >> > >>> I can create v

Re: [R] Problem with converting grib file to excel

2024-09-26 Thread Roy Mendelssohn - NOAA Federal via R-help
searched to read > > this type of file in R, I found the Raster Package. > > In python it is possible to do this using cdsapi and xarray library, but I > > am not familiar with python. > > Sincerely > > > > On Thu, Sep 26, 2024 at 2:33 AM Roy Mendelssohn - NOAA

Re: [R] Problem with converting grib file to excel

2024-09-26 Thread Roy Mendelssohn - NOAA Federal via R-help
hon it is possible to do this using cdsapi and xarray library, but I am > not familiar with python. > Sincerely > > On Thu, Sep 26, 2024 at 2:33 AM Roy Mendelssohn - NOAA Federal via R-help > mailto:r-help@r-project.org>> wrote: >> At least for me the dataset file did no

Re: [R] [Tagged] Re: Question about Date Object and time zones

2024-09-25 Thread Jeff Newmiller via R-help
... or a timezone string. ... or do what I usually do and use POSIXct for dates as well as times, self-managing use of midnight local time for "dates". On September 25, 2024 5:40:37 AM EDT, Ivan Krylov via R-help wrote: >24 сентября 2024 г. 17:10:13 GMT+03:00, Luca Brinkma

[R] plotting sf with lattice

2024-09-25 Thread Waichler, Scott R via R-help
on deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Problem with converting grib file to excel

2024-09-25 Thread Roy Mendelssohn - NOAA Federal via R-help
move rows with NA values >>> df <- na.omit(df) >>> >>> # Convert temperatures from Kelvin to Celsius >>> df$t2m <- df$t2m - 273.15 >>> df$d2m <- df$d2m - 273.15 >>> >>> # Calculate relative humidity >>> calculate_relative

Re: [R] How to install this package

2024-09-25 Thread Ivan Krylov via R-help
t work because 'termstrc' requires 'rgl' to load and your 'rgl' installation doesn't work. Try following the instructions in the README of the 'rgl' package regarding the OpenGL support. Perhaps it needs to be reinstalled after some updates you

Re: [R] Question about Date Object and time zones

2024-09-25 Thread Ivan Krylov via R-help
24 сентября 2024 г. 17:10:13 GMT+03:00, Luca Brinkmann via R-help пишет: > My current > understanding is, that a Date object does only save the days from the > origin and no more information about timezones or other information > (please correct me if I am wrong). You are correct.

[R] Question about Date Object and time zones

2024-09-24 Thread Luca Brinkmann via R-help
TE=German_Germany.utf8 LC_CTYPE=German_Germany.utf8    LC_MONETARY=German_Germany.utf8 [4] LC_NUMERIC=C    LC_TIME=German_Germany.utf8 time zone: Europe/Berlin tzcode source: internal __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more

Re: [R] store list objects in data.table

2024-09-22 Thread Ivan Krylov via R-help
Try wrapping the lm() expression in an extra list(): carsreg <- carsdt[, .(fit = list(lm(mpg ~ disp + hp + wt))), by = .(cyl)] # [1] 3 carsreg # cyl fit # # 1: 6 # 2: 4 # 3: 8 -- Best regards, Ivan ______ R-help@r-proje

Re: [R] "And" condition spanning over multiple columns in data frame

2024-09-12 Thread Ivan Krylov via R-help
# ...and finally attach to the original data.frame cbind(c10Dt) -- Best regards, Ivan __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Stop or limit to console printing large list and so on

2024-09-11 Thread Fer Arce via R-help
happening in the console, the same if you send a loong loop and want to abort...) cheers F. On 9/11/24 15:44, stephen sefick wrote: Stephen Sefick __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman

Re: [R] non-conformable arrays

2024-09-11 Thread Ivan Krylov via R-help
ot; "56" "59" "114" ... ..$ : chr "Year" - attr(*, "assign")= int 1 Browse[1]> t(random[['Year']]) |> str() num [1:19, 1:10] 0.175 0.181 0.102 0.119 0.158 ... ...and they are indeed non-conformable. Why is random$Year a matrix? -

Re: [R] Stop or limit to console printing large list and so on

2024-09-11 Thread Ivan Krylov via R-help
p in your case. -- Best regards, Ivan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented,

Re: [R] Reading a txt file from internet

2024-09-07 Thread Jeff Newmiller via R-help
7;t one to assume that it is big-endian data. But in this case there is a BOM so your machine has a buggy decoder? On September 7, 2024 2:43:24 PM PDT, Duncan Murdoch wrote: >On 2024-09-07 4:52 p.m., Jeff Newmiller via R-help wrote: >> When you specify LE in the encoding type, you are lo

Re: [R] Reading a txt file from internet

2024-09-07 Thread Jeff Newmiller via R-help
to contain embedded nulls >> >> Is there any way to read this data directly onto R? >> >> Thanks for your time >> > >The looks like a byte-order mark >(https://en.wikipedia.org/wiki/Byte_order_mark). >Try this: > >fn <- > file('https:

Re: [R] Reading a txt file from internet

2024-09-07 Thread Jeff Newmiller via R-help
= header, sep = sep, quote = quote, : > > line 4 appears to contain embedded nulls > >5: In read.table(file = file, header = header, sep = sep, quote = quote, : > > line 5 appears to contain embedded nulls > >Is there any way to read this data directly onto R? > >Thank

Re: [R] BUG: atan(1i) / 5 = NaN+Infi ?

2024-09-06 Thread Jorgen Harmse via R-help
work the answers for atan(1i) & atan(-1i) make sense. Regards, Jorgen Harmse. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read th

Re: [R] Fwd: effects() extractor for a quantile reqression object: error message

2024-09-06 Thread Gregg Powell via R-help
ive any typos please... I'm on a mobile device. Kind regards, Gregg Powell Sent from Proton Mail Android Original Message ---- On 06/09/2024 01:37, Koenker, Roger W wrote: > Apologies, forgot to copy R-help on this response. > > Begin forwarded message: > &

[R] effects() extractor for a quantile reqression object: error message

2024-09-05 Thread Christopher W. Ryan via R-help
used. Appreciate any suggestions. Thanks. --Chris Ryan __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] BUG: atan(1i) / 5 = NaN+Infi ?

2024-09-05 Thread Jeff Newmiller via R-help
vision/multiplication on a >> complex number. Which R actually does for these very particular cases; but >> not when only Im(x) is Inf. >> >> Sincerely, >> >> Leonard >> >> -- >> *From:* Bert Gunter >> *Sent:* Friday,

Re: [R] BUG: atan(1i) / 5 = NaN+Infi ?

2024-09-05 Thread Leo Mada via R-help
number. Which R actually does for these very particular cases; but not when only Im(x) is Inf. Sincerely, Leonard From: Bert Gunter Sent: Friday, September 6, 2024 1:12 AM To: Duncan Murdoch Cc: Leo Mada ; r-help@r-project.org Subject: Re: [R] BUG: atan(1i) / 5

[R] BUG: atan(1i) / 5 = NaN+Infi ?

2024-09-05 Thread Leo Mada via R-help
leted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Adding parameters for Benchmark normal distribution in shapiro.test

2024-09-02 Thread Jeff Newmiller via R-help
gt;For example in the code below, I want to test my sample to N(0, 10). > >shapiro.test(rnorm(100, mean = 5, sd = 3)) > >Is there any way to pass the information of the benchmark normal distribution? > >__ >R-help@r-project.org mail

Re: [R] Prediction from Arima model

2024-08-31 Thread David Winsemius via R-help
in the confidence interval of the new value. Is there any > method available to obtain that? > > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read th

Re: [R] aggregating data with quality control

2024-08-31 Thread Ivan Krylov via R-help
'C', 'D')) Since the factor is ordered and can be subject to comparison like status[1] < status[2], you can now use max() on your groups. If the sample contains any 'D's, max() will return a 'D', because it's larger than any 'C's. If th

[R] lattice panel layout like cross-tabs, like a 2 x 2 table.

2024-08-30 Thread Christopher W. Ryan via R-help
like strip.left() but for just one of the factors. I can't remember or find how I did it, what options to set. Can anyone remind me? Thanks. --Chris Ryan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/ma

Re: [R] Fill NA values in columns with values of another column

2024-08-29 Thread Jeff Newmiller via R-help
, 4L, 2L, 2L, NA, >> NA, NA, 5L, 9L, NA, NA, 10L, 7L, 2L, NA, 7L, NA, NA, NA, NA, >> 2L, 4L, 6L, 10L, NA, NA, NA), X2 = c(8L, 5L, 1L, 6L, 7L, 2L, >> 4L, 7L, 7L, 3L, 2L, 4L, 6L, 5L, 5L, 6L, 3L, 2L, 1L, 7L, 2L, 8L, >> 4L, 5L, 6L, 1L, 4L, 8L, 3L, 3L, 8L, 1L)), class = "

Re: [R] boxplot of raster and shapefile

2024-08-28 Thread SIBYLLE STÖCKLI via R-help
.. .. ..@ projargs: chr "+proj=somerc +lat_0=46.952405556 +lon_0=7.439583 +k_0=1 +x_0=260 +y_0=120 +ellps=bessel +units=m +no_defs" .. .. ..$ comment: chr "PROJCRS[\"unknown\",\nBASEGEOGCRS[\"unknown\",\n DATUM[\"Unknown based

[R] Fill NA values in columns with values of another column

2024-08-27 Thread Francesca PANCOTTO via R-help
ome>* -- [[alternative HTML version deleted]] ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/posting-guide.ht

Re: [R] boxplot of raster and shapefile

2024-08-27 Thread Ivan Krylov via R-help
В Mon, 26 Aug 2024 14:33:02 +0200 SIBYLLE STÖCKLI via R-help пишет: > > # Extract raster values within the shapefile > > extracted_values <- extract(raster_file, shape_file) > > # Assuming the shapefile has multiple polygons and you want to > > # create a box

[R] boxplot of raster and shapefile

2024-08-26 Thread SIBYLLE STÖCKLI via R-help
[[1]] > names(raster_file) [1] "Andrena.barbilabris_glo_ensemble" > > # Extract raster values within the shapefile > extracted_values <- extract(raster_file, shape_file) > > > > # Assuming the shapefile has multiple polygons and you want to create a boxplot

Re: [R] checking for unstated dependencies in examples ... ERROR

2024-08-25 Thread Ivan Krylov via R-help
В Sun, 25 Aug 2024 20:32:07 + Søren Højsgaard via R-help пишет: > checking for unstated dependencies in examples ... ERROR > Warning: parse error in file 'lines': > 77:20) Looks like one of your Rd \examples{} has a syntax error in it. If there is a file gRain.Rcheck

Re: [R] checking for unstated dependencies in examples ... ERROR

2024-08-25 Thread DAVID WINSEMIUS via R-help
> On 08/25/2024 1:32 PM PDT Søren Højsgaard via R-help > wrote: > > > Dear all, > When checking a package (ubuntu) I get > > checking for unstated dependencies in examples ... ERROR > Warning: parse error in file 'lines': > 77:20) > > Can

[R] checking for unstated dependencies in examples ... ERROR

2024-08-25 Thread Søren Højsgaard via R-help
memoise_2.0.1 broom_1.0.6 [76] httpuv_1.6.15 xfun_0.47 fs_1.6.4 [79] pkgconfig_2.0.3 ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do rea

Re: [R] paired raster boxplots

2024-08-25 Thread SIBYLLE STÖCKLI via R-help
rs_pres[[1]] > names(r) [1] "Andrena.barbilabris_glo_ensemble" > > > rs <- stack(r, s) > names(rs) <- c('r', 's') Error in `names<-`(`*tmp*`, value = c("r", "s")) : incorrect number of layer names > names(rs) [1] "Andrena.barbilabris

Re: [R] A question on Statistics regarding regression

2024-08-24 Thread Jeff Newmiller via R-help
ed. However I >just wanted to get some insight on above problem related to Statistics >from some of the great Statisticians in this forum. > >Thanks for your time. > >__ >R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see

Re: [R] paired raster boxplots

2024-08-24 Thread Ivan Krylov via R-help
s' is a list of rasters, 'r' could also contain more than one layer, resulting in stack(r, s) containing more than two layers. In order to avoid the error, you need to see names(rs) and either give the same number of names to the object instead of two, or additionally extract one

Re: [R] paired raster boxplots

2024-08-24 Thread SIBYLLE STÖCKLI via R-help
quot;DEBioBedeu" "FRBioBedeu" "ITBioBedeu" [15] "geometry" > s<-sf$Unterregio > r<-allrasters_pres[[1]] Kind regards Sibylle -Original Message- From: Ivan Krylov Sent: Friday, August 23, 2024 5:30 PM To: sibylle.stoec...@gmx.ch Cc: 'SIB

Re: [R] paired raster boxplots

2024-08-23 Thread Ivan Krylov via R-help
) : > 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] offer

Re: [R] paired raster boxplots

2024-08-23 Thread SIBYLLE STÖCKLI via R-help
ITBioBedeu" [15] "geometry" > 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

Re: [R] Linear regression and stand deviation at the Linux command line

2024-08-23 Thread Ivan Krylov via R-help
t from which you can extract the values you are interested in (see ?summary.lm). When in doubt, call help(name_of_the_function). -- Best regards, Ivan [1] https://cran.r-project.org/doc/manuals/R-intro.html [2] http://web.archive.org/web/20230106210646/http://ashipunov.info/shipunov/sch

Re: [R] fcaR not for latest R version

2024-08-22 Thread Ivan Krylov via R-help
n one of its dependencies, which I didn't look at. -- Best regards, Ivan __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide https://www.R-project.org/p

Re: [R] paired raster boxplots

2024-08-22 Thread Ivan Krylov via R-help
В Thu, 22 Aug 2024 08:46:03 +0200 SIBYLLE STÖCKLI via R-help пишет: > rr2s <- stack(r, r2,s) > > names(rs) <- c('r', 's', 'r2') > > Error in `names<-`(`*tmp*`, value = c("r", "s", "r2")) : > > incor

Re: [R] Force conversion of (POSIXct) time zone with base R

2024-08-22 Thread Ivan Krylov via R-help
27;) set: # time zone information successfully lost and replaced by UTC Sys.time() |> format() |> as.POSIXct(tz = 'UTC') -- Best regards, Ivan __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/m

[R] paired raster boxplots

2024-08-21 Thread SIBYLLE STÖCKLI via R-help
rr2s <- stack(r, r2,s) names(rs) <- c('r', 's', 'r2') d <- as.data.frame(rr2s) boxplot(r~s, data= d) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE an

Re: [R] allequal diff

2024-08-19 Thread SIBYLLE STÖCKLI via R-help
r2_resampled <- resample(r2, r1) compareRaster(r1, r2_resampled) CompareRaster: Evaluate whether a two or more Raster* objects have the same extent, number of rows and columns, projection, resolution, and origin (or a subset of these comparisons). Kind regards Sibylle -Original Message

Re: [R] allequal diff

2024-08-18 Thread SIBYLLE STÖCKLI via R-help
Message- From: Ivan Krylov Sent: Friday, August 16, 2024 11:45 AM To: sibylle.stoec...@gmx.ch Cc: 'SIBYLLE STÖCKLI via R-help' Subject: Re: [R] allequal diff В Fri, 16 Aug 2024 11:32:58 +0200 пишет: > # values and mask r1 > r1 <- getValues(r1) > mask1 <- is.na(r

Re: [R] Terminating a cmd windows from R

2024-08-17 Thread SIMON Nicolas via R-help
Many thanks for your help. I will consider the processx package. Best regards Nicolas Le 17 août 2024 à 16:07, Ivan Krylov a écrit :  В Sat, 17 Aug 2024 11: 47: 30 + SIMON Nicolas via R-help пишет: > nmrun <- paste0("cmd. exe /k ",nm_log,". bat ",nmi,".  c

Re: [R] Terminating a cmd windows from R

2024-08-17 Thread Ivan Krylov via R-help
В Sat, 17 Aug 2024 11:47:30 + SIMON Nicolas via R-help пишет: > nmrun <- paste0("cmd.exe /k ",nm_log,".bat ",nmi,".ctl ",nmi,".lst") You are using the /k option that instructs cmd.exe to keep the command prompt open. Does the batch file con

Re: [R] Terminating a cmd windows from R

2024-08-17 Thread SIMON Nicolas via R-help
t;.ctl ",nmi,".lst") # run system(nmrun, invisible = F, show.output.on.console = T, wait = T) Le 17 août 2024 à 13:31, Duncan Murdoch a écrit :  On 2024-08-17 6: 21 a. m. , SIMON Nicolas via R-help wrote: > I would like to stop a dos shell windows following the cmd (execute) command. Is there a

[R] Terminating a cmd windows from R

2024-08-17 Thread SIMON Nicolas via R-help
I would like to stop a dos shell windows following the cmd (execute) command. Is there a way to do that from R? __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] boxplot notch

2024-08-16 Thread Chris Evans via R-help
what you are trying to do in the geom_boxplot() call and how it uses the summarised data tibble. It may even be that if you do that, you will find what's causing the problem!  (I speak from bitter experience!!) Very best (all), Chris On 16/08/2024 17:51, SIBYLLE STÖCKLI via R-help wrote:

Re: [R] boxplot notch

2024-08-16 Thread SIBYLLE STÖCKLI via R-help
("red","darkgreen"), labels=c("ÖLN", "BIO"))+ scale_fill_manual(values=c("red","darkgreen"), labels= c("ÖLN", "BIO"))+ theme(legend.title = element_blank())+ theme(legend.text=element_text(size=20)) p1<-p

[R] boxplot notch

2024-08-16 Thread SIBYLLE STÖCKLI via R-help
olor_manual(values=c("red","darkgreen"), labels=c("�LN", "BIO"))+ scale_fill_manual(values=c("red","darkgreen"), labels= c("�LN", "BIO"))+ theme(legend.title = element_blank())+ theme(legend.text=element_text(size

Re: [R] allequal diff

2024-08-16 Thread Ivan Krylov via R-help
mask1 | mask2 # drops elements 1, 3, 4, 8 # [1] TRUE FALSE TRUE TRUE FALSE FALSE FALSE TRUE FALSE FALSE -- Best regards, Ivan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-he

Re: [R] allequal diff

2024-08-16 Thread SIBYLLE STÖCKLI via R-help
ask1 & mask2)]) [1] "'is.NA' value mismatch: 389 in current 56989152 in target" --> so there is just a mismatch in NA not in the xy pixels, right? Sibylle -Original Message- From: Ivan Krylov Sent: Friday, August 16, 2024 10:51 AM To: sibylle.stoec...@gm

Re: [R] allequal diff

2024-08-16 Thread SIBYLLE STÖCKLI via R-help
Here my idea including the error: > m1=r1[r1="NA",] > m2=r2_resampled[r2_resampled="NA",] > > > all.equal(getValues(r1)[!m1], getValues(r2_resampled)[!m2]) [1] "Numeric: lengths (80706867, 65806339) differ" -Original Message- From

Re: [R] allequal diff

2024-08-16 Thread Ivan Krylov via R-help
(), using is.na(getValues(...)) should definitely work. -- Best regards, Ivan ______ R-help@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/po

Re: [R] allequal diff

2024-08-16 Thread SIBYLLE STÖCKLI via R-help
LLE STÖCKLI via R-help Cc: sibylle.stoec...@gmx.ch Subject: Re: [R] allequal diff В Fri, 16 Aug 2024 07:19:38 +0200 SIBYLLE STÖCKLI via R-help пишет: > Is it possible to consider na.rm=TRUE? > > all.equal(getValues(r1), getValues(r2_resampled), tolerance = 0) > > [1] "

Re: [R] allequal diff

2024-08-16 Thread Ivan Krylov via R-help
В Fri, 16 Aug 2024 07:19:38 +0200 SIBYLLE STÖCKLI via R-help пишет: > Is it possible to consider na.rm=TRUE? > > all.equal(getValues(r1), getValues(r2_resampled), tolerance = 0) > > [1] "'is.NA' value mismatch: 9544032 in current 66532795 in target"

Re: [R] allequal diff

2024-08-15 Thread SIBYLLE STÖCKLI via R-help
x27;is.NA' value mismatch: 9544032 in current 66532795 in target" Kind regards Sibylle -Original Message- From: R-help On Behalf Of Ben Bolker Sent: Friday, August 16, 2024 1:06 AM To: r-help@r-project.org Subject: Re: [R] allequal diff Digging into

Re: [R] [EXTERNAL] Re: Very strange behavior of 'rep'

2024-08-15 Thread Izmirlian, Grant (NIH/NCI) [E] via R-help
Its where n.per.grp is first calculated. I rounded. Gosh do I feel stupid. Thanks to all who weighed in. Best Regards, Grant Izmirlian From: Duncan Murdoch Sent: Thursday, August 15, 2024 2:59 PM To: Izmirlian, Grant (NIH/NCI) [E] ; r-help@r-project.org

[R] allequal diff

2024-08-15 Thread SIBYLLE STÖCKLI via R-help
me values > print(result) [1] FALSE > Kind regards Sibylle [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

[R] geom_boxplot nocht=TRUE

2024-08-15 Thread SIBYLLE STÖCKLI via R-help
;ÖLN", "BIO"))+ scale_fill_manual(values=c("red","darkgreen"), labels= c("ÖLN", "BIO"))+ theme(legend.title = element_blank())+ theme(legend.text=element_text(size=20)) p1<-p1 + expand_limits(y=c(0, 80)) p1 [[alternative HT

Re: [R] [EXTERNAL] Re: Very strange behavior of 'rep'

2024-08-15 Thread Izmirlian, Grant (NIH/NCI) [E] via R-help
actorial design. node <- 29 n.per.grp <- n.per.grp.lst[conds[node, 1]] trnd <- trnds[conds[node, 2]] h <- H[conds[node, 3]] delta.j <- delta + h*shapes[[trnd]] ## covariates: ID.i <- c(outer(rep(1:n.per.grp,each=n.tt),c(0,n.per.grp),FUN="+")) arm.i <- rep(0:1, ea

[R] Very strange behavior of 'rep'

2024-08-15 Thread Izmirlian, Grant (NIH/NCI) [E] via R-help
ts, package:methods, Autoloads, package:base [[alternative HTML version deleted]] ______ R-help@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

Re: [R] geom_smooth with sd

2024-08-11 Thread SIBYLLE STÖCKLI via R-help
t: Sunday, August 11, 2024 6:30 PM To: sibylle.stoec...@gmx.ch Cc: R-help@r-project.org Subject: Re: [R] geom_smooth with sd Hi! This is probably completely off base, but your ymin and y max setup lines are different. One uses sqrt(y), while the second uses sqrt(length(y)). Could that p

[R] geom_smooth with sd

2024-08-11 Thread SIBYLLE STÖCKLI via R-help
Dear community Using after_stat() I was able to visualise ggplot with standard deviations instead of a confidence interval as seen in the R help. p1<-ggplot(data = MS1, aes(x= Jahr, y= QI_A,color=Bio, linetype=Bio)) + geom_smooth(aes(fill=Bio, ymax=after_stat(y+se*s

Re: [R] Printing

2024-08-11 Thread Ivan Krylov via R-help
(printing) branch? Then you would be able to always return(invisible(list("ei"=ap,"vi"=vap))), whether printing is TRUE or not. -- Best regards, Ivan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.et

Re: [R] a fast way to do my job

2024-08-10 Thread Yuan Chun Ding via R-help
e.rpkm$purity2) then running result2 <- residuals(lm.fit( x= pur2.1, y = dat)); now I am thinking whether an intercept is required or not. Ding From: R-help On Behalf Of Yuan Chun Ding via R-help Sent: Saturday, August 10, 2024 12:30 PM To: Bert Gunter ; Ben Bolker Cc: r-help@r-project.org Subjec

Re: [R] a fast way to do my job

2024-08-10 Thread Yuan Chun Ding via R-help
You are right. I also just thought about that, no intercept is not applicable to my case. Ding From: Bert Gunter Sent: Saturday, August 10, 2024 1:06 PM To: Yuan Chun Ding Cc: Ben Bolker ; r-help@r-project.org Subject: Re: [R] a fast way to do my job Ah, messages crossed. A no-intercept

Re: [R] a fast way to do my job

2024-08-10 Thread Yuan Chun Ding via R-help
after add intercept, all residuals are the same from lm or lm.fit. Ding From: Bert Gunter Sent: Saturday, August 10, 2024 1:00 PM To: Yuan Chun Ding Cc: Ben Bolker ; r-help@r-project.org Subject: Re: [R] a fast way to do my job Probably because you inadvertently ran different models. Without

Re: [R] a fast way to do my job

2024-08-10 Thread Yuan Chun Ding via R-help
0.99. I have not found the reason. Thanks, Ding From: Bert Gunter Sent: Friday, August 9, 2024 7:11 PM To: Ben Bolker Cc: Yuan Chun Ding ; r-help@r-project.org Subject: Re: [R] a fast way to do my job Better idea, Ben! It would work as you might expect it to to produce the same results as the

Re: [R] If loop

2024-08-10 Thread martin gregory via R-help
ng the list. Martin ______ R-help@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.

Re: [R] a fast way to do my job

2024-08-09 Thread Yuan Chun Ding via R-help
nd inform the sender that you do not wish to receive further e-mail from the sender. (LCP301) [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.eth

Re: [R] WDI package commands timing out and not working

2024-08-09 Thread Ivan Krylov via R-help
ss discriminated against by Cloudflare. Maybe petition the WDI maintainer to export a function that parses a JSON file manually downloaded with a browser? -- Best regards, Ivan ______ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https:

Re: [R] WDI package commands timing out and not working

2024-08-09 Thread Ivan Krylov via R-help
st by accessing it through Tor. -- Best regards, Ivan __ R-help@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.ht

Re: [R] An error message with the command fm<-1m

2024-08-06 Thread Jorgen Harmse via R-help
oval must be 'o' or 'O', not the digit '0'. (Digits after the first character are standard.) Regards, Jorgen Harmse. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIB

Re: [R] Difference between stats.steps() and MuMIn.dredge() to select best fit model

2024-08-01 Thread Ivan Krylov via R-help
pick d before a. A greedy algorithm for K variables out N tests N + (N-1) + ... + (N-K+1) = N*K - K(K-1)/2 models. An exhaustive search would have to test choose(N,K) = N!/(N-K)!/K! models. -- Best regards, Ivan __ R-help@r-project.org mailing

Re: [R] round and trailing zero

2024-07-30 Thread Jorgen Harmse via R-help
latter case, why stash them in a variable? For more control of the display, consider sprintf (or a wrapper that combines sprintf with cat). Regards, Jorgen Harmse. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To

Re: [R] Using optim() function to find MLE

2024-07-28 Thread Ivan Krylov via R-help
optim() doesn't support the truncated (line search) Newton-CG method. See the 'optimx' and 'nloptr' packages for an implementation of a truncated Newton method (not necessarily exactly the same one). -- Best regards, Ivan ______

Re: [R] please help generate a square correlation matrix

2024-07-27 Thread Yuan Chun Ding via R-help
O'Keefe ; r-help@r-project.org Subject: Re: [R] please help generate a square correlation matrix Your expanded explanation helps clarify your intent. Herewith some comments. Of course, feel free to ignore and not respond. And, as always, my apologies if I have failed to comprehend your intent.

Re: [R] please help generate a square correlation matrix

2024-07-27 Thread Yuan Chun Ding via R-help
mutation, then it is informative to know how strong the negative correlation when considering different combination of treatment strategies. Ding From: R-help On Behalf Of Richard O'Keefe Sent: Saturday, July 27, 2024 4:47 AM To: Bert Gunter Cc: r-help@r-project.org Subject: Re: [R] please

Re: [R] plotting nnet function....

2024-07-27 Thread Ivan Krylov via R-help
; package: https://cran.r-project.org/package=caret/vignettes/caret.html If you do encounter an error on your way or get stuck not knowing how exactly to continue, please ask a more specific question. -- Best regards, Ivan ______ R-help@r-project.org ma

Re: [R] Automatic Knot selection in Piecewise linear splines

2024-07-26 Thread Vito Muggeo via R-help
(4, 2/3), adjustcolor("orange4", 2/3)) lwds <- c(3, 2, 2) matlines(xx, fits, col = cols, lwd = lwds, lty=1) legend("topleft", c("true f(x)", "smooth.spline()", "earth()"), col=cols, lwd=lwds, bty = "n") title(paste("earth

Re: [R] please help generate a square correlation matrix

2024-07-25 Thread Yuan Chun Ding via R-help
.test(tem2[,1],tem2[,2]) r[i, j] <- tmp3$estimate P[i, j] <- tmp3$p.value } } } r<-as.data.frame(r) P<-as.data.frame(P) From: R-help On Behalf Of Yuan Chun Ding via R-help Sent: Thursday, July 25, 2024 11:26 AM To: Rui Barradas ; r-help@r-project.org Subject: Re: [

Re: [R] please help generate a square correlation matrix

2024-07-25 Thread Yuan Chun Ding via R-help
HI Rui, Thank you for the help! You did not remove a row if zero values exist in both column pair, right? Ding From: Rui Barradas Sent: Thursday, July 25, 2024 11:15 AM To: Yuan Chun Ding ; r-help@r-project.org Subject: Re: [R] please help generate a square correlation matrix Às 17: 39 de

  1   2   3   4   5   6   7   8   9   10   >