similar to: S-Plus graphs to LaTeX picture commands

This is a message for whoever maintains "chisq.test": For an outcome more extreme than 2000 simulations, a Monte Carlo p-value of "< 2.2e-16" was printed. Ripley said the proper p-value for such cases should be 1/(B+1) = 1/2001. This can be easily fixed by adding "if(PVAL==0)PVAL <- 1/(B+1)" right after the following line in the code for chisq.test (in R

Hola! I have the following in a .Rd file: \eqn{\mbox{coef} = c(\mbox{coef}[1],\ldots, \mbox{coef}[n]) } {coef = c(coef[1], coef[2], \dots, coef[n])} However, both arguments come out in the latex file! Whats happening? Kjetil Halvorsen

Hola! I ask here since I learnt from this list that the LaTeX package listings should be good for typesetting R code. I encountered one problem: \begin{lstlisting} X %*% V \end{lstlisting} in the output the * in %*% disappears! same with %/%, etc, the / disappears. Any ideas? Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. --

This is R1.9.0, windows XP I run across the following: It is admittedly silly, but it shouldn't segfault: test <- list(y=1:100, x=rep(NaN, 100)) approx(test, xout=rep(NaN, 10)) Bombs! Kjetil Halvorsen

The NEWS file for R-1.7.0 says: table() now allows exclude= with factor arguments (requested by Michael Friendly). But in a recent R-1.7.1beta: > test <- c(1,2,3,4,5,1,2,3,4,5,NA,NA,1,2,3,4,5) > table(test) test 1 2 3 4 5 3 3 3 3 3 > table(test, exclude=NULL) test 1 2 3 4 5 <NA> 3 3 3 3 3 2 > table( as.factor(test), exclude=NULL) 1 2

I just found a bug in fisher.test(). This is rw1080, on windows XP. A p-value can certainly not be Inf, but: > religion Costumbres rel orig Religion Si Algunas veces Nunca católica 2121 4700 6234 prot/evan 100 216 2461 otra C 27 67 502 otra 0 0 14 > fisher.test(religion, workspace=2000000)

Hola! Starting with rw1071, just after starting Rgui the scope is not with Rconsole, but some other place. This means I have to do a mouse click in the Rconsole before starting to write the first command. Since I always forget this, I end up writing the first command twice. This is of course only a minor nuisance, but why is focus no longer set to Rconsole window when start-up? This is on

Hola! I have some problems with gam in mgcv. Firts a detail: it would be nice igf gam would accept an na.action argument, but that not the main point. I want to have a smooth term for time over a year, the same pattern repeating in succesive years. It would be natural then to impose the condition s(0)=s(12). Is this possible within mgcv? I tried to obtain this with trigonometric terms, aca:

I should'nt have sent the last mail so fast. Same problem with \eqn{u_j = a_j + b x + c x^2, \quad j=1, \ldots, r-1} {u[j] = a[j] + b*x + c*x^2 j = 1,\dots,r-1} I thought the problem in the first case could have to do with the use use of \mbox{} (with the braces) within the arguments of \eqn, but here there are none braces in the arguments of \eqn{}{}. Another case giving the

Hola! I am using package strucchange, and encounters the following: > bp <- breakpoints(diesel90 ~ regress -1, h=NULL) Error in La.chol2inv(x, size) : element (14, 14) is zero, so the inverse cannot be computed The obvious problems have been checked, that is, the model matrix is of full rank. What can be causing this? I can send some data if that can be of help. Kjetil Halvorsen

Hola! I will want to identify pixels in an image with the mouse, for so getting the image data from the matrix(es), for use in subsequent discriminant analysis. But the following bombs R: (windows XP, rw1071) > str(baboon) list() - attr(*, "size")= int [1:2] 512 512 - attr(*, "cellres")= num [1:2] 1 1 - attr(*, "bbox")= num [1:4] 0 0 512 512 - attr(*,

I saw a reort on R-help (without answer) about not being able to load old workspaces (.RData files) in 1.9.0alpha. I have seen the same problem with some of my workspace, although most of them loads without problems. (rw1090alpha as of 17.03, windows XP). Why will some workspaces not load? Kjetil Halvorsen

I am working with the precompiled rw1090beta as of 2004-03-23. ?ts has in section "See also" "Standard package 'ts' for many additional time-series functions." td does not longer exist. ?plot.ts example section uses library(ts) ?lines.ts example section has library(ts) Kjetil Halvorsen

Hola! Any explanation of the following: > summary(as.Date(rep(NA,20))) Min. 1st Qu. Median Mean 3rd Qu. Max. NA NA NA NA NA NA NA's "1970-01-21" Where does this 1970 date come from? Kjetil Halvorsen

Hola! I got problems using an objects returned from arima (in KalmanSmooth(my.ts, ModArima$model), because my.ts showed up to have storage mode "integer" (is.integer(my.ts was TRUE). Should storage.mode() of a ts be allowed to be integer, should ts() someplace say storage.mode(ts.out) <- "double", or maybe inside arima() storage.mode(x) <- "double"

Hola! I just downloaded and installed the factory-fresh rw1081. Windows XP. In my Rprofile file in the \etc subdirectory I put (among others) options(continue=" ") options(width=65) options(scipen=5) options(show.signif.stars=FALSE) R is honoring all of this on startup, exept the width=65 option, which does'nt have any effect. I tried other values for width also (70, 90) ,

Hello friends, I would like to thank you all for your advices. Rcmdr is gread. All the best, Mehdi ********************* Mehdi Zarrei, PhDPostdoctoral fellowUniversity of Toronto Dept. of Ecology & Evolutionary Biology Royal Ontario Museum Department of Natural History Green Plant Herbarium (TRT) 100 Queen's Park Toronto, On., M5S 2C6, CanadaE-mail:

Have I been sleeping in class? rw1071 from CRAN, windows XP incidencia is made by a call to tapply > class(incidencia) [1] "array" > incidencia <- unclass(incidencia) > class(incidencia) [1] "array" Kjetil Halvorsen

Hola! I have data in the form of a symmetric distance matrix, in the file I have recorded only the upper triangular part, with diagonal. The matrix is 21x21, and the file have row and col names, and some other information. I am trying to read with the following code (I tried many variations on it, but all give the same error). The items in the data file is delimited by white space. (Part

Because of Plain text mode, my output example is split across several lines and hence is totally meaningless. If you want to see the example, follow this link : http://download.caoba.fr/r3tb7hxb2pl5/ I marked in red the places where the bug appears. ------------------------- Santiago DIEZ Quark Systems & CAOBA 23 rue du Buisson Saint-Louis, 75010 Paris ------------------------- On Thu, Feb