Displaying 20 results from an estimated 1000 matches similar to: "list arithmetic"
2001 Sep 25
1
Bug in boxplot.stats?
Is this a bug?
>"xx" <- c(50, 79, 120, 78, 90, 100, 78, 80,
+ 90, 80, 60, 39, 90, 85, 140, 100, 80, 80)
> boxplot.stats(xx)$stats
[1] 60 78 80 90 100
> boxplot.stats(0.1*xx)$stats
[1] 7.8 7.8 8.0 9.0 10.0
I suppose the way the numbers in 0.1*xx are internally stored may
cause this problem:
> formatC(0.1*xx[2],format="f",digits=16)
[1]
2001 Nov 22
1
Bug in dev.print? (PR#1179)
Dear "Debuggers",
please have look a the following:
> plot(0:1,0:1,type="n")
> points(c(0,0.5,1),c(0,0.5,1))
> dev.print(device = postscript, width=5 , height=5, file="test1.ps")
windows
2
># This worked
> dev.print(device = pdf, width=5 , height=5, file="test1.pdf")
windows
2
># That also
>
2001 Oct 05
3
3dim histogram?
Hello all,
I wonder if there is a package including a program to display a 3dim
histogram.
Thanks.
--- D.Trenkler ---
****************************************************************************
*****
Dr. Dietrich Trenkler (dtrenkler at nts6.oec.uni-osnabrueck.de)
Statistik / Empirische Wirtschaftsforschung
Universitaet Osnabrueck
Rolandstrasse 8
2000 Mar 14
2
Bug in sub? (PR#487)
I suspect that there is bug in sub when using "?":
> string_"This is a bug!"
> sub("!", ", or isn't it?", string)
[1] "This is a bug, or isn't it?"
> string_"This is a bug?"
> sub("?", ", or isn't it?", string)
[1] "This is a bug?"
Regards,
*** D.Trenkler ***
2000 Dec 17
1
AW: Permutations
Niels Waller wrote:
> Does anyone know of an R (or S-PLUS) function for delineating all possible
> combinations and permutations?
The following function delivers all permutations of 1:n.
all.perm <- function(n) {
p <- matrix(1, ncol = 1)
for (i in 2:n) {
p <- pp <- cbind(p, i)
v <- c(1:i, 1:(i - 1))
for (j in 2:i) {
v <- v[-1]
2000 Jul 26
1
Bug in stem? (PR#617)
# Your mailer is set to "none" (default on Windows),
# hence we cannot send the bug report directly from R.
# Please copy the bug report (after finishing it) to
# your favorite email program and send it to
#
# r-bugs@biostat.ku.dk
#
######################################################
"b" <- c(-5.5, -4, -2.3, -1.5, -1, -0.5, -0.41, -0.33, -0.29,
-0.26, -0.2,
2000 Dec 13
0
choose(n, k) for k>n: An inconsistency?
It took me by surprise to find that choose(4,5) delivers
[1] NaN
Warning message:
NaNs produced in: choose(n, k)
If we look at choose(4,5) as the number of ways of choosing
5 objects from 4 I would have expected 0 as result.
Furthermore dhyper(5,4,6,5) does deliver 0 and this essentially
equivalent to choose(4,5)*choose(6,0)/choose(10,5).
2000 Apr 07
4
Bug in qbinom? (PR#511)
n_10;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one works as
expected
n_100;p_0.5;jjx_0:n;qbinom(pbinom(jjx,n,p),n,p) # This one causes
severe problems
I cannot interrupt using ESC and I finally have to resort to the Windows
Task manager to kill the R session.
A friend of mine told me that he faced similar problems under Unix.
--please do not edit the information below--
Version:
2000 Mar 07
1
R-crash using cut-and-paste (PR#472)
Dear R-Team,
sorry for bugging (:-)), it's me again. In the FAQ coming along with the
R distribution it reads:
If R executes an illegal instruction, or dies with an operating system
error message that indicates a problem in the program (as opposed to
something like "disk full"), then it is certainly a bug.
This encourages me to describe a crash I encountered a few minutes ago.
I
2000 Mar 06
2
anova-bug in R-version 1.0.0? (PR#470)
# Your mailer is set to "none" (default on Windows),
# hence we cannot send the bug report directly from R.
# Please copy the bug report (after finishing it) to
# your favorite email program and send it to
#
# r-bugs@biostat.ku.dk
#
######################################################
Under R version 0.6.51 the following
A_c(13,9,15,5,25,15,3,9,6,12)
B_c(42,24,41,19,27)
2000 Mar 07
2
AW: anova-bug in R-version 1.0.0? (PR#470)
I think I've discovered what went wrong.
My workspace included a function wilcox.test formerly copied from
ctest. Now ctest ist part of the distribution and because of that I
always got the message:
[Previously saved workspace restored]
Error in autoload("wilcox.test", "ctest") :
Object already exists
I didn't take care of this message. Now I removed
2000 Mar 20
1
Elementwise c() for lists
I have two lists of equal lengths. Each one consists of vectors. Is
there an elegant way to combine them elementwise? Elegant means: without
loops.
For instance, for
> jja_list(1:3,4:7)
> jjb_list(rep(0,4),(2,5))
I would like to get:
> list(c(1:3,rep(0,4)),c(4:7,rep(2,5))).
--- D.Trenkler ---
************************************************************************
*********
2000 May 03
0
Combinatorics for nonparametric tests
For some small sample sizes I would like to exemplify the combinatorics
underlying certain nonparametric tests like Mann-Whitney-Wilcoxon,
Kruskal-Wallis and Spearman's rank correlation. I have written a
function all.perm which delivers all permutations of 1,2,...,n and
which works reasonably well. This can be used to generate P(R=r) of
Spearman's rank correlation:
2000 Apr 06
1
Generating Rd format under Windows
Excuse me if this question has been answered before.
I just came across "Writing R documentation" in the "Writing R
extensions"
manual. Under section 2.9: "Processing Rd format" I only
find hints at generating the Rd format on a Unix platform. Have I
overlooked how to generate them under Windows?
Any help is gratefully appreciated.
--- D.Trenkler ---
2001 Oct 18
1
AW: General Matrix Inverse
Thorsten is right. There is a direct formula for computing the Moore-Penrose
inverse
using the singular value composition of a matrix. This is incorporated in
the following:
mpinv <- function(A, eps = 1e-13) {
s <- svd(A)
e <- s$d
e[e > eps] <- 1/e[e > eps]
return(s$v %*% diag(e) %*% t(s$u))
}
Hope it helps.
Dietrich
2002 May 22
2
Bug in pexp (PR#1590)
I wonder if something like this has been reported
before:
> pexp(85:86,0.438)
[1] 1 -Inf
--please do not edit the information below--
Version:
platform = i386-pc-mingw32
arch = i386
os = mingw32
system = i386, mingw32
status =
major = 1
minor = 5.0
year = 2002
month = 04
day = 29
language = R
Windows NT 4.0 (build 1381) Service Pack 6
Search Path:
.GlobalEnv,
2000 Mar 07
0
AW: R-crash using cut-and-paste (PR#472)
> -----Ursprüngliche Nachricht-----
> Von: Prof Brian D Ripley
> Gesendet am: Tuesday, March 07, 2000 11:30 AM
> An: dtrenkler@nts6.oec.uni-osnabrueck.de
> Cc: r-devel@stat.math.ethz.ch; R-bugs@biostat.ku.dk
> Betreff: Re: [Rd] R-crash using cut-and-paste (PR#472)
>
> On Tue, 7 Mar 2000 dtrenkler@nts6.oec.uni-osnabrueck.de wrote:
>
> > Dear R-Team,
> >
2004 Jul 15
1
GHK simulator
Dear R-community,
not to re-invent the wheel I wonder if someone of you
has ever written a function to compute the GHK smooth recursive
simulator to estimate multivariate normal probabilities. See for instance
page 194 of
@BOOK{Greene97,
author = {William H. Greene},
year = 1997,
title = {Econometric Analysis},
edition = {3rd},
publisher = {Prentice-Hall},
address = {New Jersey
2004 Aug 06
3
Bug in qnorm or pnorm?
I found the following strange behavior using qnorm() and pnorm():
> x<-8.21;x-qnorm(pnorm(x))
[1] 0.0004638484
> x<-8.22;x-qnorm(pnorm(x))
[1] 0.01046385
> x<-8.23;x-qnorm(pnorm(x))
[1] 0.02046385
> x<-8.24;x-qnorm(pnorm(x))
[1] 0.03046385
> x<-8.25;x-qnorm(pnorm(x))
[1] 0.04046385
> x<-8.26;x-qnorm(pnorm(x))
[1] 0.05046385
> x<-8.27;x-qnorm(pnorm(x))
2004 Aug 23
1
Two factor ANOVA with lm()
The following is a data frame
> "jjd" <- structure(list(Observations = c(6.8, 6.6, 5.3, 6.1,
7.5, 7.4, 7.2, 6.5, 7.8, 9.1, 8.8, 9.1), LevelA = structure(c(1,
1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3), .Label = c("A1", "A2",
"A3"), class = "factor"), LevelB = structure(c(1, 1, 2, 2,
1, 1, 2, 2, 1, 1, 2, 2), .Label =