similar to: vectorizing a matrix computation

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

I have a 20x3 matrix as follows: > m <- replicate(3, matrix(rnorm(20),20,1)) I need to compute, say, 95th and 99th percentiles of each column such that the resulting matrix becomes 2x3 with each row representing the respective percentile. My "best effort" is to compute one column at a time as follows: > quantile(m[,1], c(0.95, 0.99)) To do the same for columns 2 and 3, I

Hi! I am installing R in non standard directory. ./configure --prefix=/non/standard/directory The installation works fine. But after starting R i get .libPaths() /usr/lib/R/library but they cant be there of course! Have I to set some additional switches during config? Eryk Dipl. bio-chem. Eryk Witold Wolski @ MPI-Moleculare Genetic Ihnestrasse 63-73 14195 Berlin 'v' tel:

Hi! >library(cluster) In this lib a function called silhoutte.default is defined than on the R prompt it type >silhouette.default Error: Object "silhouette.default" not found R1.9.0 The same error are at R1.8.1 And I knew that a function silhoutte.default are present. But the same piece of code works in R1.6.2 ??? Eryk Dipl. bio-chem. Eryk Witold Wolski @

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! I guess that it cant work. but maybee I am wrong. I would like to define a function rbind. > rbind function (..., deparse.level = 1) { which works only for objects of class Myclass Is it possible to use something like setMethod("rbind",signature(...="Myclass",deparse.level="numeric") This gives an error. Or should I use rbind.Myclass Eryk. Dipl. bio-chem.

Hi! Would like to plot two image plots with color scale (like levelplot (lattice) implements) in one graphic device. Tried to plot 2 levelplot's (lattice) (the data has not much in common) in one graphics device by setting par(mfrwo=c(2,1)) which does not work. Is there a way to force levelplot to use only a half of the graphic device? Is there a function like image which provides the color

Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with acrobat as document figure the quality is bad. Anyone knows a way to reduce the size but keep the quality? /E -- Dipl. bio-chem. Witold Eryk Wolski MPI-Moleculare

I have been struggling with this problem for a while and I hope someone could help me. Or if someone could point me to a section in the manual I would be grateful. x <- "my" plot(1:10, main=expression(paste( x, Delta, "values" ))) Q : How do I get the title to say "my (triangle symbol) values" ? The following trial-and-error produced mainly errors :

Hallo! I would like that the internal package functions does not appear in the documentation. I placed they names as aliases in the file nameofmypackage-internal.Rd as advised. After R CMD check (which runs neatly) they are still listed in 00Index.html. Is there no way to hide internal functions from the doc? Have I overlook something? Sincerely Eryk. Dipl. bio-chem. Eryk Witold Wolski

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

Hello All: This function obviously fails x <- function(z) paste("go", z, sep = ".") <- 10 x("now") But is there a way to define the name of a variable through passing a parameter in a function call? Thanks, ANDREW

my major is bioinformatics, and i'm trying to cluster ( agglomerate the closest pari of observations ) in R. i have already got my own similarities metric, but do not know how to clust it based on similarities instead of dissimilarities. since the help document of hierclust mentions the parameter "sim", which seems good to me, but it doesn't appear in the code of hierclust()

hello I am returning to some libraries that I had previously customised by amongst other things adding additional functions. I had simply typed these new functions into the file in library/R/thelibrary. These however do not seem to be loaded now as they previously were under older versions of R. What's changed and what do i need to do? Thanks Stephen

Hi, Would like to use to the hexbin package to plot a 2D hist - because it looks really _cool_. My problem is related to drawing a pdf hexbin graphic in series of other graphics. All other figures have a par(mar=c(3.2,3.2,1,1)). So the gplot.hexbin figure in this series looks a little alienated? Was trying a to specify the _mar_ using par, viewport, hexViewport etc. a little. My question. How

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 =

Let me present to you my problem : I have a character vector x and I would like to obtain the indices of the elements of this vector that yielded exactly a match. For example, x=nom, pattern="b", I would to obtain 2 because "b" is on the second position. First program : nom <- c("a","b","ab") grep("b",nom) 2 3 Then I try the

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

Hi! I encountered some things trying to build own packages under Windows 2000 for which i found no answer in the documentation. First. 1. I do not found the rcmd.exe in the bin directory of R_HOME. I installed R using the setup.exe 2. I updated the R intatllation with one of the zip files. (rw1041sp.zip) now i had the rcmd.exe in the bin. 3. I do not found any hint in the docu that perl is