Re: [R] bwplot and panel.bwplot

2004-01-04 Thread Patrick Giraudoux
I'm not aware of any such package, but I'll assume you are talking about something else. Yes sure, it was libray(lattice), and not (trellis) indeed... I wrote the code line in the email directly after two many hours work... a bit confused. Sorry. Difficult to guess without access to the

RE: [R] ade4

2004-01-04 Thread Marwan Khawaja
Dear Stephane, Yes, it worked fine -- a nice plot -- many thanks for your help. Marwan --- Marwan Khawaja http://departments.aub.edu.lb/~mk36 ---

[R] Analyzing dendograms??

2004-01-04 Thread Johan Lindberg
I have used heatmap to visualize my microarray data. I have a matrix of M-values. I do the following. #The distance between the columns. sampdist - dist(t(matrix[,]), method=euclidean) sclus - hclust(sampdist, method=average) #The distance between the rows. genedist - dist(matrix[,],

[R] array problem

2004-01-04 Thread Z P
Dear all, I define , for n=5 or any integer greater than 0. A-array((1/2)^n , c(rep(2,n))) then for any i not equal to j, and 1=i,j=n, B-apply(a,c(i,j),sum) now B is a 2 by 2 matrix, I also define another costant 2 by 2 matrix G, How can I change the values of each elements of array A,

RE: [R] array problem

2004-01-04 Thread Gabor Grothendieck
Replace your line that updates A with this: p - unique(c(i,j,1:5)) f - function(x) diag( matrix(x,4,4) ) AA - apply( outer(A,G/B), p[-(1:2)], f ) # form product A - aperm( array( AA, dim(A) ), order(p) ) # reshape Here are a couple of tests. You might want to do some more tests yourself as

Re: [R] Analyzing dendograms??

2004-01-04 Thread Thomas W Blackwell
Johan - Disclaimer: I've never used heatmap(), so probably I shouldn't be answering this. However ... the function heatmap() probably calls either plot() or image() (regular graphics) or xyplot() (lattice graphics) in order to set up axes and initialize the actual plotting. heatmap()

[R] Error running as.phylo (package ape)

2004-01-04 Thread Rogério Rosa da Silva
Hi all, I try to convert hclust object to phylo object, using as.phylo (package ape), but got an error like the following: tese- read.table (tese_guildas.txt, header=T) library(mva) library(ape) hclust.tree-hclust(dist(tese[1:156,])) phylo.tree- as.phylo(hclust.tree) Segmentation fault

[R] Density Plots

2004-01-04 Thread David Tyler
I am using an older version of R (1.6.2) to run a Monte Carlo simulation, generating 10,000 samples per 'run'. When I plot histograms I get the expected 'bins' on the x-axis and the frequency distribution on the y-axis. However when I ask R to plot the SAME data set with a density curve the

Re: [R] Density Plots

2004-01-04 Thread Jason Turner
David Tyler wrote (using an e-mail client that doesn't wrap lines): I am using an older version of R (1.6.2) to run a Monte Carlo simulation, generating 10,000 samples per 'run'. When I plot histograms I get the expected 'bins' on the x-axis and the frequency distribution on the y-axis.