similar to: AW: error-prone feature?

Hi all, once more about the old subj :-) My data has too much various distribution families and for every particular experiment I need just to decide whether the data is "quite homogeneous" or it has two or more clusters. I've revisited the following libraries: amap, clust, cclust, mclust, multiv, normix, survey. And I didn't find any ready-to-use general

Hi, Does anybody know a function to simulate a Markov chain given a probability transition matrix and an initial state ? Thanks. Philippe -- -------------------------------------------------- Philippe Hup? Institut Curie - Equipe Bioinformatique 26, rue d'Ulm - 75005 PARIS France +33 (0)1 42 34 65 29 Philippe.Hupe at curie.fr <mailto:Philippe.Hupe at curie.fr>

Hi All, while looking why the cclust(cclust) doesn't work for 1-dimensional data, I've found unpleasant behavior in semantics of R. Indeed: is.matrix(matrix(cbind(c(1,2,3,4)),ncol=2)[1:2,]) == TRUE but: is.matrix(matrix(c(1,2))[1:2,]) == FALSE kind regards, Valery A.Khamenya --------------------------------------------------------------------------- Bioinformatics

Hello, I would like to write in a variable the following string : Hello "World" If I do cat("Hello ","\"World\""), I have the good result on the screen but it can not be affected to a variable. If I do the same thing with paste paste("Hello ","\"World\"") it does not work since \" seems to be not recognized. So what

Dear R developers, I had some problem when building package: for exemple when building the package e1071 available from CRAN, I get the following message error: * checking for file 'e1071/DESCRIPTION' ... OK * preparing 'e1071': * cleaning src * running cleanup * creating vignettes ... ERROR /usr/lib/R/bin/texi2dvi: pdflatex exited with bad status, quitting.

Hi, I'm trying to debug a function I wrote to calculate the cindex for a hierarchical tree. For this it is useful to compare my calculations with those in output from the clustindex function, in the cclust library. There's no way, however, to have the cindex value for a given output of the cclust function, as a NA value is always returned. This happens almost surely because the cindex in

Hello, my data is an ordered list of observations where each observation is either of condition "efficient" (coded as 1) or "non-efficient" (0) as e.g. No. Obs Condition 1 1.1 0 2 1.2 1 3 1.4 1 4 1.5 0 5 2.5 1 etc. My goal is to plot a single line with lty='o' but different symbols for either condition whenever condition is 1 or 0, i.e. on the x-axis the

I am trying to use the function aggregate with the median function but I get the following error: Error in FUN(X[[1]], ...) : Argument "INDEX" When I replace median by mean, it works perfectly Can someone tell me where the problem comes from? Thx I am running R 2.0.0 on SunOS 5.9 -- Philippe Hup?? UMR 144 - Service Bioinformatique Institut Curie Laboratoire de Transfert (4??me

Hello, here is my problem, After looking at the mail archives, I found a description of the error I get when I use this package. At first I even tought that they were showing how to solve it. But the thing is that by saying "the programmer forgot drop=FALSE" doesn't show me how I should get rid of the problem I have looked inside the package very quickly and I found three

Hello, I have a MANOVA model and I want to test the following hypothesis: LBM = 0 where B is the parameter estimates. Is there any function to do this in R? Cheers, Philippe -- Philippe Hup? Institut Curie, CNRS UMR 144, INSERM U900 26 rue d'Ulm 75005 Paris - France Email : Philippe.Hupe at curie.fr T?l : +33 (0)1 56 24 69 91 Fax: +33 (0)1 56 24 69 11 website :

Hi all, Why subscripting a one column matrix drops one dimension? > x<- matrix(rnorm(100), ncol=1) > str(x) num [1:100, 1] -0.413 -0.845 -1.625 -1.393 0.507 ... > str(x[20:30,]) num [1:11] -0.315 -0.693 -0.771 0.448 0.204 ... > str(x[20:30]) num [1:11] -0.315 -0.693 -0.771 0.448 0.204 ... This breaks: > cov(x) [,1] [1,] 0.9600812 >

Hi- I'd like to apply some of the cluter indices, using 'clustindex' in the cclust package, but using output from hierarchical clustering using 'hclust' in the mva package. The two packages use a different format for storing the results of clustering; has anyone coerced one into the other? Thanks! -John Barnett

I have to solve a clustering problem. My first step is to determinate the number of clusters, that's why I 'm using the Calinski index ( [tr(b)/(k-1)]/[tr(w)/(k-1)] ) which i try to maximize to have the best number of clusters. A function is already implemented in R to calculate this index : clustIndex(cl,x, index="calinski") where cl is the result of a clustering method ,

Hi all, is there any package in R capable of smooth approximation of CDF basing on given sample? (Thus, I am not speaking about ecdf) In particular, I expect very much that the approximation should subject to the property: f(x0)<=f(x1) for x0<x1, where x0 and x1 belong to range of the sample given. Polynomial approximation could be OK for me as well. P.S.

Hi All is there "homals" package prepared for win32? kind regards, Valery A.Khamenya --------------------------------------------------------------------------- Bioinformatics Department BioVisioN AG, Hannover

> On only ten points, what did you expect ? Even with 1000 > observations, estimating a density is difficult, and has > been the subject of a century of research. Kernel density > estimates are among the most successful. For your immediate > application, try plot(density(rnorm(10)), type="l"), etc. wait, you misunderstood me! I'd like to see 10 or 9 points with

Hi, I just installed Ubuntu 12.04.4. I installed R by following instructions from CRAN after modification into my /etc/apt/sources.list file sudo apt-get update sudo apt-get install r-base sudo apt-get install r-base-dev When I try to load library(gdata) ou library(gplots), I have error message : library(gdata) gdata: read.xls support for 'XLS' (Excel 97-2004) files

> Nothing to do with me: you should report problems with > packages to the > maintainers, rather than R-help or a member of R-core. OK. I've sent a note about cclust patch to Evgenia Dimitriadou Thank you for your valueable comments. (No more reply needed in this thread) kind regards, Valery A.Khamenya ---------------------------------------------------------------------------

Hi All, following expression: x <- sort(rnorm(10)); e <- ecdf(x); d <- numericDeriv(e(x),"x"); makes d far from approximation of one dimensional pdf. What's wrong then here? Kind regards. --------------------------------------------------------------------------- Valery A.Khamenya Bioinformatics Department BioVisioN AG, Hannover

Hi All, We'd like to use functions provided in R in our application. Our application is written in C/C++ and currently runs on win32, Linux and Mac. We'd be happy to attach the whole R ( i.e. not just transfer some function by hand). It is important that we deal with big amount of data, so "command line"-like invocations won't be very interesting. We'd