similar to: substitution in expression

Dear all, I am having trouble trying to install R-1.8.1 on a Sun Solaris 8 (Generic_108528-23 version) machine. The configuration was successful but make check fails. I traced the the problem to the plot() function. > 1 + 1 [1] 2 > capabilities() jpeg png tcltk X11 GNOME libz http/ftp sockets TRUE FALSE TRUE TRUE FALSE TRUE TRUE TRUE

I am trying to write a wrapper for the last example in help(cv.glm) that deals with leave-one-out-cross-validation (LOOCV) for a logistic model. This wrapper will be used as part of a bigger program. Here is my wrapper funtion : logistic.LOOCV.err <- function( formu=NULL, data=NULL ){ cost.fn <- function(cl, pred) mean( abs(cl-pred) > 0.5 ) glmfit <- glm(

Hi ! I would like to understand where do affinity.spline.coefs used in function compute.affinities come from ? library(gcrma) data(affinity.spline.coefs) affinity.spline.coefs X1 X2 X3 X4 X5 X1 -0.55004171 -0.58579091 -0.08870557 -0.47774242 0.23205570 0.58002746 X2 X3 X4 X5 X1 X2

Could some kindly tell me if I am supposed to be getting the same test statistic value with var.equal=TRUE and var.equal=FALSE in t.test ? set.seed(1066) x1 <- rnorm(50) x2 <- rnorm(50) t.test(x1, x2, var.equal=FALSE)$statistic # 0.5989774 t.test(x1, x2, var.equal=TRUE)$statistic # 0.5989774 ??? Here are my own calculations that shows that perhaps the result when var.equal=TRUE is

Hello guys, I am wondering the default way of transferring arguments in R. Is it by value or by ref in default case, or could that be changed explicitly? Cheers, Xiaofan --- Xiaofan Li Department of Applied Mathematics and Theoretical Physics University of Cambridge

Hi, I want to do automatic reading of a number of tables (files) stored in ascii format without having to specify the variable name in R each time. Below is an example of how I would like to use it (I assume files pair1,...,pair8 exist in spec. dire.) for (i in 1:8){ name <- paste("pair",i,sep="") ? ? ? <-

Hi, I have a list of numbers. For each of the numbers, I take sum of squares of the numbers centered on the number chosen. If it is less than a certain constant, I will take the average of the numbers chosen. Anyone can give me a sample code. You help will be greatly appreciated. Peter [[alternative HTML version deleted]]

There was a BioConductor thread today where the poster wanted to find pairwise difference between columns of a matrix. I suggested the slow solution below, hoping that someone might suggest a faster and/or more elegant solution, but no other response. I tried unsuccessfully with the apply() family. Searching the mailing list was not very fruitful either. The closest I got to was a cryptic chunk

Hi. I'd like to do a t-test to compare the Delta values of items with Crit=1 with Delta values of items with Crit=0. What is the t.test syntax? It should produce a result like this below (I can't get in touch with the person who originally did this for me) Welch Two Sample t-test data: t1$Delta by Crit t = -3.4105, df = 8.674, p-value = 0.008173 alternative hypothesis: true

Is there a method to save a large and complex R object (either as a binary or text file) so that it can be loaded and reused at a later time? Specifically, I am creating large lists (several thousand elements), each element of which is either a vector or a matrix (with ~ 2000 rows). The dimensions of the matrices are not all the same. My ideal would be a set of functions of the form obj

Hi! Here a simple example. setClass("myclass" ,representation(info="character") ,contains="matrix" ) rownames(dd)<-c("a","b") tt<-new("myclass",dd) #the source of pain. as(tt,"matrix")<-matrix(1,3,3) Error: length of dimnames [1] not equal to array extent Is there a different way to do what I would like to do (I

Hi! Here a simple example. setClass("myclass" ,representation(info="character") ,contains="matrix" ) rownames(dd)<-c("a","b") tt<-new("myclass",dd) #the source of pain. as(tt,"matrix")<-matrix(1,3,3) Error: length of dimnames [1] not equal to array extent Is there a different way to do what I would like to do (I

Dear all, Is there any simple way in R that can I put the all the interactions of the variables in the OLS model? e.g. I have a bunch of variables, x1,x2,.... x20... I expect then to have interaction (e.g. x1*x2, x3*x4*x5... ) with some combinations(2 way or higher dimensions). Is there any way that I can write the model simpler? Thanks! Leaf

I have a vector with about 100.000 values representing a quite regular function (sinusoid like). I would like to find all local maxima of this function (should be about 4000). Is there a native routine for R? Thanks in advance Eryk. -- _|_ \|/ \|/ Eryk Witold Wolski tel :0049-(0)30-8413-1543 w w ?v? 'v? \'v'/ MPI Moleculare Genetik fax :0049-(0)30-8413-1139 |

This seems like a bug to me. Can someone verify this? First we define a function f that returns its second argument and lapply it to 1:2 using 9 as the second argument and all seems well. Note that "{" as a function does the same thing as f, as illustrated with f(1,9) and "{"(1,9); however, when we attempt to use "{" in the very same way we used f in lapply,

Hi, Executing the following code causes the prompt to enter browser() mode. >dd<-as(Npatrick1,"list") Called from: as(Npatrick1, "list") Browse[1]> R : Copyright 2004, The R Foundation for Statistical Computing Version 2.0.0 beta (2004-09-20), ISBN 3-900051-07-0 Some addtitional info about Npatrick1 > is(Npatrick1) [1] "Massvectorlist"

Hi, I clustered a distance matrix and would like to draw it using plot.hclust or plot.dendrogram. The dendrogram is not informative because I have a few extremely small dissimilarities in the distance matrix (e.g. 0), but most of the other distances are in the range 1e10+-5000. I would like to show the tree only for the height of 1e10+-5000 but unfortunately their are no parameter like

Hi, Is the ordering of the values returned something on what I can rely on, a form of a standard, that a function called unique in R (in futher versions) will return the uniq elements in order of they first occurcence. > x<-c(2,2,1,2) > unique(x) [1] 2 1 Its seems not to be the standard. E.g. matlab >> x=[2,2,1,2] x = 2 2 1 2 >> unique(x) ans =

Hi, hres <- hclust(smatr,method="single") hresd<-as.dendrogram(hres) as.dendrogram(hres) `dendrogram' with 2 branches and 380 members total, at height 2514.513 plot(hresd,leaflab="none") #<-error here. #the plotted dendrogram is incomplete. The x axis is not drawn. #The interested reader can download the save(hresd,file="hres.rda") #from the

Hi, Please consider the following trivial example. In the sample code I first define a class MVCa using composition. This class contains a slot of type list. Next I define an inheritance relation using the function setIs -> MVCa "is a" list. Then I define a function printlength which prints the list-length which is defined for the type "list". If I call this function