search for: brinchmann

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. --

I encountered something strange with read.spss (package foreign, version 0.7 with R2.0.0 and version 0.8 with R2.0.1 beta, windows XP) I made a test file test.sav with SPSS version 11.5.1 containing only one numeric variable, with a value label for one value not occuring in the file. According to ?read.spss this should result in a factor, but it results in all NA. Using the argument

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 am now downloading r-devel to compile it on windows XP. The CRAN source code page says "you will prefere to use rsync". I am googling around, and cannot find anything about rsync on windows. Anybody has any experience with rsync on windows? Kjetil

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

I was surprised by: > test <- matrix( as.character(1:4), 2) > test [,1] [,2] [1,] "1" "3" [2,] "2" "4" > apply(test, 1, paste, sep="+") [,1] [,2] [1,] "1" "2" [2,] "3" "4" > apply(test, 1, paste, sep="*") [,1] [,2] [1,] "1" "2" [2,]

I get the following > library(foreign) Error in namespaceExport(ns, exports) : undefined exports: write.foreign Error in library(foreign) : package/namespace load failed for 'foreign' with rw2000dev as of (2004-09-17 Kjetil -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi Elmandjra

I am trying to install a local package and get this unexpected error: ---------- Making package UMSA ------------ adding build stamp to DESCRIPTION installing R files installing data files installing man source files installing indices Error: couldn't find function "na.omit" Execution halted na.omit of course is in package stats, and that is listed in the Depends field in

I just discovered that when the VGAM package (not on CRAN) is loaded, glm() doesn't work. This is because VGAM defines a family function() which gets found by glm() in place of the family function from stats. Then VGAM:::family returns an object which doesn't have a $family component, (it has a component $vfamily). I thought namespaces should protect us from this happening? Kjetil --

Hola! I am writing a package with some large data files, so it would be good to use the feature of zipping data. However, I use .R files to read in the data from one (sometimes more) files, so Rcmd build --zip-data gives a package which does'nt work. How should I do it, is it better ti zip each data file individually and use the functions to read from zipped files, or are there some other

With rw2010dev I get a strange protect(): protection stack overflow error with a small data frame which otherwise is usable: If anybody wants to have a look I can provide an RData file with the problematic data frame. Doesn't seem to be necessary, the following simulated example generates the error: > testmat <- matrix(1:80, 20,4) > dim(testmat) [1] 20 4 > str(testmat) int

...but get the same error: Error in par(.split.screens[[cur.screen]]) : invalid value specified for graphics parameter "fig". The help usefully says that they are defined in NDC units, but I don't know what an NDC unit is, and there isn't any example. The code in kjetil brinchmann halvorsen's message on R-help on Mar 31 2003 (do a search for "NDC units"!) didn't work either, it gives the same message: > split.screen( matrix( c(0, 0.3, 0.5, 1, 0.3, 0.7, 0.5, 1, + + 0.7, 1, 0.5, 1, 0, 0.5, 0, 0.5, + + 0.5, 1, 0, 0.5), 5, 4, byrow=TRUE)) Error in par(.s...

Using the heatmap function in mva, it seems to be hard to use different colors in the edges leading to different groups of objects, as commonly done in many heatmaps in the microarray graphics. Any suggestions? Thanks. max

...ilk). Thanks, Tamas [1] You could ask why I don't just estimate the trans. prob. matrix just once from the whole dataset. The problem is that the process is meant to model (real) interest rates with only three levels, and involves many peculiar assumptions. By the way, many thanks to Kjetil Brinchmann Halvorsen, Rolf Turner, Ted Harding, Patrick Burns and Martin Maechler for giving detailed and helpful replies to my earlier question about Markov chains in R. -- Tam?s K. Papp E-mail: tpapp at axelero.hu (preferred, especially for large messages) tpapp at westel900.net Please try to send...

Hola! The R FAQ says: 7.12 How do I convert factors to numeric? It may happen that when reading numeric data into R (usually, when reading in a file), they come in as factors. If f is such a factor object, you can use as.numeric(as.character(f)) to get the numbers back. More efficient, but harder to remember, is as.numeric(levels(f))[as.integer(f)] In any case, do not call as.numeric()

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(*,

Hello! I was reading R-intro and I have some suggestions: R-intro.html#A-sample-session rm(fm, fm1, lrf, x, dummy) suggestion rm(fm, fm1, lrf, x, y, w, dummy) The next section will look at data from the classical experiment of Michaelson and Morley to measure the speed of light. file.show("morley.tab") mm <- read.table("morley.tab") suggestion mm <- data(morley)

Suppose I have x<-data.frame(v1=1:4, v2=c(2,4,NA,7), v3=rep(1,4), v4=LETTERS[1:4],v5=rep('Z',4)) or a much larger frame, and I wish to test for and remove the constant numeric columns. I made: is.constant<-function(x){identical(min(x),max(x))} and apply(x,2,is.constant) # Works for numerics x[,-which(apply(x,2,is.constant))] I'd really like to be able to

Dear List, We are planning a genotyping study to be analyzed using false discovery rates (FDRs) (See Storey and Tibshirani PNAS 2003; 100:9440-5). I am interested in learning if there is any consensus as to how many features (ie. how many P values) need to be studied before reasonably reliable FDRs can be derived. Does anyone know of a citation where this is discussed? Bill Dupont William D.

Hi there, Is there an implementation of the Gaveh Stehfest algorithm in R somewhere ? Or some other inversion ? Thanks, Tolga