Re: [R] Import multiple tif raster

2024-04-18 Thread Ivan Krylov via R-help
В Thu, 18 Apr 2024 11:08:33 +0200 SIBYLLE STÖCKLI via R-help пишет: > > #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 bounds `allrasters` is an

[R] Import multiple tif raster

2024-04-18 Thread SIBYLLE STÖCKLI via R-help
Dear community My aim is to import multiple .tif raster files using the help here: https://stackoverflow.com/questions/52746936/how-to-efficiently-import-multi ple-raster-tif-files-into-r Does anyone now about the error "subscripts out of bounds"? I am not sure about the check with alllrasters:

Re: [R] missing() fails to detect missing

2024-04-18 Thread Peter Dalgaard via R-help
I'm not really sure how verboten it is to have methods that omit "...". The code seems to be instrumented to handle that case. The catch is that it does so by inserting a call to .local where .local is the method function, e.g. if you insert a browser() call at the start of the method, you'll

Re: [R] missing() fails to detect missing

2024-04-17 Thread Jeff Newmiller via R-help
Fascinating. Or, well, not. Failing to use a method signature that is compatible with the generic is a no-no. So your bug seems to me to be outside the bounds of how R is supposed to be used. So don't do that. On April 17, 2024 4:25:38 PM PDT, "Boylan, Ross via R-help" wrote: >When a generic

[R] missing() fails to detect missing

2024-04-17 Thread Boylan, Ross via R-help
When a generic (S4) has an argument with a default followed by ..., missing() doesn't seem to work if the method omits the ... Sample--- foo <- function(x, y=0, ...){ "you are very generic" } # no ... in function arguments setMethod("foo",

[R-es] Ayuda: Calculo tamaño del efecto

2024-04-17 Thread Rossana Tornel Vázquez
Necesito realizar un an�lisis del tama�o de muestra para un experimento con las siguientes caracter�sticas: 1) Una variable independiente intrasujeto es el tipo de queso con dos niveles, Queso A (1) y Queso B (2). 2) La otra variable independiente intrasujeto es la ronda de degustaci�n con tres

Re: [R] Synthetic Control Method

2024-04-17 Thread Petr Pikal
Hallo Nadja Similar as Bert I do not know how the function works. From the help page - synth.data - is used in example. Check if your data structure is consistent with synth.data by comparing str(synth.data) and str( INVESTMENTVOLUME). Names of columns should be in both cases match the names in

Re: [R] Synthetic Control Method

2024-04-16 Thread Bert Gunter
Note that your unit.variable and unit.names.variable are indentical. Is this what you intended? (I have no idea how the Synth package works). Bert On Tue, Apr 16, 2024 at 12:58 AM wrote: > Good Morning > > > > I want to perform a synthetic control method with R. For this purpose, I > created

[R] Synthetic Control Method

2024-04-16 Thread nadja.delliehausen
Good Morning I want to perform a synthetic control method with R. For this purpose, I created the following code: # Re-load packages library(Synth) library(readxl) # Pfadeinstellung Excel-Blatt excel_file_path <-

[R] Reminder: Question regarding reservoir volume and water level

2024-04-15 Thread javad bayat
On Wed, Apr 10, 2024 at 11:05 PM javad bayat wrote: > Dear all; > Thank you for your reply. > David has explained an interesting method. > David I have DEM file of the region and I have extracted the xyz data from > that. > Also I can extract bathymetry data as xyz file. > I have calculated the

Re: [R] any and all

2024-04-13 Thread avi.e.gross
Yes, Lennart, I have been looking at doing something like you say by using the vectorized ways the tidyverse is now offering. For my application, if the naming was consistent, an approach like yours is good, albeit has to be typed carefully. When I cannot control the names but have to lump

Re: [R] Just for your (a|be)musement.

2024-04-13 Thread avi.e.gross
Richard, The code you show is correct and it does not include where you say ChatGTP explained it was 33/26 rather than the correct 42/216. I gather it have the proper fraction for the other two scenarios. So what would cause such a localized error? The method chosen is to try all possible

Re: [R] any and all

