Response to 1You need the log version e.g. in maximum likelihood, otherwise the
product of the densities and probabilities can become very small.
Ursprüngliche Nachricht Von: r-help-requ...@r-project.org
Datum: 04.08.21 12:01 (GMT+01:00) An: r-help@r-project.org Betreff: R-help
.
Best wishes,
Matthias
PS. Sorry for the HTML email. I’ve given up trying to fix such behavior.
Von: Martin Maechler
Gesendet: Dienstag, 13. Juli 2021 09:09
An: Matthias Gondan
Cc: r-help@r-project.org
Betreff: Re: [R] density with weights missing values
>>>>> Matthias Gonda
er
Datum: 12.07.21 18:44 (GMT+01:00) An:
r-help@r-project.org, matthias-gondan , Bert Gunter
Cc: r-help@r-project.org Betreff: Re: [R] density with
weights missing values Sure, you might think that.But most likely the reason
this code has not been corrected is that when you give weights for mi
thing?
Ursprüngliche Nachricht Von: Bert Gunter
Datum: 12.07.21 16:25 (GMT+01:00) An: Matthias
Gondan Cc: r-help@r-project.org Betreff: Re: [R]
density with weights missing values The behavior is as documented
AFAICS.na.rmlogical; if TRUE, missing values are removed from x. If FALSE
anymissing v
a.rm=TRUE)
[1] 2
Von: Richard O'Keefe
Gesendet: Montag, 12. Juli 2021 13:18
An: Matthias Gondan
Betreff: Re: [R] density with weights missing values
Does your copy of R say that the weights must add up to 1?
?density doesn't say that in mine. But it does check.
On Mon, 12 Jul 2021
Dear R users,
This works as expected:
• plot(density(c(1,2, 3, 4, 5, NA), na.rm=TRUE))
This raises an error
• plot(density(c(1,2, 3, 4, 5, NA), na.rm=TRUE, weights=c(1, 1, 1, 1, 1, 1)))
• plot(density(c(1,2, 3, 4, 5, NA), na.rm=TRUE, weights=c(1, 1, 1, 1, 1, NA)))
This seems to work (it trigge
Dear R friends,
I am currently trying to write a piece of C code that uses „embedded R“, and
for specific reasons*, I cannot keep track if R already has been initialized.
So the code snippet looks like this:
LibExtern char *R_TempDir;
if(R_TempDir == NULL)
…throw exception R not initialize
happen? [it’s documented
behavior, but still…]
So, I’ll stick with R. Still 25 years or so until retirement, but I’ll survive,
even without crossreferenced default arguments.
Best wishes,
Matthias
Von: S Ellison
Gesendet: Dienstag, 5. September 2017 16:17
An: Matthias Gondan; r-help@r-project.
mu=0.53, sigma2=4.3^2) # instead of l=u
And maybe also „in-place“ changes of values…
Best regards,
Matthias
Von: William Dunlap
Gesendet: Samstag, 2. September 2017 19:41
An: Rui Barradas
Cc: Matthias Gondan; r-help@r-project.org
Betreff: Re: [R] Strange lazy evaluation of default argume
Dear R developers,
sessionInfo() below
Please have a look at the following two versions of the same function:
1. Intended behavior:
> Su1 = function(u=100, l=u, mu=0.53, sigma2=4.3^2)
+ {
+ print(c(u, l, mu)) # here, l is set to u’s value
+ u = u/sqrt(sigma2)
+ l = l/sqrt(sigma2)
+ mu =
I found it:
quantile(ordered(1:10), probs=0.5, type=1)
works, because type=1 seems to round up or down, whatever. The default option
for is 7, which wants to interpolate, and then produces the error.
Two options come to my mind:
- The error message could be improved.
- The default type could
Dear R users,
This works:
quantile(1:10, probs=0.5)
This fails (obviously):
quantile(factor(1:10), probs=0.5)
But why do quantiles for ordered factors not work either?
quantile(ordered(1:10), probs=0.5)
Is it because interpolation (see the optional type argument) is not defined? Is
there a
The warning
1: In (ind.c == TRUE) & (ind.sgn == TRUE) :
longer object length is not a multiple of shorter object length
means that ind.c and ind.sgn have different lengths, for whatever reason.
Although R continues the routine, the warning should, in general, not be
ignored.
Try:
1:3 + 1:2
lim=c(-Inf, Inf)[as a replacement for NULL/
autoselection] and ylim=c(Inf, -Inf)[autoselection, reversed y axis] should be
handled correctly.
I would find such a feature useful. Do you think it would interfere with other
functions?
Thank you for your consideration.
Best wishes,
Matth
you might to do something like
library(SuppDists)
t = runif(100, 100, 500) # original RT
t_IG = qinvGauss(ecdf(t)(t)-0.5/length(t), 1, 16)
plot(density(t_IG))
but what is the purpose of it? Usually reaction times are thought to
follow a certain kind of distribution (e.g. an inverse Gaussian
dis
My vote:
1. Symbolic function arguments:
fn = function(a, b)
{
a/b
}
fn(b=10, a=2)
2. Names for elements of a vector and matrices
v = c(a=1, b=2)
v['a'] = v['a'] * 2
same for matrices
3. about 10,000 user-contributed packages on CRAN
4. weird things like
a = numeric(10)
a[1:10] = 1:2
interaction.plot source:
legend(xleg, yleg, legend = ylabs, col = col, pch = if (type %in%
c("p", "b"))
pch, lty = if (type %in% c("l", "b"))
lty, bty = leg.bty, bg = leg.bg) <- here I woul
Dear R team, dear Prof. Therneau,
library(survival)
data(colon)
?colon
gives me only a very rudimentary source (only a name). Is there a
possibility to get a reference to the clinical trial these data
are taken from?
Many thanks in advance. With best wishes,
Matthias Gondan
2011 20:50, schrieb Ben Bolker:
Matthias Gondan gmx.de> writes:
Hi,
Use offset variables if count occurrences of an event and you want to
model the
observation time.
glm(count ~ predictors + offset(log(observation_time)), family=poisson)
If you want to compare durations, look at library(surv
Hi,
Use offset variables if count occurrences of an event and you want to
model the
observation time.
glm(count ~ predictors + offset(log(observation_time)), family=poisson)
If you want to compare durations, look at library(survival), ?coxph
If tnoise_sqrt is the square root of tourist noise,
Dear R developers,
I want to draw an arrow in a figure with lty=2. The
lty argument also affects the edge of the arrow, which is
unfortunate. Feature or bug?
Is there some way to draw an arrow with intact edge, still
with lty=2?
Example code:
plot(1:10)
arrows(4, 5, 6, 7, lty=2)
Best wishes,
Dear R people,
On my freshly installed R-2.13.1 (winxp), the following code yields
unsatisfactory results (irregular grid lines instead of a smooth plane):
image(matrix(1, nrow=100, ncol=100))
This is working fine,
image(matrix(1, nrow=100, ncol=100), useRaster=TRUE)
but then right-click and
yrnorm(4))
[1] 0.2875775 0.7883051 0.4089769 0.8830174
Best wishes,
Matthias Gondan
--
GMX DSL Doppel-Flat ab 19,99 Euro/mtl.! Jetzt mit
gratis Handy-Flat! http://portal.gmx.net/de/go/dsl
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
1" or the equivalent for your
platform?
Allan
On 06/08/10 15:12, Matthias Gondan wrote:
Dear List,
I am aware this is slightly off-topic, but I am sure there are people
who already had the problem and who perhaps solved it.
I am running long-lasting model fits using constrOptim command.
try this (package Rsolnp)
library(Rsolnp)
g<- function(x)
{
return(x[1]^2+x[2]^2)
} # constraint
f<- function(x)
{
return(x[1]+x[2])
} # objective function
x0 = c(1, 1)
solnp(x0, fun=f, eqfun=g, eqB=c(1))
Am 10.08.2010 14:59, schrieb Gildas Mazo:
Thanks, but I still cannot get to sol
try command solnp in package Rsolnp
Am 09.08.2010 18:56, schrieb Dwayne Blind:
Hi !
Why not constrOptim ?
Dwayne
2010/8/9 Gildas Mazo
Dear R users,
I'm looking for tools to perform optimization subject to constraints,
both linear and non-linear. I don't mind which algorithm may be used, m
be helpful.
Searching around the internet was not very encourageing. Some people
wrote that it is not so simple to have Atlas fully exploit a multicore
computer.
I hope this is not too unspecific.
Best wishes,
Matthias
--
Dr. rer. nat. Matthias Gondan
Institut für Psychologie
Universität Regensbu
are some R scripts from Brunner et al. available on his website
http://www.ams.med.uni-goettingen.de/de/sof/ld/index.html
But they seem not to be working with current R versions.
Best regards,
Matthias Gondan
--
Sicherer, schneller und einfacher. Die aktuellen Internet-Browser -
jetzt koste
Dear R developers,
A great R feature is that elements of vectors, lists and dataframes can
have names:
vx = c(a=1, b=2)
lx = list(a=1, b=2)
Accessing element "a" of vx: vx['a']
Accessing element "a" of lx: lx[['a']] or lx$a
Might be a matter of taste, but I like the $ very much. Unfortunatel
. rer. nat. Matthias Gondan
Institut für Psychologie
Universität Regensburg
D-93050 Regensburg
Tel. 0941-943-3856
Fax 0941-943-3233
Email: matthias.gon...@psychologie.uni-regensburg.de
http://www.psychologie.uni-r.de/Greenlee/team/gondan/gondan.html
_
Dear R-Users,
Is anyone aware of a significance test which allows
demonstrating that one distribution dominates another?
Let F(t) and G(t) be two distribution functions, the
alternative hypothesis would be something like:
F(t) >= G(t), for all t
null hypothesis: F(t) < G(t), for some t.
Best
Dear list,
Apparently, there is no function like sscanf in R.
I have a string, "Condition: 311", and I would like
to read out the number and store it to a numeric
variable. Is there an easy way to do this?
Best wishes,
Matthias
--
__
R-help@r-project
Dear list,
I hope the topic is of sufficient interest, because it is not
R-related. I have N=100 yes/no-responses from a psychophysics
paradigm (say Y Yes and 100-Y No-Responses). I want to see
whether these yes-no-responses are in line with a model
predicting a certain amount p of yes-responses.
Dear list,
I think there is a small bug in the constrOptim documentation (R-2.7.0):
## from optim
fr <- function(x) { ## Rosenbrock Banana function
x1 <- x[1]
x2 <- x[2]
100 * (x2 - x1 * x1)^2 + (1 - x1)^2
}
grr <- function(x) { ## Gradient of 'fr'
x1 <- x[1]
x2 <- x[2]
Hi Geoff,
I think the answer to such a problem (overall survival vs. disease free
survival) depends on the regulatory
environment, for example, in a phase III clinical trial, OS would be
used, whereas in an equivalence study,
DFS would be used.
Best,
Matthias
Geoff Russell schrieb:
> Dear Us
Frank E Harrell Jr schrieb:
> Matthias Gondan wrote:
>>> data(colon)
>>> s = survfit(Surv(time, status) ~ rx, data=colon)
>>> plot(s)
>>> plot(s, col=1:3)
>>
>> By the way: Does anyone know a neat way to indicate the number of
>> patients un
Hi Eleni, hi list,
here is small sample program, the library is "survival", it is included
in the standard R distribution.
> data(colon)
> s = survfit(Surv(time, status) ~ rx, data=colon)
> plot(s)
> plot(s, col=1:3)
By the way: Does anyone know a neat way to indicate the number of
patients unde
on the Cohen book.
>
>
> On Feb 13, 2008 4:52 AM, Matthias Gondan <[EMAIL PROTECTED]> wrote:
>
>> Dear list,
>>
>> Is anyone aware of a library for sample size calculation in R, similar
>> to NQuery? I have to give a course in this area, and I would
Hi Eleni,
The problem of this approach is easily explained: Under the Null
hypothesis, the P values
of a significance test are random variables, uniformly distributed in
the interval [0, 1]. It
is easily seen that the lowest of these P values is not any 'better'
than the highest of the
P values
Dear list,
Is anyone aware of a library for sample size calculation in R, similar
to NQuery? I have to give a course in this area, and I would like to
enable the students playing around with this.
Best wishes,
Matthias
__
R-help@r-project.org mailing
(little correction below)
> Dear list,
>
> This is off-topic, but perhaps there is an expert out there who can help
> me in
> a bootstrapping test.
>
> I have two samples A, B from, say, a normal distribution. A third sample R
> is the vector of pairwise minima of the same random variables:
>
>
Dear list,
This is off-topic, but perhaps there is an expert out there who can help
me in
a bootstrapping test.
I have two samples A, B from, say, a normal distribution. A third sample R
is the vector of pairwise minima of the same random variables:
A = rnorm(100)
B = rnorm(100)
A2 = rnorm(10
Hi Diana,
Look at Edgar Brunner's Homepage
http://www.ams.med.uni-goettingen.de/en/sta/e.brunner.html
and here:
http://www.ams.med.uni-goettingen.de/en/sof/index.html
Unfortunately, the R script for a 2x2 design with two groups and two
timepoints (f1.ld.f1.r) seems to be broken. Uncommenting t
Dear List,
I have tried a stratified Cox Regression, it is working fine, except for
the "Anova"-Tests:
Here the commands (should work out of the box):
library(survival)
d = colon[colon$etype==2, ]
m = coxph(Surv(time, status) ~ strata(sex) + rx, data=d)
summary(m)
# Printout ok
anova(m, test='Ch
Peter Dalgaard schrieb:
> Matthias Gondan wrote:
>
>> Dear R users,
>>
>> I noticed a problem in the anova command when applied on
>> a single coxph object if there are missing observations in
>> the data:
>>
...
>> In the docume
Dear R users,
I noticed a problem in the anova command when applied on
a single coxph object if there are missing observations in
the data:
This example code was run on R-2.6.1:
> library(survival)
> data(colon)
> colondeath = colon[colon$etype==2, ]
> m = coxph(Surv(time, status) ~ rx + sex
ALARY" [label="133"];
"TIME" -> "PUBS" [label="38.97" dir=both];
"TIME" -> "TIME" [label="18.3" dir=both];
"PUBS" -> "PUBS" [label="196.12" dir=both];
"SALARY" -> "SA
47 matches
Mail list logo
