similar to: (no subject)

Good morning ! Today I found a strange, for my poor knowledge of R, behaviour of 'which' on a matrix: HAL9000> str(cluster.matrix) num [1:227, 1:6300] 2 2 2 2 2 2 2 2 2 2 ... HAL9000> class(cluster.matrix) [1] "matrix" HAL9000> ase <- cluster.matrix[1:5,1:5] HAL9000> ase [,1] [,2] [,3] [,4] [,5] [1,] 2 2 2 0 -2 [2,] 2 2 2 0 -2 [3,]

I don't know if there is some fault in compiling or a bug of the new R-1.7.0 version: cl.pam.2 <- pam(as.dist(1-cor(mel.data)),2) plot(cl.pam.2) perform a right partitioning and silhouette plot on the old R-1.6.2 instead "Error in clusplot.default(x$diss,...... ; x is not numeric" is the output on the new R-1.7.0. Same platform: RH8.0 i686. Some suggestions? A.S.

Hi dear R-list. I would superimpose many lines of class 'table' over a plot(table.0,type ="l",...), I tried 'lines(table.n,type="l",..)' but its behaviour is quite different from those of 'plot' : >plot(table,type="l"...) # start graph from x = 0 >lines(table,type="l"...) # start graph from x = 1 Some sugestions? Thanks

Hi! Probably I perform this question because I did'nt still understand the R-philosophy. I have to build many matrix, with different dimensions, and I would to assign them a same 'prefix' name, i.e. "aval", but, obviuosly different suffix, something like: for (i in 1:n) { aval%i% <- matrix(scan(data....),nrow=nr[i],ncol=nc[i] ... } where "%i%"

I would write in tab-text file a table like this one: TAB colname1 TAB colname2 TAB TAB... colnameN rowname1 # # # rowname2 # # # rownameM # # # then I wrote something like: mymatrix <- matrix(nrow=M,ncol=N,byrow=T) rownames(mymatrix) <- chvector1 colnames(mymatrix) <- chvector2

Hello! What kind of calculation perform 'chisq.test', if its argument is a matrix? Thanks for any suggestion! ------------------------------------------------------------ Sincerely yours. Dr. Alessandro Semeria Tel. +39 544 536811 Models and Simulation Lab of Fax. +39 544 538663 The Environment Research Center - Montecatini

Look at 'lm' in the 'base' package for linear fit, in the manual "An Introduction to R", distributed with R, page 55.... Best! A. S. ---------------------------- Sincerely yours. |------------------------------------+------------------------------------| |Alessandro Semeria |Tel. +39 544 536811 |

RE-Hello dear R-list! I have to test if N (non linear) experimantal distributions fit some of known one. I would proceed in this way: 1. Fit with 'nlm' and e.g. power law one ( non outlier) distribution of mine. 2. Perform 'chisq.test' for each pairs (power law ,distribution(i)) but I'm neither a statistician nor a skilled on R, so I ask you: ther'is some other way, more

This is not the right place for this question, but however: when you connect to a remote machine with ssh with a user name "user" you aren't able to use the specific shell profile ($PATH,....) of "user", then you have to specify the complete PATH of all command (i.e. /usr/local/bin/R). Best A.S. ---------------------------- Alessandro Semeria Models and Simulations

On gentoo: emerge R and dowload,configuration, compilation of all you need to working with base R is performed. Best! A.S. ---------------------------- Alessandro Semeria Models and Simulations Laboratory Montecatini Environmental Research Center (Edison Group), Via Ciro Menotti 48, 48023 Marina di Ravenna (RA), Italy Tel. +39 544 536811 Fax. +39 544 538663 E-mail:

Hello! I would plot more than 10 (this number is always different) graphics on a figure with layout like 'mfrow=c(4,2)'. There is some existing function able to open the right number of device if number of graphics is done or I have to build it myself with some 'IF.....ELSE.....'? Thanks. A.S. ----------------------------

I have written a function to plot data which will be used for various different chemistries. A simplified version is: plot_data <- function(risk,levels,chem,sd2,measure){ plot(risk, levels,main=paste ("per", sd2, measure, "\n in usual", chem)) } The problem is with the title. This works fine if the variable "chem" is just text, but if it is an expression then

Dear all, Is there any generic algorithm code for optimization implemented in R? I searched without success. Thanks, -- Zhu Wang Statistical Science Department Southern Methodist University Phone: (214)768-2453 Fax: (214)768-4035 Email: zhuw at mail.smu.edu

History files saved using nFile Save History do not open when using File Load History. I notice that when they are saved they do not seem to have a file extension and there is no option to chose an extension. Is this a problem. Other than this it seems so simple that I can't see what the problem is. I've tried using loadhistory and savehistory with the same problems. The syntax guide

Hello, I am analysing several samples whose sizes are from 9 to 110. I would like to test their distribution with R, whether they are normal or not. I wonder which test for normality from R should I use . Thank you for help. Samuel Samuel BERTRAND Doctorant Laboratoire de Biomecanique LBM - ENSAM - CNRS UMR 8005 151, bd de l'Hopital 75013 PARIS Tel. +33 (0) 1 44 24 64 53 Fax +33 (0) 1

Hello! I don't know if there is some R-package able to associate ontology to a long list of GeneBank Name (a txt-tab file or an XML file), i.e. I would as output a formatted file with 4 columns (1:GeneBank Name 2,3,4:ontology). I know that I have to perform a mapping of genes, I got a look on AnnBuilder pkg, but I 've not idea from where to start. Some suggestion? Thanks in advance!

If some can help to find under wich package can i find the commando's wich i schould to use to work with the empirical cummulative (ecdf) distribution and also the QQplot. Zmarrou Hicham Univesity of Amsterdam -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

I can''t start r inside of emacs. When I try it says that it cannot find .Rhistory and then it says " This program ( pressumably Rterm) has performed an illegal operation and will be shut down.". How can I fix this problem? [[alternate HTML version deleted]]

everytime i try to read in data i get the following error: Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : line 3 did not have 5 elements it was my understanding it wasn''t necessary to specifiy all defaults. The bizarre thing is this error comes from using read.table. When I use scan I get the following error: Error in

Hi all, After a short hiatus away from R I have found that it's changed a bit. I used to keep a definition of .First in .Rprofile that did a couple of things on startup (load a couple of libraries). Now, I've discovered that when I change the definition of .First in .Rprofile it doesn't change anything when I start up, because .First is held over in .Rdata from the last session.