2024-04-13 Thread Lennart Kasserra
Hi Avi, As Dénes Tóth has rightly diagnosed, you are building an "all or nothing" filter. However, you do not need to explicitly spell out all columns that you want to filter for; the "tidy" way would be to use a helper function like `if_all()` or `if_any()`. Consider this example (I hope I

[R] Just for your (a|be)musement.

2024-04-13 Thread Richard O'Keefe
I recently had the chance to read a book explaining how to use ChatGPT with a certain programming language. (I'm not going to describe the book any more than that because I don't want to embarrass whoever wrote it.) They have appendix material showing three queries to ChatGPT and the answers.

Re: [R] any and all

2024-04-12 Thread avi.e.gross
Thanks everyone and any/all reading this. I think I got my answer. And, no, I suspected I did not need to provide a very specific example, at least not yet. The answer is that my experiment was not vectorized while using dplyr verbs like mutate do their work implicitly in a vectorized way.

Re: [R] any and all

2024-04-12 Thread Dénes Tóth
Hi Avi, As Duncan already mentioned, a reproducible example would be helpful to assist you better. Having said that, I think you misunderstand how `dplyr::filter` works: it performs row-wise filtering, so the filtering expression shall return a logical vector of the same length as the

Re: [R] any and all

2024-04-12 Thread Duncan Murdoch
On 12/04/2024 3:52 p.m., avi.e.gr...@gmail.com wrote: Base R has generic functions called any() and all() that I am having trouble using. It works fine when I play with it in a base R context as in: all(any(TRUE, TRUE), any(TRUE, FALSE)) [1] TRUE all(any(TRUE, TRUE), any(FALSE, FALSE))

[R] any and all

2024-04-12 Thread avi.e.gross
Base R has generic functions called any() and all() that I am having trouble using. It works fine when I play with it in a base R context as in: > all(any(TRUE, TRUE), any(TRUE, FALSE)) [1] TRUE > all(any(TRUE, TRUE), any(FALSE, FALSE)) [1] FALSE But in a tidyverse/dplyr environment, it

Re: [R] Debugging functions defined (locally) inside another functions

2024-04-12 Thread Iago Giné Vázquez
Thanks a lot both Duncan and Ivan, I will keep that example in mind, Duncan, great! Best regards, Iago De: Duncan Murdoch Enviat el: divendres, 12 d�abril de 2024 15:36 Per a: Iago Gin� V�zquez ; r-help@r-project.org Tema: Re: [R] Debugging functions defined

Re: [R] Debugging functions defined (locally) inside another functions

2024-04-12 Thread Duncan Murdoch
On 12/04/2024 8:15 a.m., Iago Giné Vázquez wrote: Hi all, I am trying to debug an error of a function g defined and used inside another function f of a package. So I have f <- function(whatever){ ... g <- function(whatever2){ ... } ... } If I wanted to debug some thing

Re: [R] Debugging functions defined (locally) inside another functions

2024-04-12 Thread Ivan Krylov via R-help
В Fri, 12 Apr 2024 12:53:02 + Iago Giné Vázquez пишет: > How should I call trace() if f was a function? Let the tracer be quote(debug(g)) and use as.list(body(f)) to determine where it should be injected: f <- function() { message('exists("g") so far is ', exists('g')) g <- function() {

Re: [R] Debugging functions defined (locally) inside another functions

2024-04-12 Thread Iago Giné Vázquez
Thank you Ivan, your example solves my issue this time through debug(environment(Adder$add)$add) Just for the future, you say Moreover, `g` doesn't exist at all until f() is evaluated and reaches this point. If `f` was a function, it would be possible to trace() it, inserting a

Re: [R] Debugging functions defined (locally) inside another functions

2024-04-12 Thread Ivan Krylov via R-help
В Fri, 12 Apr 2024 12:15:07 + Iago Giné Vázquez пишет: > f <- function(whatever){ >... >g <- function(whatever2){ > ... >} >... > } > > If I wanted to debug some thing directly inside f I would do > debug(f). But this does not go inside g code. On the other hand, >

Re: [R] Debugging functions defined (locally) inside another functions

2024-04-12 Thread Iago Giné Vázquez
To be precise, in the case I am looking this time f is not a function, but f <- ggplot2::ggproto(...) So debug(f) produces Error in debug(f) : argument must be a function Iago De: R-help de part de Iago Gin� V�zquez Enviat el: divendres, 12 d�abril de 2024

[R] Debugging functions defined (locally) inside another functions

2024-04-12 Thread Iago Giné Vázquez
Hi all, I am trying to debug an error of a function g defined and used inside another function f of a package. So I have f <- function(whatever){ ... g <- function(whatever2){ ... } ... } If I wanted to debug some thing directly inside f I would do debug(f). But this does not

Re: [R] Regexp pattern but fixed replacement?

2024-04-11 Thread Duncan Murdoch
On 11/04/2024 12:58 p.m., Iris Simmons wrote: Hi Duncan, I only know about sub() and gsub(). There is no way to have pattern be a regular expression and replacement be a fixed string. Backslash is the only special character in replacement. If you need a reference, see this file:

Re: [R] Regexp pattern but fixed replacement?

2024-04-11 Thread Duncan Murdoch
On 11/04/2024 12:57 p.m., Dave Dixon wrote: Backslashes in regex expressions in R are maddening, but they make sense. R string handling interprets your replacement string "\\" as just one backslash. Your string is received by gsub as "\" - that is, just the control backslash, NOT the character

Re: [R] Regexp pattern but fixed replacement?

2024-04-11 Thread Iris Simmons
Hi Duncan, I only know about sub() and gsub(). There is no way to have pattern be a regular expression and replacement be a fixed string. Backslash is the only special character in replacement. If you need a reference, see this file:

Re: [R] Regexp pattern but fixed replacement?

2024-04-11 Thread Dave Dixon
Backslashes in regex expressions in R are maddening, but they make sense. R string handling interprets your replacement string "\\" as just one backslash. Your string is received by gsub as "\" - that is, just the control backslash, NOT the character backslash. gsub is expecting to see \0,

[R] Regexp pattern but fixed replacement?

2024-04-11 Thread Duncan Murdoch
I noticed this issue in stringr::str_replace, but it also affects sub() in base R. If the pattern in a call to one of these needs to be a regular expression, then backslashes in the replacement text are treated specially. For example, gsub("a|b", "\\", "abcdef") gives "def", not

Re: [R] seeking help with splicing in R

2024-04-11 Thread Ebert,Timothy Aaron
Whatever you had as HTML was deleted. If that was data we did not get it. 1) manipulate wpi_def_nic2004 and wpi_def_nic2008 first so that the data are compatible, then join them. 2) The full_join statement should explicitly state the columns to join by. Using by=NULL joins by all the columns

[ESS] Parameterized quarto qmd in ESS?

2024-04-11 Thread Micha Silver via ESS-help
I'm working on a quarto template with ESS that runs 'R' code chunks . When using RStudio, if there are parameters set in the yaml header, the default values are available in code chunks (for testing the template). However in ESS the `params$...` variables are not available. How would I create

Re: [ESS] Displaying R plots within an Emacs buffer

2024-04-11 Thread Stephen J. Eglen via ESS-help
I've tidied up the code for showing R plots as SVG in an Emacs buffer; it is now on melpa: https://melpa.org/#/essgd so package.el can grab it easily enough. Best wishes, Stephen On Wed, Mar 27 2024, Stephen J. Eglen wrote: I've got it now working locally so that the aspect ratio is

[R] seeking help with splicing in R

2024-04-11 Thread Manya Budhiraja
Hello! I was trying to splice two wholesale price index deflators series which have different base years. One of them is called wpi_def_nic2004(from 2005 to 2012), and another is called wpi_def_nic2008(from 2012 to 2019). I am trying to create a single series such that the base year prices are

[R] [R-pkgs] New package - fastTS

2024-04-11 Thread Peterson, Ryan
Hi R enthusiasts, I am happy to announce a new package available on CRAN: fastTS (https://cran.r-project.org/web/packages/fastTS/). fastTS is especially useful for large time series with exogenous features and/or complex seasonality (i.e. with multiple modes), allowing for possibly

Re: [R] Question regarding reservoir volume and water level

2024-04-10 Thread javad bayat
Dear all; Thank you for your reply. David has explained an interesting method. David I have DEM file of the region and I have extracted the xyz data from that. Also I can extract bathymetry data as xyz file. I have calculated the storage (volume) of reservoir at the current elevation. But the

Re: [R] Exceptional slowness with read.csv

2024-04-10 Thread avi.e.gross
Dave, Your method works for you and seems to be a one-time fix of a corrupted data file so please accept what I write not as a criticism but explaining my alternate reasoning which I suspect may work faster in some situations. Here is my understanding of what you are doing: You have a file in

Re: [R] Exceptional slowness with read.csv

2024-04-10 Thread Dave Dixon
That's basically what I did 1. Get text lines using readLines 2. use tryCatch to parse each line using read.csv(text=...) 3. in the catch, use gregexpr to find any quotes not adjacent to a comma (gregexpr("[^,]\"[^,]",...) 4. escape any quotes found by adding a second quote (using str_sub from

Re: [R] Exceptional slowness with read.csv

2024-04-10 Thread avi.e.gross
It sounds like the discussion is now on how to clean your data, with a twist. You want to clean it before you can properly read it in using standard methods. Some of those standard methods already do quite a bit as they parse the data such as looking ahead to determine the data type for a

Re: [R] Exceptional slowness with read.csv

2024-04-10 Thread Rui Barradas
Às 06:47 de 08/04/2024, Dave Dixon escreveu: Greetings, I have a csv file of 76 fields and about 4 million records. I know that some of the records have errors - unmatched quotes, specifically. Reading the file with readLines and parsing the lines with read.csv(text = ...) is really slow. I

Re: [R] Problem with base::order

2024-04-10 Thread Sigbert Klinke
Hi, you are unfortunately right. Executing x <- sample(c(1,2,NA), 26, replace=TRUE) y <- sample(c(1,2,NA), 26, replace=TRUE) o <- order(x, y, decreasing = c(T,F), na.last=c(F,T)) cbind(x[o], y[o]) shows that the second entry of na.last is ignored without warning. Thanks Sigbert Am 10.04.24

Re: [R] Problem with base::order

2024-04-10 Thread Ivan Krylov via R-help
В Wed, 10 Apr 2024 09:33:19 +0200 Sigbert Klinke пишет: > decreasing=c(F,F,F) This is only documented to work with method = 'radix': >> For the ‘"radix"’ method, this can be a vector of length equal to >> the number of arguments in ‘...’ and the elements are recycled as >> necessary. For the

[R] Problem with base::order

2024-04-10 Thread Sigbert Klinke
Hi, when I execute order(letters, LETTERS, 1:26) then everything is fine. But if I execute order(letters, LETTERS, 1:26, na.last=c(T,T,T), decreasing=c(F,F,F)) I get the error message Error in method != "radix" && !is.na(na.last) : 'length = 3' in constraint to 'logical(1)' Shouldn't both

Re: [R] cpgram: confidence band with a chirp signal

2024-04-09 Thread Laurent Rhelp
I answer to myself: I have to generate my signal up to f1 = Fs/2 in order to have a flat response over the entire frequency range. Le 09/04/2024 à 19:48, laurentRhelp a écrit : Dear RHelp-list,    I generate a swept sine signal using the signal library. I can see using the spec.pgram command

Re: [R] Building R-4,3,3 fails

2024-04-09 Thread Rich Shepard
On Tue, 9 Apr 2024, Ivan Krylov wrote: That's fine, R will run straight from the build directory. It has to do so in order to compile the vignettes. Ivan, That's good to know. Thanks. But let's skip this step. Here's reshape.tex from R-4.3.3: https://0x0.st/XidU.tex/reshape.tex (Feel free

Re: [R] Building R-4,3,3 fails

2024-04-09 Thread Rich Shepard
On Tue, 9 Apr 2024, Ivan Krylov wrote: At this point in the build, R already exists, is quite operable and even has all the recommended packages installed. The build system then uses this freshly compiled R to run Sweave on the vignettes. Let me break the build in a similar manner and see what

[R] cpgram: confidence band with a chirp signal

2024-04-09 Thread laurentRhelp
Dear RHelp-list,    I generate a swept sine signal using the signal library. I can see using the spec.pgram command that the spectrum of this signal is white. But when I am calculating the cumulative periodogram using the command cpgram the process doesn't stay inside the confidence band (cf.

Re: [R] Question regarding reservoir volume and water level

2024-04-09 Thread David Stevens via R-help
Water engineer here. The standard approach is to 1) get the storage vs. elevation data from the designers of the reservoir or, barring that, 2) get the bathymetry data from USBR or state DWR, or, if available, get the DEM data from USGS if the survey was done before the reservoir was built or

