For interested people I got an answer on this link, thank you.
https://stackoverflow.com/questions/79573847/ggplot-trouble-when-i-combine-a-layer-with-group-and-a-layer-with-group-1?noredirect=1#comment140334381_79573847
Le 12/04/2025 à 17:24, Laurent Rhelp a écrit :
Dear RHelp-list,
I
Dear RHelp-list,
I want to combine with ggplot2 a layer where I compare
three regressions and a layer where I just display three
new samples to evaluate their locations to the regressions lines.
(the mock data is placed at the end of the mail)
# 1. comparison of the three regressions
20 a 105.4257 145.2059 20 b 166.4094 233.7112 20 c
Le 16/06/2024 à 20:54, Gabor Grothendieck a écrit :
> This can be vectorized. Try
>
> ix <- seq_along(vec2)
> S_diff2 <- sapply(seq_len(N1-(N2-1)*ratio_sampling), \(j)
> sum((vec1[(ix-1)*ratio_sampling+j] - vec2[ix])**2))
>
>
e into a set of instructions in more native R.
>> Just for comparison, in python, things like comprehensions for list or
> dictionaries or tuples often are syntactic sugar and the interpreter may
> simply rewrite them more like the first program you typed and evaluates
> that. The
Dear RHelp-list,
I try to use the package comprehenr to replace a for loop by a list
comprehension.
I wrote the code but I certainly miss something because it is very
slower compared to the for loops. May you please explain to me why the
list comprehension is slower in my case.
Here
p On Behalf Of Bert Gunter
Sent: Friday, May 24, 2024 6:52 PM
To: Laurent Rhelp
Cc: r-help@r-project.org
Subject: Re: [R] dplyr, group_by and selective action according to each group
Laurent:
As I don't use dplyr, this won't help you, but I hope you and others may
find it entertaining
Dear RHelp-list,
Using dplyr and the group_by approach on a dataframe, I want to be
able to apply a specific action according to the group name. The code
bellow works, but I am not able to write it in a more esthetic way using
dplyr. Can somebody help me to find a better solution ?
Thank
urent
Le 21/05/2024 à 12:15, Ivan Krylov a écrit :
> В Tue, 21 May 2024 11:29:33 +0200
> Laurent Rhelp пишет:
>
>> pos <- seek(con_in,2,origin="start")
>> # We have to repeat the command to return the good amount of read
>> # bytes
>> print(paste0(&qu
Dear RHelp-list,
I want to write at a specific position in a file without reading all
the file
because it is very large and I cannot read it in my RAM. But I miss
something
about the use of the command writeBin in conjunction with seek. In the
example bellow the seek
commands works well with
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
ls on the right as well, you will just need
to uncomment the line
ans$right <- TRUE
in your my_yscale.components.subticks function.
-Deepayan
On Tue, May 30, 2023 at 12:39 AM Laurent Rhelp wrote:
Dear R-Help-list,
I want to display many sensors on the same page so I have to
Dear R-Help-list,
I want to display many sensors on the same page so I have to adapt
the size of the y-axis labels and I woul like to adapt the number of
ticks according to the sensor. I use the yscale.components argument with
the function yscale.components.subticks: see the code below.
I
It is great !
I had not thought to use the value 0 for the variables style1 and
style2. I did not well understand when reading the documentation.
Thank you
Laurent
Le 13/02/2023 à 19:04, Deepayan Sarkar a écrit :
On Mon, Feb 13, 2023 at 3:48 AM Laurent Rhelp wrote:
Dear R-Help-list,
I
Dear R-Help-list,
I want to use the doubleYScale function from latticeExtra to overlap two
lattice graphs (cf. code below). The overlapping works but I lose the
groups of
every lattice, there are only two colors. Reading the documentation,
the arguments style1 and style2 give me the impression
labels do not
appear on the y-axis if I keep ylim = c(0,max(lim)). But nevermind
because your solution to use prepanel works very well and I will think
next time to use prepanel.
Best regards
L.
Le 22/01/2023 à 18:43, Deepayan Sarkar a écrit :
>
>
> On Sun, Jan 22, 2023 at 6:48 PM Laur
Dear RHelp-list,
I want to choice my scale for every panel in a lattice graph according to
the ylim range of every panel, with 10 ticks and with a start from 0.
Also I want to plot a grid
according to the y ticks (I did that in the panel argument with the
panel.abline function) .
So I
Thank you very much Deepayan, I will see the help of panel.bwplot.
Best regards
L.
Le 11/12/2022 à 18:53, Deepayan Sarkar a écrit :
On Sun, Dec 11, 2022 at 11:05 PM Laurent Rhelp wrote:
Excellent, it works.
But, may you please explain me how xyplot knows that it has to apply
panel.bwplot on
.
Le 11/12/2022 à 17:48, Deepayan Sarkar a écrit :
On Sun, Dec 11, 2022 at 2:33 PM Laurent Rhelp wrote:
I understand the idea but I did not succeed.
Here is what I tried:
## 1.middles of classes calculation
m <- tapply(DF$x, groups, mean)
## 2. create a new factor columns with
om
> , y = list(log=T, tck = c(1,0))
> , x = list(at = 1:5, tck = 1)
> )
>
> If this is not what you meant, you'll need to clarify ... or just move on.
>
> -- Bert
>
> On Sun, Dec 11, 2022 at 1:06 AM Laurent Rhelp
> wrote:
>
>
m). The potential problem with that is the box
> widths,
> which panel.bwplot() will not know how to compute.
>
> See if the following gives you a reasonable starting point:
>
> DF <- within(DF, m <- tapply(y, groups, mean))
> xyplot(
in(DF, m <- tapply(y, groups, mean))
xyplot(y ~ m, DF, scales = list(log = TRUE),
panel = panel.bwplot, horizontal = FALSE,
box.width = .0001)
Best,
-Deepayan
On Sat, Dec 10, 2022 at 7:46 PM Laurent Rhelp wrote:
Dear R-Help list,
I would like to use bwplot from the lattice
Dear R-Help list,
I would like to use bwplot from the lattice package with a log scale
both on
the x-axis and the y-axis but I do not know how to do that because I do
not know
how to change the factor x-axis in a numeric x-axis.
Here is my example:
library(lattice)
# the mock data
y <-
Very effective solution, I hope I remember that for the nex time.
Thank you
Le 09/07/2021 à 19:50, David Winsemius a écrit :
On 7/9/21 10:40 AM, Laurent Rhelp wrote:
Dear R-Help-list,
I have a list init_l containing 16 dataframes and I want to create
a matrix 4 x 4 from this list with
,"X2","X3","X4")))
> It doesn't give exactly what your code does, but your code introduces
> an extra level of "list", which you may not want.
>
> -Bill
>
> On Fri, Jul 9, 2021 at 10:40 AM Laurent Rhelp <mailto:laurentrh...@fre
Dear R-Help-list,
I have a list init_l containing 16 dataframes and I want to create a
matrix 4 x 4 from this list with a dataframe in every cell of the
matrix. I succeeded to do that but my loop is very uggly (cf. below).
Could somebody help me to write nice R code to do this loop ?
Thank
It is good trick and we have to know that R arrays use the 'column major
mode' to store the value.
Thank you
Laurent
Le 05/07/2021 à 19:03, David Winsemius a écrit :
On Jul 5, 2021, at 7:56 AM, Laurent Rhelp wrote:
Dear R-Help,
I have an array x made up of three matrices of 5
[i,,])$d ))
HTH,
Eric
On Mon, Jul 5, 2021 at 5:56 PM Laurent Rhelp
wrote:
Dear R-Help,
I have an array x made up of three matrices of 5 rows and 3 columns
of
complex numbers (the complex numbers are not the problem)
## my array
x <- structure(c(5.6196790161893828+0i, 5.7565523942393364
ope this helps!
>
> On Mon, Jul 5, 2021 at 10:56 AM Laurent Rhelp <mailto:laurentrh...@free.fr>> wrote:
>
> Dear R-Help,
>
> I have an array x made up of three matrices of 5 rows and 3
> columns of
> complex numbers (the complex numbers are not the prob
= FALSE)
>
> seems to work well for me. I hope this helps!
>
> On Mon, Jul 5, 2021 at 10:56 AM Laurent Rhelp <mailto:laurentrh...@free.fr>> wrote:
>
> Dear R-Help,
>
> I have an array x made up of three matrices of 5 rows and 3
> columns of
( svd(x[i,,])$d ))
>
> HTH,
> Eric
>
> On Mon, Jul 5, 2021 at 5:56 PM Laurent Rhelp <mailto:laurentrh...@free.fr>> wrote:
>
> Dear R-Help,
>
> I have an array x made up of three matrices of 5 rows and 3
> columns of
> complex numbers (the complex
Dear R-Help,
I have an array x made up of three matrices of 5 rows and 3 columns of
complex numbers (the complex numbers are not the problem)
## my array
x <- structure(c(5.6196790161893828+0i, 5.7565523942393364+0i,
8.5242834298729342+0i,
10.304766710160479+0i, 11.4129670101082
.w) <- as.POSIXlt(time(DF.w))$sec
# now run plot or xyplot as before
On Mon, Feb 22, 2021 at 11:56 AM Gabor Grothendieck
wrote:
I assume that this is a lattice problem. Replacing xyplot with plot
and using all the same arguments there is no error.
On Mon, Feb 22, 2021 at 11:26 AM Laurent Rhe
o suppress the X axis and then issue
an axis command. plot.zoo has examples.
On Mon, Feb 22, 2021 at 12:14 PM Laurent Rhelp wrote:
You are right, if I replace xyplot by plot there is no error message but
the POSIXct format is not taken into account for the labels on the
x-axis. Do I have to wri
Le 22/02/2021 à 18:01, Deepayan Sarkar a écrit :
On Mon, Feb 22, 2021 at 9:56 PM Laurent Rhelp wrote:
Dear R-Help-List,
I have to process time series with a sampling frequency of 1 MHz.
I use the POSIXct format for the date-times with microsecond in a zoo
object and the xyplot.zoo
Dear R-Help-List,
I have to process time series with a sampling frequency of 1 MHz.
I use the POSIXct format for the date-times with microsecond in a zoo
object and the xyplot.zoo function to do the graphs.
As I show in the below example I had a trouble to plot the labels on the
x-axis with
;,"","","","","","",""), flush=TRUE, multi.line=FALSE,
skip=nextElem(iter_index1)-1, nlines=1, quiet=TRUE))
[1] "N163" "-0.054023" "-0.049345" "-0.037158" "-0.04112"
"-0.044612&quo
| こんにちは | مرحبا | Привет /, :v ).put;'
hello_7lang.txt", intern = TRUE, ignore.stderr = FALSE))
[1] """"""
"你好 Chinese"
[5] "こんにちは Japanese" "مرحبا Arabic""Привет Rus
Laurent Rhelp wrote:
# sensors to keep
sensors <- c("N053", "N163")
# filter on the beginning of the line
i <- pipe("grep -E '^(N053|N163)' test.txt")
# or:
# filter on the beginning of the given column
# (use $2 for the second column, etc.)
i <- pi
;,"","","","","","",""), flush=TRUE, multi.line=FALSE,
skip=nextElem(iter_index1)-1, nlines=1, quiet=TRUE))
[1] "N163" "-0.054023" "-0.049345" "-0.037158" "-0.04112"
"-0.044612&quo
;","","","","",""),
flush=TRUE, multi.line=FALSE, skip=nextElem(iter_index1)-1, nlines=1, quiet=TRUE))
[1] "N163" "-0.054023" "-0.049345" "-0.037158" "-0.04112"
"-0.044612" &quo
15094 -0.012104
5 N163 -0.054023 -0.049345 -0.037158 -0.041120 -0.044612 -0.036953
-0.036061 -0.044516 -0.046436
Best, Bill.
W. Michels, Ph.D.
On Sun, May 17, 2020 at 5:43 PM Laurent Rhelp wrote:
Dear R-Help List,
I would like to use an iterator to read a file filtering some
selected lin
15094 -0.012104
5 N163 -0.054023 -0.049345 -0.037158 -0.041120 -0.044612 -0.036953
-0.036061 -0.044516 -0.046436
Best, Bill.
W. Michels, Ph.D.
On Sun, May 17, 2020 at 5:43 PM Laurent Rhelp wrote:
Dear R-Help List,
I would like to use an iterator to read a file filtering some
selected lin
Dear R-Help List,
I would like to use an iterator to read a file filtering some
selected lines according to the line name in order to use after a
foreach loop. I wanted to use the checkFunc argument as the following
example found on internet to select only prime numbers :
|
I have a spatial cross-sectional dataset. I think I have an instrument for
my dependent variable, but there's reason to suspect I need a spatial lag
with spatial error model. The dataset has about 100,000 observations.
I am looking for a package which will estimate my model in R (or Stata or
P
Dear List members,
I would be glad to get some tips on the following topic.
I use rgl package to generate a 3D-cube with labeled dots at certain
coordinates. After that I try to export the scene as a WebGL
application. For this I use the syntax from rgl documentation
"writeWebGL" mentioned in
Hi,
I?m fairly new to R and have a problem mentioned in the subject ...
I want to draw a scatterplot in 3d - either with scatterplot3d or -
preferably - with the rgl package - but instead of points or text
(text3d command of rgl) I would like to draw either histograms or pie
charts to visualiz
= ... code to each plot,
but I didn't try it.) If the legend is irrelevant to you, then this
should be fine.
HTH,
Dennis
On Sun, Jun 22, 2014 at 2:09 PM, Laurent Rhelp wrote:
Le 20/06/2014 21:50, Christoph Scherber a écrit :
Dear Laurent
for numeric x variables, you could try jitter:
xy
Le 20/06/2014 21:50, Christoph Scherber a écrit :
Dear Laurent
for numeric x variables, you could try jitter:
xyplot(y~jitter(x,0.5))
Cheers
Christoph
Am 20.06.2014 21:45, schrieb Laurent Rhelp:
Hi,
I like to use with xyplot (package lattice) the groups argument and
superpose.symbol to
Hi,
I like to use with xyplot (package lattice) the groups argument and
superpose.symbol to compare several curves. But, when there are a great
many points, the symbols are very close and the graph becomes
unreadable. Would there be an argument or a tip not to draw all the
symbols, for exa
Hi,
I have a weird issue when creating a package with latin1 encoding.
All works fine except the fact that I can't put non-ASCII characters
in the title tag. There is NA in the place of the title in the pdf file
of the documentation.
I can put non-ASCII characters in the other tags (desc
ards,
Ben
>
> From: Uwe Ligges
>To: Ben Rhelp
>Cc: "r-help@r-project.org"
>Sent: Tuesday, 24 April 2012, 19:30
>Subject: Re: [R] Compiling RMySQL on Win7 64bits RS-DBI.c:1:0: sorry,
>unimplemented: 64-bit mode not compiled in
>
Save\TempFolder\RtmpS8df20\downloaded_packagesâ
Warning messages:
1: running command 'D:/BenSave/R/R-2.15.0/bin/i386/R CMD INSTALL -l
"D:/BenSave/R/R-2.15.0/library"Â Â
C:\Users\BenSave\TempFolder\RtmpS8df20/downloaded_packages/RMySQL_0.9-3.tar.gz'
had status 1
2: In install.packages("RMySQL", type = &quo
Hi,
On a Win7 64 bit, I have installed:
*R-15.0 and I am running it in 64bits
*Rtools215.exe*MySQL 5.5.23 64bits
My PATH start with
D:\BenSave\Rtools\bin;D:\BenSave\Rtools\MinGW64\bin;D:\BenSave\R\R-2.15.0\bin;D:\BenSave\Rtools\MinGW\bin;D:\BenSave\Rtools\gcc-4.6.3\bin;
Note that I place MinGW
Hi Prof Brian,
Thank you for your email and for writing MASS. This book is brilliant.
- Original Message
> From: Prof Brian Ripley
> To: Ben Rhelp
> Cc: r-help@r-project.org
> Sent: Thu, 16 June, 2011 14:48:00
> Subject: Re: [R] Porting "unmaintained" pack
- Original Message
> From: Uwe Ligges
> To: Ben Rhelp
> Cc: r-help@r-project.org
> Sent: Thu, 16 June, 2011 14:38:12
> Subject: Re: [R] Porting "unmaintained" packages to post R 2.10.0 era
>
[...]
> >
> > What about "--binary is deprecate
Hi all,
I am trying to re-compile some "unmaintained" (it seems) packages, namely
rSoNIA
and dynamicnetwork from:
http://csde.washington.edu/~skyebend/rsonia/rsoniaDemo/
These packages predates R 2.10.0 so they need to be recompile.
After split the single big file in /man in each packages into
Le 07/05/2011 06:17, Penny Bilton a écrit :
I am trying to find a confidence band for a fitted non-linear curve. I
see that the predict.nls function has an interval argument, but a
previous post indicates that this argument has not been implemented.
Is this still true? I have tried various way
Hi Martin,
Thanks a lot for your help. this was spot on. With your instructions, I manage
to successfully compile Rgraphiz.
thanks again.
regards,
Ben
- Original Message
> From: Martin Morgan
> To: Ben Rhelp
> Cc: r-help@r-project.org
> Sent: Wed, 4 May, 2011 1:01:
Hi all,
I am trying to compile Rgraphiz on Windows 7 64bit with R-2.13.0. I have
installed
Rtools213.exe from [1]. The 64bit packages in [2] provided me with the 64 bit
version
of graphviz. After intalling the binary version Rgraphviz 1.30 (in 32bit) it
complains (as
expected) that:
> libra
Hi All,
I am trying to use a neural network for my work, but I am not sure about my
approach to select a parsimonious model. In R with nnet, the IAC has
not been defined for a feed-forward neural network with a single hidden layer.
Is this because it does not make sens mathematically in this cas
Hi,
I am trying to install the "ncdf" package on my Ubuntu 8.04 (Karmic)
box in the basement. I have installed the most recent distribution of
R, 2.12.0 without too many problems. However, when it comes to install
the "ncdf" netcdf package, I unfortunately get the following error:
Error in dyn.lo
Hi Terry, David, and Thomas,
Thank you for all your emails and the time you to took to clarify my
misunderstanding on survival analysis. I will need a bit of time to digest all
this information and to do some more reading.
Best regards,
Ben
> From: Terry Therneau
>
> 1. survreg() does N
Hi David,
Thank you for your reply. See below for more information.
> From: David Winsemius
>
> On Nov 25, 2010, at 7:27 AM, Ben Rhelp wrote:
>
> > I manage to achieve similar results with a Cox model as follows but I don't
> > really understand why we ha
I manage to achieve similar results with a Cox model as follows but I don't
really understand why we have to take the inverse of the linear prediction with
the Cox model and why we do not need to divide by the number of days in the
year
anymore?
Am I getting a similar result out of pure luck?
Hi all,
Is there an equivalent to predict(...,type="linear") of a Proportional hazard
model for a Cox model instead?
For example, the Figure 13.12 in MASS (p384) is produced by:
(aids.ps <- survreg(Surv(survtime + 0.9, status) ~ state + T.categ +
pspline(age, df=6), data = Aidsp))
zz <- predi
Hi all,
I am looking at the tutorial/appendix from John Fox on “Cox
Proportional-Hazards Regression for Survival Data” available here:
http://cran.r-project.org/doc/contrib/Fox-Companion/appendix-cox-regression.pdf
I am particularly interested in modelling survival with time-dependent
covariates
Dear R-List,
I try to use accented chars in a plot using odfWeave (see the example
below). I found the solution using the iconv command to write the
accented char "é". But, is somebody know another solution not to have to
use iconv command every time ?
Thanks
Example :
<>=
par(mfrow=c(1,
Le 25/03/2010 20:28, Laurent Rhelp a écrit :
Dear R-List,
I am working with binary data that I want to store in a PostgreSQL
DataBase. I decided to use a TXT field. I read my binary file with
readBin function, I succeed in my data storage in the database but I
have some trouble to extract
Dear R-List,
I am working with binary data that I want to store in a PostgreSQL
DataBase. I decided to use a TXT field. I read my binary file with
readBin function, I succeed in my data storage in the database but I
have some trouble to extract the data : the correct amount of bytes is
sto
of 0 secs
y <- x+.01 # 1 us
identical(x,y)
[1] FALSE
y <- x+.001 # 0.1 us
identical(x,y)
[1] TRUE
On Sat, Feb 6, 2010 at 9:41 AM, Laurent Rhelp wrote:
Gabor Grothendieck a écrit :
zoo is independent of time and date class so it does not restrict y
, 2010 at 8:25 AM, Laurent Rhelp wrote:
Dear R-List,
I have the habit of using R for my data processing and I like to use the
power of the lattice package. Presently, I have to manage time series. So,
in order to work properly I want to discover the zoo package and the related
methods (since
Dear R-List,
I have the habit of using R for my data processing and I like to use
the power of the lattice package. Presently, I have to manage time
series. So, in order to work properly I want to discover the zoo package
and the related methods (since lattice can work with zoo class). But m
Dear All I hope that someone can help.
I am working with sp pakage and akima
library("akima")
library(sp)
imagine lots of different dataframes, of row = 100 columns = 3 of x and y
coordinates with z values I will call these data frames for the sake of this
example akima
akima<-as.list(1:100)
David Kane a écrit :
I have a pdf file that I would like to parse into R:
http://www.williams.edu/Registrar/geninfo/faculty.pdf
For now, I open the file in Acrobat by hand, then save it "as text"
and then use readLines(). That works fine but a) I am concerned that
some information may be lost
the x axis corresponds to row
> number and the y axis to column number, with column 1 at the
> bottom, i.e. a 90 degree counter-clockwise rotation of the
> conventional printed layout of a matrix.
>
> [2] http://wiki.r-project.org/rwiki/doku.php?id=tips:spatial-data
>
> Rhelp wanted wrote:
Dear all.
I am using the akima function to produce 3d contour plots using interp based
on irregular data.
using the eg in the akima manual
library("akima")
data(akima)
plot(y ~ x, data = akima, main = "akima example data")
with(akima, text(x, y, formatC(z,dig=2), adj = -0.1))
## linear interpola
I have to start by saying that I am new to R, so I might miss something crucial
here. It seems to me that the results of friedman.test and ks.test are "wrong".
Now, obviously, the first thing which crossed my mind was "it can't be, this is
a package used by so many, someone should have observed"
Dear all,
The following code is inspired by the help file for the relist()
function (see?relist), which explicitly details how you can use a
relistable object in conjunction with optim to pass and reconstruct
complex parameter structures/groupings. The idea is that the optim()
function can only wo
Kiaora Rolf,
Thanks for the reply. The syntax that I am using is based on examples
from ?optim, where a very similar system is setup with the rosenbrook
banana function. Given that nls() is basically a wrapper for optim(),
it seems reasonable that the logic should carry across. Furthermore,
the fo
Dear R-help,
Could you please examine the following code, and see if I have discovered a bug
or not, or am just doing something silly.
I am trying to create a package to do fish stock assessment using the nls() function to
fit the modelled stock size to the various pieces of information that w
Dear R gurus,
I am trying to plot a legend in the margins of a figure. The basic idea is to
have two (or more) plots on the same figure, but then to have a common legend
at the bottom of the plot. The approach that I have taken is to setup the
figures, then do a "dummy" plot of the legend to de
81 matches
Mail list logo
