similar to: chisq.test on a matrix

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! 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. ----------------------------

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

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.

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

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

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 |

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

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

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

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