similar to: KS test in R.1.3.0 has incorrect p-values. (PR#1004)

I'm using the dump command to pass data to WinBUGS/OpenBUGS/JAGS and have run into a problem. Here is some R-code: foo <- array(1:6, dim=c(2,3)) dump('foo', file='dumpdata.R') dump('foo', file='dumpdata.R', append=TRUE, evaluate=TRUE) foo2 <- array(c(2,3,5,7,9,7,5,3), dim=c(2,4)) dump('foo2', file='dumpdata.R', append=TRUE) And here is

Hi All, I just used "dump" on Splus to transfer a pile of survival objects from Splus 3.4 on Solaris 7 to R 0.63.3 on Intel. The only trick is that survival objects contain an element holding the original call that generated the object. When Splus writes these out, it doesn't mark them in any way, so when R tries to read them in, it ends up trying to reevaluate the call. Not

Hi All, I just used "dump" on Splus to transfer a pile of survival objects from Splus 3.4 on Solaris 7 to R 0.63.3 on Intel. The only trick is that survival objects contain an element holding the original call that generated the object. When Splus writes these out, it doesn't mark them in any way, so when R tries to read them in, it ends up trying to reevaluate the call. Not

Hi, Sorry if this is already known... contour() sometimes draws contour lines that intersect. Is there a temporary fix? A dataset which causes problems is at http://www.maths.uwa.edu.au/~adrian/dumpdata.R If you try just source("dumpdata.R") image(huh) contour(huh) the 100 x 100 matrix 'huh' contains an hourglass-shaped region of values around 0.8. The contour plot

Greetings all, I am calculating two t-test values for each of many files then save it to file calculate another set and append, repeat. But I can't figure out how to write it to file and then append subsequent t-tests. (maybe too tired ;} ) I have tried to use "dump" and "file.append" to no avial. ttest_results = tempfile() two_sample_ttest <- t.test (tempA, tempB,

The note to 1004 says "fixed for 1.3.1" Uh. No. It ain't. The problem was more serious than guessed as even the simplest testing would show. For example, Example 5.4 in Hollander and Wolfe (Nonparametric Statistical, Methods, 2nd ed., Wiley, 1999, pp. 180-181) R Version 1.3.1 (SuSE Linux 7.1) > X <-

> tmp <- data.frame(y=rnorm(12), a=factor(rep(letters[1:4],3))) > tmp y a 1 -0.60866099 a 2 0.55500538 b 3 0.12231693 c 4 -0.24613790 d 5 -0.09253593 a 6 -1.54652581 b 7 0.17204210 c 8 -1.22778942 d 9 1.22151194 a 10 -0.43982577 b 11 -1.25444287 c 12 -0.97251060 d > tmp.aov <- aov(y ~ a, data=tmp) > summary(tmp.aov) Df Sum Sq Mean Sq F value

I've been trying to get plot.hclust to work with a hclust object I created and have not had much success. It seems that there is some "hidden" characteristic of a hclust object that I can't see. This is most easily seen in the following example, where plot.hclust works on one object, but when this object is "dumped" and then re-read, plot.hclust no longer works. Is

On Sun, 1 Jul 2001 mcdowella@mcdowella.demon.co.uk wrote: > Full_Name: Andrew Grant McDowell > Version: R 1.1.1 (but source in 1.3.0 looks fishy as well) > OS: Windows 2K Professional (Consumer) > Submission from: (NULL) (194.222.243.209) Please upgrade: we've found a number of Win2k bugs and worked around them since then, let alone teh bug fixes and improvements in R .... >

Full_Name: Martin Kober Version: 2.6.0 OS: Vista & Linux Submission from: (NULL) (137.208.185.169) Hi! I just stumbled upon a problem with file source()ing: R will fail to source a file if it contains lines longer than about 8192 bytes ("input buffer overflow"). While it's save to say that human-written code won't contain lines that long, dump()ed data structures may

In message <Pine.GSO.4.31.0107010731110.7616-100000@auk.stats>, Prof Brian D Ripley <ripley@stats.ox.ac.uk> writes > >You do realize that the Kolmogorov tests (and the Kolmogorov-Smirnov >extension) assume continuous distributions, so the distribution theory >is not valid in this case? > >S-PLUS does stop you doing this: > >> ks.gof(o,

for a silly question, wondering how to test fit with the one sample as follow. I have read _fitting distributions with R_, but that doesn't answer my specific question. inclined to use Kolmogorov-Smirnov D, and its associative p value. much appreciation! X20.001 232 93 84 185 336 417 228 199 2110 1411 612 1913 1314 3015

Hello! I need to compare 2 datasets whether they come from the same distribution. I use function ks.boot{Matching}. And what is the confidence level of the p-value, returned by ks.boot function? The code is: set=read.table("http://stella.sai.msu.ru:8080/~gala/data/testsets.csv", header=T,sep=',') set1=set[!is.na(set$set1),'set1']

Not an R question directly, but has anyone got a method for placing a moderately large number of (near) equi-spaced points on a sphere? I have a nasty feeling platonic solids are needed for exact solutions and I'm thinking of samplings involving around 200 - 1000 regularly-spaced points, Thanks, Richard Rowe Richard Rowe Senior Lecturer Department of Zoology and Tropical Ecology, James

I think this was asked recently but I can't find the item. Is it possible to find the colour at a point on a graphics page? I want to estimate the area of a complex shape by colouring it then random sampling from an enclosing box ... Richard Rowe Richard Rowe Senior Lecturer Department of Zoology and Tropical Ecology, James Cook University Townsville, Queensland 4811, Australia fax (61)7

I am attempting to plot frequency of occurrence against (calendar) time. As the output is to be produced as a Powerpoint picture on a dark background I wish to make the plot output yellow > par(col="yellow", col.axis="yellow", col.lab="yellow", col.main="yellow") > plot(t2, t1, ylim=c(0,40), main="Episynlestes run 2",

I am working with transition matrices of sequences of animal behaviours. What I would like to do is parse the original matrices, adjusting row/column order so that the matrix has its main values in blocks surrounding the diagonal. This would cause behaviours involved in functional groupings (e.g. grooming, resting, foraging etc) to appear as blocks. This can be done manually by applying

I am having problems with the 'if' syntax. I have an n x 4 matrix, X say. The first two columns hold x, y values and I am attempting to fill the second two columns with the quadrant in which the datapoint (x, y) is and with the heading angle. So I have two problems 1) how to do this elegantly (at which I've failed as I can't seem to vectorize the problem) and 2) how to

I'm trying to use the biological female and male signs in R2.5.1 under Windows XP. I can access and insert these symbols using word-processors. In general these should be available as \u2640 and \u2642 but I can't make them happen in R using (say) text(5,5, "\u2640") message "invalid \uxxxx sequence" I've replaced Arial with Arial unicode MS in Rdevga and

This message is in MIME format --_=XFMail.1.3.p0.Linux:990422185944:4782=_ Content-Type: text/plain; charset=us-ascii R-0.64.0 on RedHat Linux 5.1 Some problems with multivariate fast fourier transform. I have attached the dump of a 30 x 3 matrix that seems to reliably reproduce these problems 1) mvfft doesn't like vectors. It complains about not having enough memory: R>