Re: [R] CEoptim problems

2024-04-09 Thread Adelchi Azzalini
Thanks for the suggestion, Ivan. The issue has been overcome with a simple change of the code to the form (is.null(A) || any(is.na(A)) following advice from Peter Dalgard. However, I have kept a note of the R Inferno reference, for future problems. Best wishes, Adelchi > On 9 Apr 2024, at

Re: [R] CEoptim problems

2024-04-09 Thread Adelchi Azzalini
> On 9 Apr 2024, at 14:54, peter dalgaard wrote: > > Hi, Adelchi, > > Depends on what you want help with... > > The proximate cause would seem to be that the code ought to have "is.null(A) > || any(is.NA(A))", which I presume you could fairly easily fix for yourself > in the package

Re: [R] Question regarding reservoir volume and water level

2024-04-09 Thread peter dalgaard
So, you know how to get volume for given water level. For the reverse problem, you get in trouble because of the nonlinearity inherent in the dependence of surface area on the level. I don't think there is a simple solution to this, save for mapping out the volume as a function of water

Re: [R] CEoptim problems

2024-04-09 Thread peter dalgaard
Hi, Adelchi, Depends on what you want help with... The proximate cause would seem to be that the code ought to have "is.null(A) || any(is.NA(A))", which I presume you could fairly easily fix for yourself in the package sources or even locally in an R session. Vector-valued logicals in flow

Re: [R] How to set the correct libomp for R

2024-04-09 Thread Ivan Krylov via R-help
В Tue, 9 Apr 2024 09:55:06 +0200 gernop...@gmx.net пишет: > If I only move away /usr/local/lib/libomp.dylib, I can still install > it. So it seems that also here the internal libomp.dylib from R is > used. Just the bundled omp files at /usr/local/include (omp-tools.h, > omp.h, ompt.h) seem to be

Re: [R] CEoptim problems

2024-04-09 Thread Ivan Krylov via R-help
В Tue, 9 Apr 2024 12:04:26 +0200 Adelchi Azzalini пишет: > res <- CEoptim(sumsqrs, f.arg = list(xt), continuous = list(mean = > c(0, 0, 0), sd = rep(1,3), conMat = A, conVec = b), discrete = > list(categories = c(298L, 298L), smoothProb = 0.5),N = 1, rho > = 0.001) > > Error in

[R] CEoptim problems

2024-04-09 Thread Adelchi Azzalini
In the attempt to explore the usage of package CEoptim, I have run the code listed at the end of this message. This code is nothing but the one associated to example 5.7 in the main reference of the package, available at https://www.jstatsoft.org/article/view/v076i08 and is included in the

Re: [R] How to set the correct libomp for R

2024-04-09 Thread gernophil--- via R-help
Sorry, if have to correct this. If I only move away /usr/local/lib/libomp.dylib, I can still install it. So it seems that also here the internal libomp.dylib from R is used. Just the bundled omp files at /usr/local/include (omp-tools.h, omp.h, ompt.h) seem to be used. So maybe this is caused

Re: [R] How to set the correct libomp for R

2024-04-09 Thread gernophil--- via R-help
Sorry fort he late reply, your mail ended up in my spam and I've just seen it recently. > Does the behaviour change if you temporarily move away > /usr/local/lib/libomp.dylib? It does not change the behavior after loading (or attaching) data.table using "library(data.table)". It is still

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread jim holtman
Try reading the lines in (readLines), count the number of both types of quotes in each line. Find out which are not even and investigate. On Mon, Apr 8, 2024, 15:24 Dave Dixon wrote: > I solved the mystery, but not the problem. The problem is that there's > an unclosed quote somewhere in those

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Eberhard W Lisse
I find QSV very helpful. el On 08/04/2024 22:21, Dave Dixon wrote: > I solved the mystery, but not the problem. The problem is that > there's an unclosed quote somewhere in those 5 additional records I'm > trying to access. So read.csv is reading million-character fields. > It's slow at that.

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Dave Dixon
Right, I meant to add header=FALSE. And, it looks now like the next line is the one with the unclosed quote, so read.csv is trying to read million-character headers! On 4/8/24 12:42, Ivan Krylov wrote: В Sun, 7 Apr 2024 23:47:52 -0600 Dave Dixon пишет: > second_records <-

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Dave Dixon
Good suggestion - I'll look into data.table. On 4/8/24 12:14, CALUM POLWART wrote: > data.table's fread is also fast. Not sure about error handling. But I > can merge 300 csvs with a total of 0.5m lines and 50 columns in a > couple of minutes versus a lifetime with read.csv or readr::read_csv >

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Dave Dixon
Thanks, yeah, I think scan is more promising. I'll check it out. On 4/8/24 11:49, Bert Gunter wrote: > No idea, but have you tried using ?scan to read those next 5 rows? It > might give you a better idea of the pathologies that are causing > problems. For example, an unmatched quote might

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Dave Dixon
I solved the mystery, but not the problem. The problem is that there's an unclosed quote somewhere in those 5 additional records I'm trying to access. So read.csv is reading million-character fields. It's slow at that. That mystery solved. However, the the problem persists: how to fix what is

Re: [R] Building R-4,3,3 fails

2024-04-08 Thread Rich Shepard
On Mon, 8 Apr 2024, Ivan Krylov wrote: A Web search suggests that texi2dvi may output this message by mistake when the TeX installation is subject to a different problem: https://web.archive.org/web/20191006123002/https://lists.gnu.org/r/bug-texinfo/2016-10/msg00036.html Ivan, That thread is

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Rui Barradas
Às 19:42 de 08/04/2024, Ivan Krylov via R-help escreveu: В Sun, 7 Apr 2024 23:47:52 -0600 Dave Dixon пишет: > second_records <- read.csv(file_name, skip = 2459465, nrows = 5) It may or may not be important that read.csv defaults to header = TRUE. Having skipped 2459465 lines, it may

Re: [R] Building R-4,3,3 fails

2024-04-08 Thread Rich Shepard
On Mon, 8 Apr 2024, Ivan Krylov wrote: Questions about building R do get asked here and R-devel. Since you're compiling a released version of R and we don't have an R-SIG-Slackware mailing list, R-help sounds like the right place. Ivan, Okay: What are the last lines of the build log,

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Ivan Krylov via R-help
В Sun, 7 Apr 2024 23:47:52 -0600 Dave Dixon пишет: > > second_records <- read.csv(file_name, skip = 2459465, nrows = 5) It may or may not be important that read.csv defaults to header = TRUE. Having skipped 2459465 lines, it may attempt to parse the next one as a header, so the second call

[R] Building R-4,3,3 fails

2024-04-08 Thread Rich Shepard
I've been building R versions for years with no issues. Now I'm trying to build R-4.3.3 on Slackware64-15.0 (fully patched) with TeXLive2024 (fully patched) installed. The error occurs building a vignette. Is this mail list the appropriate place to ask for help or should I post the request on

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread CALUM POLWART
data.table's fread is also fast. Not sure about error handling. But I can merge 300 csvs with a total of 0.5m lines and 50 columns in a couple of minutes versus a lifetime with read.csv or readr::read_csv On Mon, 8 Apr 2024, 16:19 Stevie Pederson, wrote: > Hi Dave, > > That's rather

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Bert Gunter
No idea, but have you tried using ?scan to read those next 5 rows? It might give you a better idea of the pathologies that are causing problems. For example, an unmatched quote might result in some huge number of characters trying to be read into a single element of a character variable. As your

Re: [R] duplicated() on zero-column data frames returns empty

2024-04-08 Thread Jorgen Harmse via R-help
I appreciate the compliment from Ivan and still share the puzzlement at the empty return. What is the policy for changing something that is wrong? There is a trade-off between breaking old code that worked around a problem and breaking new code written by people who make reasonable

[R] Questions about ks.test function {stats}

2024-04-08 Thread Jin, Ziyan
Dear R-help, Hope this email finds you well. My name is Ziyan. I am a graduate student in Zhejiang University. My subject research involves ks.test in stats-package {stats}. Based on the code, I have two main questions. Could you provide me some more information? I download different

Re: [R] Exceptional slowness with read.csv

2024-04-08 Thread Stevie Pederson
Hi Dave, That's rather frustrating. I've found vroom (from the package vroom) to be helpful with large files like this. Does the following give you any better luck? vroom(file_name, delim = ",", skip = 2459465, n_max = 5) Of course, when you know you've got errors & the files are big like that

[R] Exceptional slowness with read.csv

2024-04-08 Thread Dave Dixon
Greetings, I have a csv file of 76 fields and about 4 million records. I know that some of the records have errors - unmatched quotes, specifically.  Reading the file with readLines and parsing the lines with read.csv(text = ...) is really slow. I know that the first 2459465 records are good.

Re: [R] How to set the correct libomp for R

2024-04-08 Thread Ivan Krylov via R-help
В Mon, 8 Apr 2024 10:29:53 +0200 gernophil--- via R-help пишет: > I have some weird issue with using multithreaded data.table in macOS > and I am trying to figure out, if it’s connected to my libomp.dylib. > I started using libomp as stated here: > https://mac.r-project.org/openmp/ Does the

[R] How to set the correct libomp for R

2024-04-08 Thread gernophil--- via R-help
Hey everyone, I have some weird issue with using multithreaded data.table in macOS and I am trying to figure out, if it’s connected to my libomp.dylib. I started using libomp as stated here: https://mac.r-project.org/openmp/   Everything worked fine till beginning of this year, but all of a

Re: [R] Question regarding reservoir volume and water level

2024-04-07 Thread javad bayat
Dear all; Many thanks for your replies. This was not homework. I apologize. Let me explain more. There is a dam constructed in a valley with the highest elevation of 1255 m. The area of its reservoir can be calculated by drawing a polygon around the water and it is known. I have the Digital

Re: [R] Question regarding reservoir volume and water level

2024-04-07 Thread avi.e.gross
John, Your reaction was what my original reaction was until I realized I had to find out what a DEM file was and that contains enough of the kind of depth-dimension data you describe albeit what may be a very irregular cross section to calculate for areas and thence volumes. If I read it

Re: [R] duplicated() on zero-column data frames returns empty

2024-04-07 Thread Mark Webster via R-help
With respect to duplicated.data.frame taking account of row names to return all the rows as unique: thinking about this some more, I can see that making sense in isolation, but it's at odds with the usual behaviour of duplicated for other classes, e.g. primitive vectors, where it doesn't take

Re: [R] Question regarding reservoir volume and water level

2024-04-07 Thread avi.e.gross
Chris, since it does indeed look like homework, albeit a deeper looks suggests it may not beI think we can safely answer the question: >Is there any way to write codes to do this in R? The answer is YES. And before you ask, it can be done in Python, Java, C++, Javascript, BASIC, FORTRAN and

Re: [R] Question regarding reservoir volume and water level

2024-04-07 Thread Sorkin, John
Aside from the fact that the original question might well be a class exercise (or homework), the question is unanswerable given the data given by the original poster. One needs to know the dimensions of the reservoir, above and below the current waterline. Are the sides, above and below the

Re: [R] Question regarding reservoir volume and water level

2024-04-07 Thread Rui Barradas
Às 13:27 de 07/04/2024, javad bayat escreveu: Dear all; I have a question about the water level of a reservoir, when the volume changed or doubled. There is a DEM file with the highest elevation 1267 m. The lowest elevation is 1230 m. The current volume of the reservoir is 7,000,000 m3 at 1240

Re: [R] Question regarding reservoir volume and water level

2024-04-07 Thread Chris Ryan via R-help
Homework? -- Sent from my Android device with K-9 Mail. Please excuse my brevity. On April 7, 2024 8:27:18 AM EDT, javad bayat wrote: >Dear all; >I have a question about the water level of a reservoir, when the volume >changed or doubled. >There is a DEM file with the highest elevation 1267 m.

[R] Question regarding reservoir volume and water level

2024-04-07 Thread javad bayat
Dear all; I have a question about the water level of a reservoir, when the volume changed or doubled. There is a DEM file with the highest elevation 1267 m. The lowest elevation is 1230 m. The current volume of the reservoir is 7,000,000 m3 at 1240 m. Now I want to know what would be the water

Re: [R] duplicated() on zero-column data frames returns empty

2024-04-07 Thread Ivan Krylov via R-help
В Fri, 5 Apr 2024 16:08:13 + Jorgen Harmse пишет: > if duplicated really treated a row name as part of the row then > any(duplicated(data.frame(…))) would always be FALSE. My expectation > is that if key1 is a subset of key2 then all(duplicated(df[key1]) >= > duplicated(df[key2])) should

[R-es] Ayuda con kproto- Desviacion tipica

2024-04-06 Thread Antonio Rodriguez Andres
Hola usuarios de R Cuando aplico k prototype clustering <- kproto(x = data, k = 3, verbose = TRUE, lambda = 2) El objeto clustering no es un data frame, y summary(clustering) me da las medias para variables numéricas pero me gustaría obtener la desviación típica. Una forma que podría ser es

Re: [R] duplicated() on zero-column data frames returns empty

2024-04-05 Thread Jorgen Harmse via R-help
(I do not know how to make Outlook send plain text, so I avoid apostrophes.) For what it is worth, I agree with Mark Webster. The discussion by Ivan Krylov is interesting, but if duplicated really treated a row name as part of the row then any(duplicated(data.frame(�))) would always be FALSE.

Re: [R] duplicated() on zero-column data frames returns empty vector

2024-04-05 Thread Mark Webster via R-help
Hello Ivan, thanks for this. > Part of the problem is that it's not obvious what should be a > zero-column but non-zero-row data.frame mean. >  > On the one hand, your database relation use case is entirely valid. On > the other hand, if data.frames are considered to be tables of data with >

Re: [R] duplicated() on zero-column data frames returns empty vector

2024-04-05 Thread Ivan Krylov via R-help
Hello Mark, В Fri, 5 Apr 2024 03:58:36 + (UTC) Mark Webster via R-help пишет: > I found what looks to me like an odd edge case for duplicated(), > unique() etc. on data frames with zero columns, due to duplicated() > returning a zero-length vector for them, regardless of the number of >

[R] duplicated() on zero-column data frames returns empty vector

2024-04-05 Thread Mark Webster via R-help
Hello, I found what looks to me like an odd edge case for duplicated(), unique() etc. on data frames with zero columns, due to duplicated() returning a zero-length vector for them, regardless of the number of rows: df <- data.frame(a = 1:5)df$a <- NULLnrow(df) # 5 (row count preserved by

Re: [R] split a factor into single elements

2024-04-02 Thread Bert Gunter
Note: > levels(factor(c(0,0,1))) ## just gives you the levels attribute [1] "0" "1" > as.character(factor(c(0,0,1))) ## gives you the level of each value in the vector [1] "0" "0" "1" Does that answer your question or have I misunderstood. Cheers, Bert On Tue, Apr 2, 2024 at 12:00 AM Kimmo

Re: [R] split a factor into single elements

2024-04-02 Thread Ebert,Timothy Aaron
Using levels rather than length might cause problems. 2024 1, 1, 0, 0 will have a different number of levels than 2024, 3, 8, 0, 0 and I cannot assume that the two tailing zeros are zero for all records. The code can be simplified if you can assume more. It might require more work if I have

Re: [R] How to tweak genomic plot with genoPlotR?

2024-04-02 Thread Luigi Marongiu
Already did... On Tue, Apr 2, 2024 at 10:45 AM Eric Berger wrote: > > According to https://cran.r-project.org/web/packages/genoPlotR/index.html > the maintainer of genoPlotR is > > Lionel Guy > > Send your question also to him. > > On Tue, Apr 2, 2024 at 11:27 AM Luigi Marongiu > wrote: > > >

Re: [R] How to tweak genomic plot with genoPlotR?

2024-04-02 Thread Eric Berger
According to https://cran.r-project.org/web/packages/genoPlotR/index.html the maintainer of genoPlotR is Lionel Guy Send your question also to him. On Tue, Apr 2, 2024 at 11:27 AM Luigi Marongiu wrote: > > I would like to use your genoPlotR package > (doi:10.1093/bioinformatics/btq413) to

[R] How to tweak genomic plot with genoPlotR?

2024-04-02 Thread Luigi Marongiu
I would like to use your genoPlotR package (doi:10.1093/bioinformatics/btq413) to compare the genomes of two isolates of E. coli K-12 that I have. One is a K-12 that was in my lab's fridge; the other is a derivative of K-12 bought some time ago, HB101. I tried to use genoPlotR, but I could not

Re: [R] split a factor into single elements

2024-04-02 Thread Kimmo Elo
Hi, why would this simple procedure not work? --- snip --- mydf <- data.frame(id_station = 1234, string_data = c(2024, 12, 1, 0, 0), rainfall_value= 55) mydf$string_data <- as.factor(mydf$string_data) values<-as.integer(levels(mydf$string_data)) for (i in 1:length(values)) { 

[ESS] showing ESS session input/output and plots in a web browser

2024-04-01 Thread Jason E. Aten via ESS-help
My rbook project may be of interest to other ESS users. It is a project I wrote in Go (golang) that mirrors the R session (text commands, text outputs, and plots) to a web browser. This makes it easy for me to work on a remote linux box but see all the graphics (without having to run a local X

Re: [ESS] ess-set-working-directory

2024-03-31 Thread Kevin Coombes via ESS-help
That did it. Thanks, Kevin On Sun, Mar 31, 2024, 7:17 AM Stephen J. Eglen wrote: > > > hi Kevin, Rodney, > > I think the variable and setting you are looking for is: > > (setq ess-startup-directory 'default-directory) > > Best wishes, > > Stephen > > > > > > > > On Sun, Mar 31 2024, Kevin

Re: [ESS] ess-set-working-directory

2024-03-31 Thread Stephen J. Eglen via ESS-help
hi Kevin, Rodney, I think the variable and setting you are looking for is: (setq ess-startup-directory 'default-directory) Best wishes, Stephen On Sun, Mar 31 2024, Kevin Coombes via ESS-help wrote: I don't know where this "feature" lives in the code. But it might help find it to

Re: [ESS] ess-set-working-directory

2024-03-31 Thread Kevin Coombes via ESS-help
I don't know where this "feature" lives in the code. But it might help find it to know that it does the same thing in any git project or subversion project, even if it isn't an R package. Hopefully, someone can track it down, since I view it as a bug I want to remove, and not a feature. Best,

<    1   2   3   4   5   6   7   8   9   10   >