similar to: binom.test bug?

I suspect that the recent queries about reading database files into R implicitely answered my question here, but I'll ask anyway. Are there any native binary file read routines in R? I need to a binary data file consisting of 2 byte records, strip off the upper 6 bits (using a logical AND, leaving a ten bit word), and then write the ten bit number in text form for R. Assuming that R

> I want to combine (rbind) many dataframes into a single data frame, but "automatically" > specifying the names of the dataframes as listing of object names. > E.g., combine these 18 df objects into one big df using something conceptually like this : > rbind(objects(pattern="*.df")) Brian Ripley suggested that something along the lines of:

I am having trouble understanding what must be a quite simple rule of scope. In the code example below, the argument to the main function, my.subset, would like to be passed to the function LME. This does not happen however. Curiously, the other argument (my.data) is passed to LME. This problem goes away with this: assign("my.subset", my.subset, env = sys.frame()). It works, but it

I would like a function that return the X,Y position of the mouse cursor whenever it is called PERIOD, i.e., without the button presses of locator() or identify(). E.G., > where.is.mouse() [1] 123 34 # X and Y coordinates of present mouse position in current graphics window Is this a conceivable R function or does it need to be addressed with a call to the OS (Windows 2000 in my case)?

How is exact pattern matching achieved in GREP (and GREPlike) functions ? # Want: listing of all object names that end in *.lm > objects(pattern="*.lm",pos=1) # ... but get: all objects that partially match *.lm, e.g., *.lme [1] "j3.lm" "J3.lme" "j8.lm" "J8.lme" # Want: position of string "4jan2002" in vector >

Windows-32 has a time structure called FILETIME, a 64-bit value representing the number of 100-nanosecond intervals since January 1, 1601 (UTC). That is not a typo, the year is 1601. Does anyone have a clue(or algorhithm)for how this is converted to something a little more POSIX-like ? Thank you, Derek -- Derek N. Eder Gothenburg University VINKLA - Vigilance and Neurocognition

I want to migrate my R workspace databases (e.g., .Rdata) from Windows to Linux. The R help file for "save" says that "[a]ll R platforms use the XDR representation of binary objects in binary save-d files, and these are portable across all R platforms." Is it is simple as that? Copy the Windows .Rdata file to the Linux .ext3 file system? Thanks, Derek Derek N. Eder

My question is so basic that I am (almost too) embarrassed to admit that I could not find an answer after an hour's worth of homework. What is the Trellis / Lattice analog for the axis(graphics) function that enables the creation of axes in locations other than the default (i.e., bottom for X axis and right for Y axis) ? For example when plotting mileage against weight (in American

Isn't there a cryptic footnote somewhere in the blue book that describes regular expressions as "a funny little language ..." ? Well, I'm not laughing! I've googled and Man paged, but the GREP-style solution to the following problem still eludes me: ... Using a regular expression in the "pattern" argument to list.files() to restrict the return to filenames which

I have a time stamp in UTC (GMT) time: > format(ISOdatetime(1970,1,1,0,0,0)+1165398135729/1000,"%Y-%m-%d %H:%M:%OS3") "2006-12-06 09:42:18.823" (note millisecond accuracy, but not relevant to question here) Now, this time stamp actually "happened" at local (Swedish) time one hour later (10:42). Regarding summer/winter adjustments in time ("spring

Has anyone out there compiled R for the Umbutu Linux* (ne? Debian) v. 5.04 distribution for Intel-type platforms (32 and 64 bit) ? Thank you, Derek Eder * Umbutu, a popular new Linux distribution, not a Nigerian scam, I promise! http://www.ubuntulinux.org/ -- Derek Eder SDS KLINIKEN Vasaplatsen 8 SE 411 34 G?teborg (Gothenburg) Sweden phone: +46 (31)* - 10 77 80 fax: +46 (31)* - 10 77

Full_Name: Petr Savicky Version: 2.4.0 OS: Fedora Core release 2 Submission from: (NULL) (62.24.91.47) the error is > binom.test(0.56*10000,10000) Error in binom.test(0.56 * 10000, 10000) : 'x' must be nonnegative and integer while > binom.test(5600,10000) yields correct result. The same error occurrs for > binom.test(0.57*10000,10000)

Hello all, I am trying to understand the different results I am getting from the following 3 commands: chisq.test(c(62,50), p = c(0.512,1-0.512), correct = F) # p-value = 0.3788 binom.test(x=62,n=112, p= 0.512) # p-value = 0.3961 2*(1-pbinom(62,112, .512)) # p-value = 0.329 Well, the binom.test was supposed to be "exact" and give the same results as the pbinom, while the chisq.test

Hi all, I am wondering whether there is a small bug in the binom.test function of the ctest library (I'm using R 1.6.0 on windows 2000, but Splus 2000 seems to have the same behaviour). Or perhaps I've misunderstood something. the command binom.test(11,100,p=0.1) and binom.test(9,100,p=0.1) give different p-values (see below). As 9 and 11 are equidistant from 10, the mean of the

Hello everybody. Does anyone else find the last test in the following sequence odd? Can anyone else reproduce it or is it just me? > binom.test(100,200,0.13)$p.value [1] 2.357325e-36 > binom.test(100,200,0.013)$p.value [1] 6.146546e-131 > binom.test(100,200,0.0013)$p.value [1] 1.973702e-230 > binom.test(100,200,0.00013)$p.value [1] 0.9743334 (R 1.5.1, Linux RedHat 7.1) --

No, since I'd like to test null: p <= p0 alternative: p > p0. and my understanding is that binom.test tests null: p = p0 (can only be a "simple" null hypothesis according to help(binom.test)) alternative: p > p0 (or p < p0 or p != p0). Thanks, Mirko. > -----Urspr?ngliche Nachricht----- > Von: Douglas Bates [mailto:bates at stat.wisc.edu] >

Hi Folks, The recent correspondence about "strange fisher.test result", and especially Peter Dalgaard's reply on Tue 03 April 2007 (which I want to investigate further) led me to take a close look at the code for binom.test(). I now have a query! The code for the two-sided case computes the p-value as follows: if (p == 0) (x == 0) else if (p == 1) (x == n)

R-experts: A quick question, please. >From a lab exp, I got 12 positives out of 50. To get 90% CI for this , I think binom.test might be the one to be used. Is there a better way or function to calculate this? > binom.test(x=12, n=50, p=12/50, conf.level = 0.90) Exact binomial test data: 12 and 50 number of successes = 12, number of trials = 50, p-value = 1 alternative

R 0630 for windows > library(ctest) > binom.test(7,10,p=0.3, alternative="two.sided") returns a p-value of =< 2.2e-016 and a warning In Splus 3.4 > binom.test(7,10,p=0.3, alternative="two.sided") returns a p-value of 0.0106 I think it is the max(v[v<=(1+eps)*PVAL]) causing the problem... max() of an empty vector....... Mai Z

running R 1.4.1 on MAC and 1.2.2 on Linux When I use run binom.test with small N the results are a little perplexing to me >binom.test(9,20,p=0.5) gives the below plus other stuff 95 percent confidence interval: 0.2305779 0.6847219 Now: >pbiom(9,20,0.6847219) [1] 0.02499998 # i.e., lower 2.5% of distribution >pbinom(9,20,0.2305779) [1] 0.9923132 >pbinom(8,20,0.2305779)