-Oorspronkelijk bericht-
Van: Enrico Schumann
Verzonden: woensdag 25 juni 2025 16:44
Aan: Heuvel, E.G. van den (Guido)
CC: 'r-help@R-project.org'
Onderwerp: Re: [R] Potential bug in readLines when reading empty lines
[Externe email]
Quoting "Heuvel, E.G. van de
, it appears to not be working
correctly.
-Oorspronkelijk bericht-
Van: Jeff Newmiller
Verzonden: woensdag 25 juni 2025 16:10
Aan: Heuvel, E.G. van den (Guido) ; Heuvel, E.G. van den
(Guido) via R-help ; 'r-help@R-project.org'
Onderwerp: Re: [R] Potential bug in readLines when re
ow should I read test.txt in such a way that the
empty 2nd line is left intact?
Best regards,
Guido van den Heuvel
Statistics Netherlands
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/
Hi,
In an MCMC algorithm, we noticed that our C++ code runs a lot slower on Windows
than on macOs or Linux (the algorithm takes ~20 seconds on Mac, and easily 2.5
minutes on Windows).
After profiling, it appeared that the main time difference was in the time
spent in calls to std::exp and std::
Hi R-help,
I am a researcher fitting linear mixed models using the package nlme. I am
wondering whether the sigma value that is outputted as part of the model
object is standard deviation or variance? for example, I might fit a model
lme_mod <- nlme::lme(response ~ predictor1 + predictor2, random
Many thanks to all who helped out! What a wonderful community is this!
Cheers!
Ben van den Anker
On Friday, December 25, 2020, 10:52:38 PM GMT+1, Abby Spurdle
wrote:
I'm not sure what the official status of Rmetrics is.
However, as far as I can see most of the Rme
Hello everyone,
Could anyonre recommend some good resources for finance applications in R? I
find the packages quantmod, TTR and PerformanceAnalytics a bit outdated. There
must be something more recent on the market. Any suggestions will be much
appreciated!
Cheers,
Ben van den Anker
On Mon, 20 Apr 2015 18:10:10 +0200, Marc Schwartz
wrote:
On Apr 20, 2015, at 10:01 AM, j. van den hoff
wrote:
hi,
consider the following example:
8<-
x <- matrix(1:6, 3, 2)
layout(1:2)
barplot(x, beside = TRUE, col = 1:6)
barplot(x, beside =
hi,
consider the following example:
8<-
x <- matrix(1:6, 3, 2)
layout(1:2)
barplot(x, beside = TRUE, col = 1:6)
barplot(x, beside = FALSE, col = 1:6)
8<-
it seems, it is not possible to make `beside=FAlSE' plots behave the
On Sat, 09 Nov 2013 18:18:23 +0100, Duncan Murdoch
wrote:
On 13-11-09 11:57 AM, j. van den hoff wrote:
On Sat, 09 Nov 2013 17:16:28 +0100, Duncan Murdoch
wrote:
On 13-11-09 8:50 AM, j. van den hoff wrote:
I'd very much appreciate some help here: I'm in the process of
clarifyi
On Sat, 09 Nov 2013 17:16:28 +0100, Duncan Murdoch
wrote:
On 13-11-09 8:50 AM, j. van den hoff wrote:
I'd very much appreciate some help here: I'm in the process of
clarifying
whether I can use `computeContour3d' to derive estimates of the surface
area of a single closed
I'd very much appreciate some help here: I'm in the process of clarifying
whether I can use `computeContour3d' to derive estimates of the surface
area of a single closed isosurface (and prospectively the enclosed
volume). getting the surface area from the list of triangles returned by
`comp
On Fri, 08 Nov 2013 15:32:00 +0100, Barry Rowlingson
wrote:
On Fri, Nov 8, 2013 at 1:01 PM, j. van den hoff
wrote:
ps: it of course would still be nice, if the surface area (or a vector
containing the individual triangle areas)
were returned to the caller as well ...
Does the
regarding my previous mail for this topic, I have in the meantime
identified my misconception.
actually, `computeContour3d' returns the vertices just fine in the correct
coordinate frame.
the misconception was caused basically by assuming that the `level'
argument was a fractional
thresh
I want to compute the total surface area of an isosurface in 3D space as
approximated with
`computeContour3d' by adding up the areas of the triangles making up the
contour surface.
problem:
the vertex matrix returned by `computeContour3d' in general seems to
provide the vertices
not in the fr
Dear R experts,
I am trying to simulate a list containing data matrices. Unfortunately, I don't
manage to get it to work.
A small example:
n=5
nbootstrap=2
subsets<-list()
for (i in 1:n){
subsets[[i]] <- rnorm(5, mean=80, sd=1)
for (j in 1:nbootstrap){
test<-list()
Dear R-users,
For obtaining a weighted correlation matrix with (some) missing data I use the
sweep operator (swp function in the ggm-package) to find the maximum likelihood
estimates of the standardized slopes and error variances. In the end I would
like to use a reverse sweep operator but unfo
nsectSprays, col = "lightgray", yaxt='n')
>> axis(2, at=seq(0,25,2.5))
>> boxplot(count ~ spray, data = InsectSprays, col = "lightgray", yaxt='n')
>> axis(2, at=seq(0,25,1))
>>
>> HTH,
>> Pascal
>>
>>
>> On 29/03/1
ias
-Oorspronkelijk bericht-
Van: Pascal Oettli [mailto:kri...@ymail.com]
Verzonden: vrijdag 29 maart 2013 8:44
Aan: Berg, Tobias van den
CC: 'r-help@r-project.org'
Onderwerp: Re: [R] changing y-axis intervals in a boxplot
Hello,
?par (particularly the 'yaxt' argumen
Dear R-users,
I'm a bit of a rooky in R and do not know how to change the y-axis intervals of
a boxplot.
The y-axis in my boxplot ranges from 0 to 100 which is fine. The problem is the
interval level in between. Now it increases with steps of 20 but I'd like to
have it changed by steps of 5.
[dwinsem...@comcast.net]
Verzonden: vrijdag 7 september 2012 18:17
Aan: Berg, Tobias van den
CC: PIKAL Petr; r-help
Onderwerp: Re: [R] error: in catg (xi, name=nam, label=lab): "LO2" has <2
category levels
On Sep 7, 2012, at 8:03 AM, Berg, Tobias van den wrote:
> Dear all,
>
&g
quot;nee geen atrofie",..: 1 1 1 1 1 1 1 1 1 1 ...
[[1]]
NULL
[[2]]
NULL
[[3]]
NULL
[[4]]
NULL
[[5]]
NULL
[[6]]
NULL
[[7]]
NULL
[[8]]
NULL
[[9]]
NULL
[[10]]
NULL
It would be great to receive ideas, comments or questions about my challenge.
Kind regards, Tobias
-Oorspronkelijk ber
hi list,
used versions: 2.12.1 and 2.14.0 under ubuntu and macosx.
I recently stumbled over a problem with `nls', which occurs if the model
is not specified explicitly but via an evaluation of a 'call' object.
simple example:
8<--
or message:
Error in "DGL growth (cm.yr"^{ : non-numeric argument to binary operator
I used the same line for the plot function in another graphic and it worked
perfectly well.
Any suggestions?
--
Prof. Eduardo van den Berg
Dpto. de Biologia - Universidade Federal de Lavras
Cx. P. 303
Hi Hadley,
have a look at:
http://www.riani.it/pub/zrc-csda98.pdf
and some Gauss code:
http://www.riani.it/Gauss/procedures/BOXPLOTB.G
Best regards,
Daniele
2011/9/21 Hadley Wickham :
> Hi all,
>
> Does anyone have an R implementation of the queplot (K. M. Goldberg
> and B. Iglewicz. Bivariate
to get a more realistic timing.
Any comment or idea to further speed up multivariate time series
creation with classes xts or timeSeries starting from a data.frame
like the one reported here is welcome.
Best regards,
Den
a data.frame example (code below to generate it)
ID DATE
to get a more realistic timing.
Any comment or idea to further speed up multivariate time series
creation with classes xts or timeSeries starting from a data.frame
like the one reported here is welcome.
Best regards,
Den
a data.frame example (code below to generate it)
IDDATE
lt;- daply(X, "ID", buildXtsFromDataFrame,
env = e1)
return(get("xx", e1))
}
Time04 <- replicate(100,
system.time(xtsDaply(X, X$ID))[[1]])
2011/4/4 Joshua Ulrich :
> Hi Dan,
>
> On Mon, Apr 4, 2011 at 7:49 AM, Den Alpin wrote:
>> I retrieve
Daply <- function(...)
{
library(plyr)
e1 <- new.env(parent = baseenv()) #create a new env
res <- daply(X, "ID", buildTimeSeriesFromDataFrame,
env = e1)
return(get("xx", e1)) # return xx from env
}
##replicate 100 times
#Time03 <- replicate(100,
# sy
Hate to say that, but it looks like Stata is way above R, considering
exact logistic regression.
To use elrm() I have to aggregate my data,which is really time consuming
when I look for the way out through many variables. But the worst thing
is that I am not not sure if I can trust to p-values in
'Aggregate multiple functions into a single function. Combine multiple
functions to a single function returning a named vector of outputs'
This is a short description of each() function from plyr package
Here is an example from help
each(min, max)(1:10)
Hope this helps
Regards
Denis
У Пят
logistic regression by MCMC. I found this as the first set of hits I
> got from
>
> # install.packages('sos') # if necessary
> library(sos)
> findFn('exact logistic regression')
>
> HTH,
> Dennis
>
> On Wed, Feb 2, 2011 at 8:43 AM, Den wrote:
&g
Hello to R people
Does anybody know to calculate exact logistic regression in R?
Does such option exist anywhere?
Surprisingly, could not find it using search engine. It is hard to
believe, however, that such useful function is not implemented in R yet?
Could you help, please
Thank you
Denis
max(subset(dataName, condition=="GPR119a")$responce)
У Аўт, 01/02/2011 у 19:05 +, A. Ramzan піша:
> Hello
>
> I am trying to find a way to find the max value, for only a subset of a
> dataframe, depending on how the data is grouped for example,
>
> How would I find the maxmium responce, for
Maybe that:
su <- lapply(dat[2:5],function(x)table(x))
su
mode(su)
myBYdata <- data.frame( do.call(cbind,lapply(su, as.data.frame)) )
myBYdata
У Няд, 23/01/2011 у 07:43 -0500, Wei Yang піша:
> Dear all,
>
> I would like to ask whether anyone has experience with the problem below.
>
>
> I want
And afterwards those NAs are removed by sort command.
It is a lot easier for me to deal with data when I know what I am doing.
Thank you again for help. Sorry for annoying naive questions.
With best regards
Denis Kazakiewicz
Belarus
Няд, 23/01/2011 у 05:38 -0800, P Ehlers піша:
> Den wrot
Hi
So let's call first data frame df1 and second df2
#replace NA with 2.36
df2y <- lapply(df2, function(x) (ifelse( is.na(x), 2.36,x)))
#combibe
both <- rbind(df1, df2)
With best regards
Denis
У Няд, 23/01/2011 у 12:41 +0700, Jumlong Vongprasert піша:
> Dear All.
> I have some problem with com
ut not broken Belarus
У Суб, 22/01/2011 у 17:55 -0800, Dennis Murphy піша:
> Hi:
>
> I wouldn't pretend to speak for Henrique, but I'll give it a shot.
>
> On Sat, Jan 22, 2011 at 4:44 AM, Den wrote:
> Dear R community
> Recently, dear Henriqu
...One more important thing:
After updating your packages with sudo
if you have additional packages installed with
'install.packages("", dependencies=TRUE)' in R
you will have to reinstall them
Regards
Denis
У Суб, 22/01/2011 у 12:20 +0530, Ajay Ohri піша:
> Dear List
>
> I use synaptic to dow
Dear R community
Recently, dear Henrique Dallazuanna literally saved me solving one
problem on data transformation which follows:
(n_, _n, j_, k_ signify numbers)
SOURCE DATA:
id cycle1 cycle2 cycle3 … cycle_n
1 c c c c
1 m m m
erely yours
Denis Kazakiewicz
Belarus
У Пят, 21/01/2011 у 18:37 -0200, Henrique Dallazuanna піша:
> Just change the FUN function:
>
> aggregate(.~ id, lapply(df, as.character), FUN =
> function(x)paste(sort(x), collapse = ''), na.action = na.pass)
>
> On Fri, Jan 21, 20
Hello
Consider following dataframe named df
var1var2var3
3771354 565
654654 963 6677
775 147 657754
df <- read.table('clipboard', header = TRUE)
df
#find indexes with '77' in var 1
myIndexes <- grep( glob2rx("*77*"), df$var1)
myIndexes
#find actual values of seach ab
ot;", na.action =
> na.pass)
>
>
>
> On Fri, Jan 21, 2011 at 5:45 PM, Den
> wrote:
> Dear Henrique
> Thank you again for helping me
> Unfortunately, your code see
y this:
>
> aggregate(.~ id, lapply(test, as.character), FUN = paste, collapse =
> "")
>
> On Fri, Jan 21, 2011 at 10:25 AM, Den wrote:
> Dear [R] people
> Could you please help with following data transformation.
> Any suggestions, hints, refer
Dear [R] people
Could you please help with following data transformation.
Any suggestions, hints, references and even guessing on performing any
of the following steps are highly appreciated. Those transformations are
crucial for my work.
(n_, _n, j_, k_ signify numbers)
SOURCE DATA:
id
Dear R people
Could you please help.
Basically, there are two variables in my data set. Each patient ('id')
may have one or more diseases ('diagnosis'). It looks like
id diagnosis
1 ah
2 ah
2 ihd
2 im
3 ah
3 stroke
4 ah
4 ihd
4 angina
5
Hi,
I'd like to select one row in a data frame per subset which is maximal for a
particular value. I'm pretty close to the solution in the sense that I can
easily select the maximal values per subset using "aggregate", but I can't
really figure out how to select the rows in the original data fram
I get the error
Error : evaluation nested too deeply: infinite recursion /
options(expressions=)?
during a 'R CMD check ...'
on one of my packages. The reason seems to be that this package is
mutually dependent on another one (i.e. the DESCRIPTION files of package
A lists package B under "Depen
this is probably not really a R specific question, if so apologies for
off-topic posting:
I'm interested in the probability density function of the maximum values
from repeated samples of size N from a normal distribution:
smp <- rnorm(N, meanval, stdev)
with some mean 'meanval' and standard dev
__
R-help@r-project.org mailing list
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.
lue","red"),
cex = .75),
space = "bottom", border = T),main=species.vector[i]) )
}
print(Lp, split=c(1,1,2,3), more=TRUE)
print(Pp, split=c(2,1,2,3), more=TRUE)
print(Pl, split=c(1,2,2,3), more=TRUE)
print(Ra, split=c(2,2,2,3), more=TRUE)
print(Lc, split=c(1,
Dear R users,
I am trying to reverse the orientation of axis labels and title in the right
margin of a plot, so that they read from top to bottom. I know that this can be
done using text() as follows:
par(mar=c(5,4,4,4)+.1)
plot(1:4,las=0)
par(new=T)
y <- rnorm(4)
plot(y,axes=FALSE,ann=FALSE,
On Mon, Jul 07, 2008 at 11:15:57AM +0200, Thiemo Schreiber wrote:
> Hello everyone,
>
> I have biological data from a competition experiment where a free ligand is
> titrated against the binding of a protein.
>
> Now, I would like to fit a standard on-site binding curve to this data in
> order
On Tue, Apr 29, 2008 at 01:49:01PM +0200, Philipp Pagel wrote:
>
> > When exporting to PDF a graph with a legend, in the final PDF, the
> > text is going beyond the legend box.
> > > dev2bitmap("test.pdf", type="pdfwrite", h=6, w=6)
> > The legend looks OK on the screen. I noticed that the size
dear list,
my question concerns the use of `eval' in defining the model formula
for `nls' (version 2.6.2.).
consider the following simple example, where the same model and data
are used to perform unweighted and weighted fits. I intentionally
used very uneven weights to guarantee large difference
On Wed, Apr 16, 2008 at 03:56:26PM -0600, Matthew Keller wrote:
> Yes Chuck, you're right.
>
just a comment:
> Thanks for the help. It was a data.frame not a matrix (I had called
> as.matrix() in my script much earlier but that line of code didn't run
> because I misnamed the object!). My bad. T
have not followed the thread completely, but:
have you tried `bitmap' with `type = pdfwrite' (or psgrb)
for comparison? at least with `pdf' there are some issues which
can be avoided by using ghostscript via `bitmap'.
joerg
On Mon, Mar 31, 2008 at 04:17:50PM -0400, Francois Pepin wrote:
> Prof B
On Wed, Mar 26, 2008 at 05:27:22PM +1300, Frank Scherr wrote:
> Hello together,
>
> Is there a tool to test the statistical differences between parameter
> estimates of a nlsList fit?
>
> I fitted degradation data using the nlsList method and want to find out
> whether derived rate constants
On Sun, Dec 02, 2007 at 11:08:01AM -0800, Milton Cezar Ribeiro wrote:
> Dear all,
> I am still fighting against my "power model".
> I tryed several times to use nls() but I can??t run it.
> I am sending my variables and also the model which I would like to fit.
> As you can see, this "power model"
On Mon, Nov 26, 2007 at 04:35:53PM +0100, Niccolò Bassani wrote:
> Dear R-users,
> I'm posting a problem I already asked help for some time ago, because I'm
> facing that problem once again and even because now, reading that old
> e-mail, and the answer recevied, I understand I've not made myself c
On Tue, Nov 13, 2007 at 09:04:25AM -0500, Duncan Murdoch wrote:
> On 11/13/2007 8:39 AM, Joerg van den Hoff wrote:
> >my understanding of `eval' behaviour is obviously
> >incomplete.
> >
> >myquestion: usually `eval(expr)' and `eval(expr
my understanding of `eval' behaviour is obviously
incomplete.
myquestion: usually `eval(expr)' and `eval(expr,
envir=parent.frame())' should be identical. why does the
last `eval' (yielding `r3') in this code _not_ evaluate in
the local environment of function `f' b
On Mon, Nov 12, 2007 at 06:59:56PM -0500, Duncan Murdoch wrote:
> On 12/11/2007 2:56 PM, Joerg van den Hoff wrote:
> >On Mon, Nov 12, 2007 at 11:09:21AM -0500, Duncan Murdoch wrote:
> >>On 11/12/2007 9:14 AM, Joerg van den Hoff wrote:
> >>>On Mon, Nov 12, 2007 at 07
On Mon, Nov 12, 2007 at 04:27:26PM -0300, Horacio Castellini wrote:
> Hi all, I'd like fit this function:
>
> G(t)=A*1/(1+t/B)*1/sqrt(1+t/C)
>
> where A, B and C are fitting constants that one like search. I have
> got a fcs<-(t,G(t)) value table in a data frame.
>
> Any one can help me? Tahnks
On Mon, Nov 12, 2007 at 11:09:21AM -0500, Duncan Murdoch wrote:
> On 11/12/2007 9:14 AM, Joerg van den Hoff wrote:
> >On Mon, Nov 12, 2007 at 07:36:34AM -0500, Duncan Murdoch wrote:
> >>On 11/12/2007 6:51 AM, Joerg van den Hoff wrote:
> >>>I initially thought, thi
On Mon, Nov 12, 2007 at 03:25:38PM +0100, Martin Maechler wrote:
> >>>>> "DM" == Duncan Murdoch <[EMAIL PROTECTED]>
> >>>>> on Mon, 12 Nov 2007 07:36:34 -0500 writes:
>
> DM> On 11/12/2007 6:51 AM, Joerg van den Hoff wrote:
&g
On Mon, Nov 12, 2007 at 07:36:34AM -0500, Duncan Murdoch wrote:
> On 11/12/2007 6:51 AM, Joerg van den Hoff wrote:
> >I initially thought, this should better be posted to r-devel
> >but alas! no response.
>
> I think the reason there was no response is that your example
I initially thought, this should better be posted to r-devel
but alas! no response. so I try it here. sory for the
lengthy explanation but it seems unavoidable. to quickly see
the problem simply copy the litte example below and execute
f(n=5)
which crashes. called with n != 5 (and of co
On Sat, Oct 13, 2007 at 03:08:58PM +0300, Martin Ivanov wrote:
>
> Dear R users,
>
> I am using the vars package to calculate the impulse response functions and
> the forecast error variance decomposition of a VAR model. Unfortunately I do
> not know whether these functions assume unit or one s
69 matches
Mail list logo
