similar to: converting factor to numeric

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

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

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

Hi, Is there a function in R that I can use to fit the data with skew t distribution? Speaking in detail, I first used the kernel density estimation to fit my data, then I drew the skew t using my specified location, scale, shape, and df to make it close to the kernel density. Now I want to get the parameter estimations of the skew t which give me the closet density to the kernel density.

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

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

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

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

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

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

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)

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.

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

Hola! Este mail llegó a r-help hoy, yo lo mandó también a esta lista. Es interesante con competiciones de predicción/clasificación! ¿Alguien que quiere cooperar? Kjetil ---------- Forwarded message ---------- From: Kjetil Halvorsen <kjetilbrinchmannhalvorsen@gmail.com> Date: Fri, Jan 29, 2010 at 12:19 Subject: Fwd: [R] Classification of supernovae - a challenge To: Kjetil Halvorsen

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

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,]