similar to: A basic question

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

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

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.

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

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

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

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

May be you aren't on the right working dir: getwd(),setwd(). Your R-script are able to find the thyphoon.txt path? 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: alessandro.semeria at

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

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

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.

Question: How can a repeat a dataframe n times in the direction of the rows? (Or: How I can rbind a dataframe n times ?) Example Data: x <- data.frame( alpha=letters[1:3], num=1:3 ) n <- 4 Complicated solution: To rbind a dataframe n times I can program a loop: xr <- x for( i in 2:n ) xr <- rbind( xr, n ) Observation: To cbind a dataframe n times I

Dear all, I want to use GraspeR (on R 1.6.0) and have to install tcltk library. When I try to do this, I get back this error message: Error in firstlib(which.lib.loc, package) : TCL_LIBRARY is not set Error in library(pkg, character.only = TRUE) : .First.lib failed But I have a tcl/tk package (Active TCL 8.4.1.0) correctly install on my PC (Win95) and in the autoexec.bat the

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